Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 756,886 T→A 55.9% L327Q (CTG→CAG)  sdhA → succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913756,8860TA55.9% 69.5 / 216.2 163L327Q (CTG→CAG) sdhAsuccinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base A (43/48);  total (76/87)
Fisher's exact test for biased strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GATCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCGGT  >  NC_000913/756788‑756983
                                                                                                  |                                                                                                 
 aTCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGc                                                                                                 >  1:537231/1‑100 (MQ=255)
  tCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGcc                                                                                                >  2:2692970/1‑100 (MQ=255)
  tCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGcc                                                                                                >  2:4263220/1‑100 (MQ=255)
   cGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCg                                                                                               <  2:386004/100‑1 (MQ=255)
   cGAAATCCGTGAAGGCCGCGGCTGTGATGGTCCGTGGGGGCCACACCCGAAACTGAAACGCGATCACCTGGGTAAAGAAGTTCTCGACTCCCGTCAGc                                                                                                 <  1:2713234/98‑1 (MQ=255)
    gAAAGCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCgg                                                                                              <  1:1567742/100‑1 (MQ=255)
     aaaTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCggg                                                                                             <  2:1296074/100‑1 (MQ=255)
       aTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTa                                                                                           >  2:909559/1‑100 (MQ=255)
       aTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTa                                                                                           >  2:1187188/1‑100 (MQ=255)
         ccGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATc                                                                                         >  2:2374361/1‑100 (MQ=255)
         ccGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATc                                                                                         <  1:2553145/100‑1 (MQ=255)
          cGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATcc                                                                                        >  2:4107891/1‑100 (MQ=255)
           gTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCt                                                                                       >  2:3828526/1‑100 (MQ=255)
           gTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCt                                                                                       >  2:4130801/1‑100 (MQ=255)
           gTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCt                                                                                       <  2:653843/100‑1 (MQ=255)
           gTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAAATGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCt                                                                                       >  2:3761524/1‑100 (MQ=255)
            tGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTg                                                                                      >  2:4136688/1‑100 (MQ=255)
                   cgcgGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCttt                                                                               >  2:3932644/1‑100 (MQ=255)
                     cgGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTcc                                                                             <  2:3418576/100‑1 (MQ=255)
                     cgGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTcc                                                                             <  1:3670253/100‑1 (MQ=255)
                      ggCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGCTCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGCGCTTTccc                                                                            >  2:1216672/1‑100 (MQ=255)
                      ggCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTccc                                                                            <  2:3245120/100‑1 (MQ=255)
                      ggCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTccc                                                                            >  2:614121/1‑100 (MQ=255)
                       gCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCg                                                                           <  1:2060319/100‑1 (MQ=255)
                       gCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCg                                                                           >  1:646375/1‑100 (MQ=255)
                         tgtgATGGACCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTa                                                                         <  2:1104570/100‑1 (MQ=255)
                           tgATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTAcc                                                                       <  1:1045255/100‑1 (MQ=255)
                           tgATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTAcc                                                                       >  1:1230490/1‑100 (MQ=255)
                           tgATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTAcc                                                                       <  1:748336/100‑1 (MQ=255)
                             aTGGGCCGTGGGGGCTACACGCGAAACTGAAACTCGATCACGTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCtt                                                                     <  1:388996/100‑1 (MQ=255)
                              tGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTc                                                                    >  1:3950551/1‑100 (MQ=255)
                              tGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTc                                                                    >  2:1659730/1‑100 (MQ=255)
                              tGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTc                                                                    >  2:2811725/1‑100 (MQ=255)
                               ggTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGCGCTTTCCCGTACCTTCg                                                                   >  2:2062415/1‑100 (MQ=255)
                               ggTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCg                                                                   >  1:3460089/1‑100 (MQ=255)
                                 tCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCt                                                                 >  2:872596/1‑100 (MQ=255)
                                 tCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCt                                                                 <  1:1688856/100‑1 (MQ=255)
                                  ccGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTc                                                                >  1:2242183/1‑100 (MQ=255)
                                     tGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACg                                                             <  2:1917227/100‑1 (MQ=255)
                                       ggggCCACACGCGAAACTGAAACTCGCTCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTc                                                           >  2:2598086/1‑100 (MQ=255)
                                       ggggCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTc                                                           <  1:3128801/100‑1 (MQ=255)
                                           ccACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATc                                                       >  2:1373167/1‑100 (MQ=255)
                                           ccACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATc                                                       <  1:1834114/100‑1 (MQ=255)
                                           ccACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATc                                                       <  1:1341187/100‑1 (MQ=255)
                                           ccACACGCGAAACTGAAACTCGATCACCTGGGCAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCAGGAGCTTTCCCGTGCCTTCGCTCACGTCGATc                                                       >  1:1871469/1‑100 (MQ=255)
                                             acacGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCg                                                     <  2:1069646/100‑1 (MQ=255)
                                               acGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGt                                                   <  2:1505403/100‑1 (MQ=255)
                                                  cgAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGaa                                                <  2:3295209/100‑1 (MQ=255)
                                                   gaaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGaaa                                               >  2:1439518/1‑100 (MQ=255)
                                                    aaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAag                                              >  2:1453823/1‑100 (MQ=255)
                                                    aaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAag                                              <  1:3483450/100‑1 (MQ=255)
                                                     aactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAaga                                             >  2:1113274/1‑100 (MQ=255)
                                                      actgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAagag                                            >  1:1375124/1‑100 (MQ=255)
                                                      actgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAagag                                            <  1:690361/100‑1 (MQ=255)
                                                      actgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGCGAAagag                                            <  2:332715/100‑1 (MQ=255)
                                                         gtacctCGAGCAGCTGGGTATAGAAGTGCTCGAGTCCCGGCAGCCGGGTATCCTGGAGCTTGCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                         <  1:3828509/96‑1 (MQ=255)
                                                         gaaactCGATCAGCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                         >  2:2805096/1‑100 (MQ=255)
                                                         gaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                         >  2:3202865/1‑100 (MQ=255)
                                                           aactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGAt                                       <  1:992299/100‑1 (MQ=255)
                                                            actCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGAtt                                      >  2:3186702/1‑100 (MQ=255)
                                                             cgcGATCACCTGGATAATGAAGTTCCCGAATCCCGTCTGCCGAGTATCCTGGAGCTTTCCCGTAACTTCGCGCACGTCGATCTGGTGAAAGAGCCg                                         <  1:3693001/94‑1 (MQ=255)
                                                               cgagtacatGGGGAAAGAAGTGCTCGAACCCCGTCTGCCGGGTAAACTGGAGCGTCCCCATACCTTCGCTCACGGCGATCCGGTGAAagag                                            <  1:1406235/83‑1 (MQ=255)
                                                                gATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTAACTTCGCTCACATCGATCCGGTGAAAGAGCCGATTCCgg                                  >  1:2674860/1‑100 (MQ=255)
                                                                gATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCgg                                  >  1:1092338/1‑100 (MQ=255)
                                                                gATCACCTGCGCAAAAAAGTTATAGAATCCCGTCAGCCGGGTATCCTGGAGCTTTACCGTACCTTCACTCACGTCGATCCGGTGAAAGCGCCGATTCCg                                   >  1:860539/1‑99 (MQ=255)
                                                                 aTCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGt                                 <  2:2969345/100‑1 (MQ=255)
                                                                   cACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTa                               >  2:2441274/1‑100 (MQ=255)
                                                                    aCCTGGTGAAAGAAGTACTACAATCCCTACTGCCAGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGAGAAAGAGCCGATTCCgg                                  <  1:2797231/96‑1 (MQ=255)
                                                                    aCCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTAt                              >  2:2297305/1‑100 (MQ=255)
                                                                    aCCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTAt                              >  1:4229411/1‑100 (MQ=255)
                                                                     ccTGGGTAAATAAGTTCTCGAATCCCGGCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATc                             <  2:362767/100‑1 (MQ=255)
                                                                      cTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATcc                            <  2:3289112/100‑1 (MQ=255)
                                                                       tGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATccc                           >  2:2016828/1‑100 (MQ=255)
                                                                       tGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATccc                           >  2:1865691/1‑100 (MQ=255)
                                                                        gggTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCa                          >  1:3190584/1‑100 (MQ=255)
                                                                          gTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAc                        <  2:1073316/100‑1 (MQ=255)
                                                                           tAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                       >  2:445439/1‑100 (MQ=255)
                                                                           tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                       <  2:2243061/100‑1 (MQ=255)
                                                                           tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                       <  1:3226187/100‑1 (MQ=255)
                                                                               gaagTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCTTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTc                   <  1:2738610/100‑1 (MQ=255)
                                                                               gaagTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTc                   <  1:891991/100‑1 (MQ=255)
                                                                                 agTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCAc                 <  2:2704714/100‑1 (MQ=255)
                                                                                 agTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCAc                 <  2:1142655/100‑1 (MQ=255)
                                                                                   ttCTCGAATCCCGTGAGCCGGGCATCGTGGAGCTTTCCCGTACCTTCGCGCACGTCGACCCGGTGAACGAGTCGATTCCAGCTATCCCGACCTGCCa                  >  1:2690016/1‑97 (MQ=255)
                                                                                   ttCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACt                <  1:3045254/99‑1 (MQ=255)
                                                                                   ttCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCCCACCTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACt                <  2:891189/99‑1 (MQ=255)
                                                                                   ttCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTa               >  1:1064304/1‑100 (MQ=255)
                                                                                    tctcGAATCCCGTCTGCCGGGTATCCTGGGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACa             <  2:424972/100‑1 (MQ=255)
                                                                                     ctcGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACa             <  2:1643260/100‑1 (MQ=255)
                                                                                     ctcGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACa             >  2:732671/1‑100 (MQ=255)
                                                                                      tcGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTaaa             >  2:1848025/1‑97 (MQ=255)
                                                                                       cGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatg           >  1:2869229/1‑100 (MQ=255)
                                                                                        gAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatga          >  1:1142156/1‑100 (MQ=255)
                                                                                          aTCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatgatg        <  1:3186690/100‑1 (MQ=255)
                                                                                          aTCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatgatg        <  2:1768301/100‑1 (MQ=255)
                                                                                           tCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATgg       <  2:3540898/100‑1 (MQ=255)
                                                                                             ccGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGc     <  1:1551704/100‑1 (MQ=255)
                                                                                             ccGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGc     >  2:2461161/1‑100 (MQ=255)
                                                                                             ccGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGc     >  1:719107/1‑100 (MQ=255)
                                                                                                tCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCGGt  <  2:1046201/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
GATCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCGGT  >  NC_000913/756788‑756983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: