Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 460,991 G→A 67.9% A702T (GCA→ACA)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913460,9910GA67.9% 109.9 / 85.1 100A702T (GCA→ACA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base A (43/25);  total (58/42)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 103 total aligned reads displayed.

CCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGC  >  NC_000913/460892‑461090
                                                                                                   |                                                                                                   
ccGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTgg                                                                                                     >  1:359631/1‑100 (MQ=255)
   gAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACaa                                                                                                  <  1:1886031/100‑1 (MQ=255)
     tcttGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATg                                                                                                <  2:1100103/97‑1 (MQ=255)
       gTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGAc                                                                                              >  1:648107/1‑100 (MQ=255)
        tGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGaga                                                                                             >  1:1100101/1‑98 (MQ=255)
        tGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGAcc                                                                                             >  1:2245015/1‑100 (MQ=255)
         gCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCg                                                                                            >  2:448202/1‑100 (MQ=255)
          cGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCgg                                                                                           >  1:2399475/1‑100 (MQ=255)
          cGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCg                                                                                            >  1:2269013/1‑99 (MQ=255)
            aCGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTg                                                                                         <  2:2863395/100‑1 (MQ=255)
            aCGCCGAAAGAGGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTg                                                                                         <  1:2638482/100‑1 (MQ=255)
             cGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTga                                                                                        <  2:2732299/100‑1 (MQ=255)
               ccGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTgaga                                                                                      >  2:1605914/1‑100 (MQ=255)
                cGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGAt                                                                                     <  2:2205829/100‑1 (MQ=255)
                cGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGAt                                                                                     >  2:2466570/1‑100 (MQ=255)
                cGAAAGATGGTCCCAGTGCCGGTATTGCTATGTGCACCTCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGAt                                                                                     >  1:1798133/1‑100 (MQ=255)
                 gAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATc                                                                                    <  1:1434738/100‑1 (MQ=255)
                     gATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCActc                                                                                >  2:879483/1‑100 (MQ=255)
                      aTGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCActct                                                                               >  2:619733/1‑100 (MQ=255)
                         gTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCg                                                                            <  2:1988820/100‑1 (MQ=255)
                          tCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGt                                                                           <  1:1731796/100‑1 (MQ=255)
                           ccGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTg                                                                          <  2:326939/100‑1 (MQ=255)
                           ccGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTg                                                                          <  1:360385/100‑1 (MQ=255)
                           ccGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTg                                                                          >  1:2347354/1‑100 (MQ=255)
                            cGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTgg                                                                         <  2:2435082/100‑1 (MQ=255)
                             gAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGt                                                                        <  2:418612/100‑1 (MQ=255)
                               gTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCa                                                                      >  2:197971/1‑100 (MQ=255)
                                tGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTgg                                                                         <  1:1210918/96‑1 (MQ=255)
                                  ccGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGt                                                                   >  1:1096742/1‑100 (MQ=255)
                                    ggTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTAc                                                                 <  2:692700/100‑1 (MQ=255)
                                       aTTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCATGTACTGc                                                              <  2:1527251/100‑1 (MQ=255)
                                       aTTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGc                                                              >  1:2403165/1‑100 (MQ=255)
                                        ttGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGcc                                                             <  2:2107127/100‑1 (MQ=255)
                                        ttGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGcc                                                             <  1:1866455/100‑1 (MQ=255)
                                         tGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCg                                                            >  1:1443629/1‑100 (MQ=255)
                                         tGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCg                                                            >  1:970126/1‑100 (MQ=255)
                                          gCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGccag                                                           >  2:2462629/1‑98 (MQ=255)
                                          gCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGCGATGACTCTGCGTGGTCAGGTACTGCCg                                                            >  1:2408970/1‑99 (MQ=255)
                                          gCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGa                                                           >  2:2714339/1‑100 (MQ=255)
                                              tgtgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCgg                                                       >  1:353445/1‑100 (MQ=255)
                                              agtgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATc                                                         >  1:552769/2‑98 (MQ=255)
                                               gtgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCggt                                                      >  1:2195837/1‑100 (MQ=255)
                                               gtgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCggt                                                      <  1:2791083/100‑1 (MQ=255)
                                                tgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCCCTCTGCGTGGTCAGGTACTGCCGATCggtg                                                     <  2:2765150/100‑1 (MQ=255)
                                                tgCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCggtg                                                     >  1:1985976/1‑100 (MQ=255)
                                                 gCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCggtgg                                                    >  1:2787715/1‑100 (MQ=255)
                                                   aCCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGtt                                                  <  2:773537/100‑1 (MQ=255)
                                                   aCCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGtt                                                  >  1:2192256/1‑100 (MQ=255)
                                                   aCCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTACGTGGTCAGGTACTGCCGATCGGTGGtt                                                  >  2:1627323/1‑100 (MQ=255)
                                                    ccGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGttt                                                 <  1:1951394/100‑1 (MQ=255)
                                                      gcgcTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTga                                               >  2:456603/1‑100 (MQ=255)
                                                        gcTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaa                                             <  2:341155/100‑1 (MQ=255)
                                                        gcTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaa                                             <  1:581452/100‑1 (MQ=255)
                                                        gcTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaa                                             <  1:448200/100‑1 (MQ=255)
                                                         cTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaag                                            <  1:2006860/100‑1 (MQ=255)
                                                         cTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTACGTGGTCAGGTACTGCCGATCGGTGGTTTgaaag                                            <  1:1627319/100‑1 (MQ=255)
                                                            gTTTCTTACCTAACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTgaaagaaa                                         >  1:1521017/1‑100 (MQ=255)
                                                             tttCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGaaaa                                        >  2:59075/1‑100 (MQ=255)
                                                                cTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACt                                     >  2:2665330/1‑100 (MQ=255)
                                                                 ttGCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTTAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTc                                    >  2:2625065/1‑100 (MQ=255)
                                                                   gCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCt                                  <  2:2095802/100‑1 (MQ=255)
                                                                   gCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCt                                  >  1:824991/1‑100 (MQ=255)
                                                                   gCCTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCt                                  >  2:696944/1‑100 (MQ=255)
                                                                    ccTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTg                                 <  1:1954234/100‑1 (MQ=255)
                                                                     cTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGTTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTg                                 >  1:2697255/1‑99 (MQ=255)
                                                                     cTGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTgg                                >  2:1620823/1‑100 (MQ=255)
                                                                      tGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAGCTCCTGGc                               <  2:1764307/100‑1 (MQ=255)
                                                                      tGACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGc                               >  2:1721700/1‑100 (MQ=255)
                                                                       gACCGGTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCa                              >  2:190084/1‑100 (MQ=255)
                                                                        aCCGGTAACCCGGTTCGTGCCGCTGTGACAATGACCGGAGACATCACTCTTCGTGGTCAGGAACTGCCGCTCCGTGGTTTGAAAGGAAAACTCCTGGCAg                             >  1:2764372/1‑100 (MQ=255)
                                                                         ccGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAgc                            >  1:2060646/1‑100 (MQ=255)
                                                                          cGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGATCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAgcg                           >  2:731952/1‑100 (MQ=255)
                                                                           ggTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAgcgc                          >  1:2083663/1‑100 (MQ=255)
                                                                            gTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCa                         >  2:1691530/1‑100 (MQ=255)
                                                                            gTAACCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGTAgcg                           >  1:1435862/1‑98 (MQ=255)
                                                                               aCCCGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcg                      >  1:261553/1‑100 (MQ=255)
                                                                                cccGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                     <  2:674366/100‑1 (MQ=255)
                                                                                cccGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                     >  2:2270049/1‑100 (MQ=255)
                                                                                cccGGTTCGTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATcgc                     >  2:1356130/1‑100 (MQ=255)
                                                                                      tCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCggg               <  1:10715/100‑1 (MQ=255)
                                                                                      tCGTGCCGATGTGGCAATGACCGGTGAGATCACACTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCggg               <  1:2751719/100‑1 (MQ=255)
                                                                                        gTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGAt             >  1:383769/1‑100 (MQ=255)
                                                                                        gTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGAt             >  2:1875307/1‑100 (MQ=255)
                                                                                        gTGCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGAt             >  1:204874/1‑100 (MQ=255)
                                                                                        gTGCCGAGGTGGCAATGACCGGTGAGATCAATCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGACAAACTCCTGGCAGCGCATCGCGGCGGGAt             >  1:1525636/1‑100 (MQ=255)
                                                                                          gCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTa           <  2:2195845/100‑1 (MQ=255)
                                                                                          gCCGATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTa           >  2:1783770/1‑100 (MQ=255)
                                                                                             gATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaaaa        <  2:175654/100‑1 (MQ=255)
                                                                                             gATGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaaaa        >  1:2312548/1‑100 (MQ=255)
                                                                                              aTGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTaaa         <  1:1766459/98‑1 (MQ=255)
                                                                                              aTGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAAc       <  2:2136462/100‑1 (MQ=255)
                                                                                              aTGTGACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAAc       <  2:2275280/100‑1 (MQ=255)
                                                                                              aTGTGACAATGACCGGTGAGATCACTCTGCGTCGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAAc       <  1:1135574/100‑1 (MQ=255)
                                                                                               tgtgACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACa      >  2:568556/1‑100 (MQ=255)
                                                                                                gtgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAg     <  2:1443632/100‑1 (MQ=255)
                                                                                                gtgACAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAg     <  1:164521/100‑1 (MQ=255)
                                                                                                 tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAg     >  1:1188894/1‑99 (MQ=255)
                                                                                                 tgGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGt    >  2:2834432/1‑100 (MQ=255)
                                                                                                   gCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGc  >  2:2151198/1‑100 (MQ=255)
                                                                                                   gCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGc  <  2:78919/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTGCCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGC  >  NC_000913/460892‑461090

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: