Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 618,040 A→T S1228S (TCA→TCT)  SAUSA300_RS02920 → serine‑aspartate repeat‑containing protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793618,0400AT82.6% 85.4 / 40.8 132S1228S (TCA→TCT) SAUSA300_RS02920serine‑aspartate repeat‑containing protein D
Reads supporting (aligned to +/- strand):  ref base A (22/1);  new base T (104/5);  total (126/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-02

Only 100 of 132 total aligned reads displayed.

TCAGATAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTGACTCAGA‑TTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGACTCAGATAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGATAGTGATTCAGACTCAGA  >  NC_007793/617954‑618138
                                                                                        |                                                                                                  
ccAGATAGCGATTCAGATTCAGACAGCGCCTCAGACGCAGATAGCGATTCAGGTGCGGACAGCGATTCAGA‑CGCAGATAGCGATT‑CTGATGTGGACAgat                                                                                       >  1:3579915/2‑98 (MQ=0)
 cAGATAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTGACTCAGA‑TTCAGATAGCGACTGCAGAtcggaagagca                                                                                       >  1:1542661/1‑90 (MQ=14)
      aGCGACTCTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAAACGCCGATAGCGACTCTGC‑TGCGGATAACGACT‑CAGATGCAGACAGCGAttacg                                                                                 >  1:3130173/1‑97 (MQ=2)
                  tcggacagcgactcagattcagacagcgactcagacgcaGATAGCGACTCAGA‑TGCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCg                                                                     >  1:3829705/4‑100 (MQ=11)
                  gcggacagcgactcagacgcagatagcgactcagacgcaGATAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCg                                                                     >  1:5955169/4‑100 (MQ=11)
                    ctctagcgactctgatgcggacagcgactcagacgcagaTAGCGACTCAGA‑CGCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGat                                                                   >  1:5097565/5‑99 (MQ=2)
                           gactcagacgcagatagcgactcagacgcagatagcGACTCTGA‑TGCGGACAGCGACT‑CAGATTCCGACAGCGATTCAGACGCAGATAGCGACtctgatt                                                            >  1:6015302/1‑100 (MQ=14)
                             ctcagattcagacagcgactcagacgcagaTAGCGACTCAGA‑CGCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGATtctgatgcg                                                          >  1:3731227/1‑97 (MQ=11)
                              cctgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGattcagattcag                                                         >  1:124014/10‑100 (MQ=11)
                               ttgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcaga                                                        >  1:6222744/9‑100 (MQ=11)
                               gcgaggcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcaga                                                        >  1:3230065/9‑100 (MQ=11)
                               ctgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACtcagacgcaga                                                        >  1:4408402/9‑100 (MQ=2)
                               ctgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACtcagacgcaga                                                        >  1:5304036/9‑100 (MQ=2)
                               ctgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACtcagacgcaga                                                        >  1:3143221/9‑100 (MQ=2)
                               ctgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACtcagacgcaga                                                        >  1:3848184/9‑100 (MQ=2)
                               cctgggcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACtcagacgcaga                                                        >  1:382393/9‑100 (MQ=2)
                               cagactcagacagcgactcagactcagaTAGTGACTCAGA‑TTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGACTCAGATAGCGACtcagattcaga                                                        >  1:2792668/1‑100 (MQ=255)
                               cagactcagacagcgactcagactcagaTAGTGACTCAGA‑TTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGACTCAGATAGCGACtcagattcaga                                                        >  1:127730/1‑100 (MQ=255)
                               agcttgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATtctgattcgga                                                        >  1:3149353/9‑100 (MQ=2)
                               acgaagcagatagcgactctgatgcggaCAGTGACTCAGA‑CGCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAAtctgaagcaga                                                        >  1:6652682/7‑100 (MQ=2)
                                tgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCCGACGCAGATAGCGAttcagattcagac                                                       >  1:758795/8‑99 (MQ=2)
                                tgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATtccgattcagac                                                       >  1:5102163/8‑99 (MQ=11)
                                tgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACagcgactcagacgcagatagcgactcagacgcagat                                                       >  1:4610633/8‑100 (MQ=11)
                                tgatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGCCtcagacgcagat                                                       >  1:5727756/8‑100 (MQ=2)
                                ctgtgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACACCGATTCAGACGCAGATAGCGATtcagattcagat                                                       >  1:410638/8‑100 (MQ=11)
                                ccctgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACagcgactcagacgcagatagcgactcagacgcagat                                                       >  1:5421114/8‑100 (MQ=11)
                                agatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACagcgactcagacgcagatagcgactcagacgcagat                                                       >  1:3520794/8‑100 (MQ=11)
                                 gctgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:414679/7‑98 (MQ=11)
                                 gattcagacagcgactcagacgcagaTAGCGACTCAGA‑CGCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGATtctgatgcggaca                                                      <  1:2223079/97‑8 (MQ=0)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATGGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:6224015/7‑98 (MQ=2)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCCGATAGCGAttcagattcagaca                                                      >  1:2239772/7‑98 (MQ=2)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGCttcagattcagaca                                                      >  1:4885022/7‑98 (MQ=2)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:1176726/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:4442970/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:3683718/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:5030851/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:6466110/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:6339632/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:2368281/7‑98 (MQ=11)
                                 gatgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGAttcagattcagaca                                                      >  1:2175714/7‑98 (MQ=11)
                                 gatgcagacagcgactcagacgcagaTAGCGATTCAGA‑TGCGGACAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGAAAGCGATtctgatgcagaca                                                      >  1:3086470/5‑98 (MQ=2)
                                  gtgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGcgactcagacgcagatagcgactcagacgcagatag                                                     >  1:1795222/6‑100 (MQ=11)
                                  ctgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGcgactcagacgcagatagcgactcagacgcagatag                                                     >  1:4910679/6‑100 (MQ=11)
                                  ctgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGcgactcagacgcagatagcgactcagacgcagatag                                                     >  1:2142470/6‑100 (MQ=11)
                                  atgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAg                                                     >  1:506086/6‑100 (MQ=11)
                                  atgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAg                                                     >  1:2369391/6‑100 (MQ=11)
                                  atgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGcgactcagacgcagatagcgactcagacgcagatag                                                     >  1:3028127/6‑100 (MQ=11)
                                   tgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGc                                                    >  1:4969799/5‑100 (MQ=11)
                                   tgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCgactcagacgcagatagcgactcagacgcagatagc                                                    >  1:2185964/5‑100 (MQ=11)
                                   tgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCgactcagacgcagatagcgactcagacgcagatagc                                                    >  1:2980657/5‑100 (MQ=11)
                                   cgcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCGGATAGCGATTCAGATTCAGACAGc                                                    >  1:4183365/5‑100 (MQ=11)
                                   cacggatagcgactctgatgcggaCAGCGACTCAGA‑CGCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCAGACAGc                                                    >  1:336843/8‑100 (MQ=2)
                                    gcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCg                                                   >  1:5619826/4‑100 (MQ=11)
                                    gcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCg                                                   >  1:6530453/4‑100 (MQ=11)
                                    gcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGACTCAGATTCAGACAGCg                                                   >  1:6234034/4‑100 (MQ=11)
                                    gcggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCAATTCAGACGCAGATAGCGATTCAGATTCAGACAGCg                                                   >  1:4113938/4‑100 (MQ=11)
                                    acggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGactcagacgcagatagcgactcagacgcagatagcg                                                   >  1:4243821/4‑100 (MQ=11)
                                     gggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGAGGACAGCGATTCAGACGAAGATAGCGATTCAGATTCAGACAGCGa                                                  >  1:5580588/3‑100 (MQ=2)
                                     cggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGa                                                  >  1:904622/3‑100 (MQ=11)
                                     cggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGActcagacgcagatagcgactcagacgcagatagcga                                                  >  1:3716803/3‑100 (MQ=11)
                                     cagatagcgactcagattcagaTAGTGATTCAGACTTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGGCTCAGATAGTGACTCAGATTCAGATAGCg                                                   >  1:3159701/1‑100 (MQ=11)
                                     aggccagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGActcagacgcagatagcgactcagacgcagatagcga                                                  >  1:1814697/5‑100 (MQ=11)
                                      taacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAAACGCAGATAGCGac                                                 >  1:4817571/3‑99 (MQ=11)
                                      ggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGac                                                 >  1:178558/2‑99 (MQ=11)
                                      ggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGac                                                 >  1:632321/2‑99 (MQ=11)
                                      ggacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGac                                                 >  1:6036896/2‑99 (MQ=11)
                                      ccacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCCGACAGCGATTCAGACGCAGATAGCGATTCCGATTCGGTTAACGac                                                 >  1:4060829/3‑99 (MQ=2)
                                      agatagcgactcagacgcagaTAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGAt                                                 >  1:3313709/1‑100 (MQ=11)
                                      agacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGAt                                                 >  1:4109666/1‑100 (MQ=2)
                                      agacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGac                                                 >  1:2039525/1‑99 (MQ=11)
                                      agacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGACTCAGACGCAGATAGCGACTCAGACGCCGATAGCgcc                                                 >  1:5823165/1‑98 (MQ=11)
                                      agacagcgactcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGActcagacgcagatagcgactcagacgcagatagcgac                                                 >  1:4481661/1‑99 (MQ=11)
                                       gacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATacagacagatcgg                                                >  1:2665582/1‑95 (MQ=2)
                                       gacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGcactc                                                <  1:6291527/100‑5 (MQ=11)
                                       gacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGact                                                >  1:5017869/1‑98 (MQ=11)
                                       gacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGacc                                                <  1:4769031/100‑3 (MQ=11)
                                       gacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGActcagacgcagatagcgactcagacgcagatagcgact                                                >  1:3628143/1‑98 (MQ=11)
                                       cccagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGAtt                                                >  1:4364540/3‑100 (MQ=11)
                                       cacagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGACT‑CTGATGCGGACAGCGActcagacgcagatagcgactcagacgcagatagcgacg                                                >  1:6453800/2‑98 (MQ=11)
                                       aacggcgactcagacgcagaTAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGAtt                                                >  1:2301942/5‑100 (MQ=11)
                                        ctagcgactctgatgcggaCAGCGACTCAGA‑CGCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCAGACAGCGACTc                                               >  1:2219100/3‑100 (MQ=2)
                                        ccagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTc                                               >  1:1402954/2‑100 (MQ=11)
                                        acagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATacagacagatcgga                                               >  1:75982/1‑94 (MQ=2)
                                         cagcgattcagacgcagaTAGCGATTCAGA‑TTCAGATAACGACT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGAATca                                              >  1:3452218/1‑99 (MQ=2)
                                          agcgactcagacgcagaTAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGGTAGCGATTcag                                             >  1:5854696/1‑98 (MQ=11)
                                          agcgactcagacgcagaTAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTcag                                             >  1:4887167/1‑98 (MQ=11)
                                          agcgactcagacgcagaTAGCGACTCTGA‑TGCGGACAGCGACT‑CAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGATAGCGATTcag                                             >  1:1395631/1‑98 (MQ=11)
                                              attcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCGGAcgc                                         <  1:2704333/98‑3 (MQ=2)
                                               ggcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCGGAcaca                                        >  1:3761248/3‑100 (MQ=11)
                                                tcagacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTcaggcgcag                                       >  1:2569269/1‑92 (MQ=2)
                                                  tgacgcagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCGGACGCAGAt                                     >  1:71438/2‑100 (MQ=2)
                                                      ccagaCAGCGATTCAGA‑CGCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTTAGCTTCAGATAGCGATTCTGATGCGGATAGCg                                 >  1:2944327/2‑100 (MQ=2)
                                                       cagaTAGTGACTCAGA‑TTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGACTCAGATAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGa                                >  1:3742993/1‑100 (MQ=255)
                                                       cagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCTGACGCAGATAGCGa                                >  1:1108026/1‑100 (MQ=2)
                                                       cagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATGCGGACAGCGATTCGGACGCAGAtcggaa                                >  1:4824215/1‑95 (MQ=2)
                                                       cagaTAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGATTCGAACGCAGAtcggaa                                >  1:3747470/1‑95 (MQ=2)
                                                          ctAGCGATTCAGA‑TTCAGATAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGCTTCAGACAGCGACTCGGACGCAGATAGCGATtc                             >  1:1923951/2‑100 (MQ=11)
                                                              tGACTCAGA‑CGCAGATAGCGATT‑CAGATGCGGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATtcagat                         >  1:4609483/1‑100 (MQ=2)
                                                                     agacgcagatAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGAGAGCGATTCTGATGCAGATAGCGATtcagacgcagat                   >  1:6468681/6‑100 (MQ=2)
                                                                       acccagatAGCGATT‑CTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCCGACGCAGATAGCGTTTCAGACGCAGATAg                 >  1:2027464/4‑100 (MQ=2)
                                                                                        |                                                                                                  
TCAGATAGCGACTCAGATTCAGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGTGACTCAGA‑TTCAGATAGCGACT‑CAGATTCGGACAGCGATTCAGACTCAGATAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGATAGTGATTCAGACTCAGA  >  NC_007793/617954‑618138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑