Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,231,004 G→C 8.1% P217P (CCG→CCC rluF → 23S rRNA pseudouridine(2604) and tRNA(Tyr) pseudouridine(35) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,231,0040GC 8.1% 478.2 / 25.8 149P217P (CCG→CCCrluF23S rRNA pseudouridine(2604) and tRNA(Tyr) pseudouridine(35) synthase
Reads supporting (aligned to +/- strand):  ref base G (63/73);  new base C (5/7);  total (68/81)
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.84e-01

Only 100 of 151 total aligned reads displayed.

AAAGTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGC  >  NC_000913/4230855‑4231153
                                                                                                                                                     |                                                                                                                                                     
  aGTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCt                                                                                                                                                     <  1:1068929/150‑1 (MQ=255)
   gTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTg                                                                                                                                                    >  1:1290331/1‑150 (MQ=255)
     tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                  >  1:2203053/1‑150 (MQ=255)
     tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                  <  2:1907599/150‑1 (MQ=255)
      aaaaaaGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTgggg                                                                                                                                                 >  1:582261/1‑150 (MQ=255)
       aaaaaGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggggg                                                                                                                                                >  1:1412125/1‑150 (MQ=255)
              agCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccg                                                                                                                                         <  2:1967259/150‑1 (MQ=255)
              agCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCAGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAATCTCGAACGCACGCGAATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccg                                                                                                                                         <  2:2144538/150‑1 (MQ=255)
               gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGgcgc                                                                                                                                        <  1:1775673/150‑1 (MQ=255)
               gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGgcgc                                                                                                                                        <  2:733552/150‑1 (MQ=255)
               gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgc                                                                                                                                        <  1:321743/150‑1 (MQ=255)
               gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgc                                                                                                                                        <  1:1402204/150‑1 (MQ=255)
               gcgcCGTTTGGCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGgcgc                                                                                                                                        <  2:2485656/150‑1 (MQ=255)
                cgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                       <  2:797683/150‑2 (MQ=255)
                cgcCGGTTGTCTTCCGCATTACCCTGGTGCACGGGCTGAACCGTCAGATCCGGCGCATGGGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                       <  2:2013602/150‑2 (MQ=255)
                cgcCATTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGcgcg                                                                                                                                       <  2:146175/150‑1 (MQ=255)
                  ccGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                     >  1:1526023/1‑150 (MQ=255)
                  ccGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                     <  1:240835/150‑1 (MQ=255)
                   cGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAtt                                                                                                                                    <  2:1290331/150‑1 (MQ=255)
                       tGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAc                                                                                                                                >  2:663814/1‑150 (MQ=255)
                        gTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAcc                                                                                                                               <  2:1755664/150‑1 (MQ=255)
                          cTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcga                                                                                                                             >  2:2411404/1‑150 (MQ=255)
                           ttCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgac                                                                                                                            <  1:96635/150‑1 (MQ=255)
                           ttCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgac                                                                                                                            <  2:2328247/150‑1 (MQ=255)
                            tCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgacg                                                                                                                           <  1:973928/150‑1 (MQ=255)
                              cGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGAt                                                                                                                         <  2:2375060/150‑1 (MQ=255)
                               gCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATg                                                                                                                        >  2:644396/1‑150 (MQ=255)
                                 aTTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAg                                                                                                                      <  1:379064/150‑1 (MQ=255)
                                 aTTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAg                                                                                                                      >  1:1005760/1‑150 (MQ=255)
                                  ttACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGt                                                                                                                     <  2:496263/150‑1 (MQ=255)
                                   tACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAAGGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGtt                                                                                                                    <  2:561685/150‑1 (MQ=255)
                                     cccTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTaa                                                                                                                  >  1:1750244/1‑150 (MQ=255)
                                       cTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATc                                                                                                                <  2:2451681/150‑1 (MQ=255)
                                        tGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCg                                                                                                               <  1:1981878/150‑1 (MQ=255)
                                           tGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                                                                                            <  2:1878927/150‑1 (MQ=255)
                                               ggggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCtt                                                                                                        >  1:637718/1‑150 (MQ=255)
                                                   cTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAg                                                                                                    >  1:277232/1‑150 (MQ=255)
                                                    tGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGc                                                                                                   >  2:2267837/1‑150 (MQ=255)
                                                       aCCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCa                                                                                                <  1:1522709/150‑1 (MQ=255)
                                                         cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                              <  2:1707659/150‑1 (MQ=255)
                                                         cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                              >  1:439873/1‑150 (MQ=255)
                                                          gTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                             <  1:1204757/150‑1 (MQ=255)
                                                            cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaa                                                                                           >  1:523218/1‑150 (MQ=255)
                                                             aGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaa                                                                                          <  1:1144615/150‑1 (MQ=255)
                                                             aGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaa                                                                                          >  1:585436/1‑150 (MQ=255)
                                                              gATCCGGCGCATGTGCGAGCATTTCCGCTATGAAGTCAAAAAGCTGCAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCGGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaaa                                                                                         <  2:1670446/150‑1 (MQ=255)
                                                               aTCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAt                                                                                        <  1:471079/150‑1 (MQ=255)
                                                               aTCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAt                                                                                        >  1:1314866/1‑150 (MQ=255)
                                                                  cGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTcc                                                                                     <  2:1641075/150‑1 (MQ=255)
                                                                         tgtgCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCg                                                                              <  1:1275739/150‑1 (MQ=255)
                                                                            gCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAgg                                                                           <  1:246556/150‑1 (MQ=255)
                                                                            gCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAgg                                                                           <  2:589278/150‑1 (MQ=255)
                                                                             cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCTGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                              <  2:1154894/130‑1 (MQ=255)
                                                                             cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                          >  1:858668/1‑150 (MQ=255)
                                                                             cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCAGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                              >  1:1154894/1‑130 (MQ=255)
                                                                               aGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaa                                                                        >  2:543699/1‑150 (MQ=255)
                                                                                 cATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaaaa                                                                      >  1:225283/1‑150 (MQ=255)
                                                                                     tCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAGTTAACCGACGATGAGTTAATCGCCCTCGTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTa                                                                  >  2:2071189/1‑150 (MQ=255)
                                                                                      cGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa                                                                 <  1:1620665/150‑1 (MQ=255)
                                                                                       ggCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaaa                                                                <  1:2416920/150‑1 (MQ=255)
                                                                                         cTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGc                                                              <  1:858811/150‑1 (MQ=255)
                                                                                            tGAAGTGAAAAAGCTGGAACGCACTAGCATCATGAACGTTAGCTTAAGCGGCTTTCCTGTGGGGGAATGGCGCGATTTAACCGCCGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGATGTAAAACCTAAAGCGaa                                                           <  2:1015252/150‑1 (MQ=255)
                                                                                                 tGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCa                                                      >  1:1038622/1‑150 (MQ=255)
                                                                                                 tGAAAAAGCTGGAACGCACGAGCATCATGAACGTTAGCTTAAGCGGCATTCCGGTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACATAAAGCGAAGGCCa                                                      >  2:821901/1‑150 (MQ=255)
                                                                                                  gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaa                                                     >  1:1017184/1‑150 (MQ=255)
                                                                                                    aaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAc                                                   >  1:515361/1‑150 (MQ=255)
                                                                                                    aaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCCGCATTCCGCTGGGGGAATCGCGCGATTTAACCGAAGATGAGTTAATCGACATCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAc                                                   >  1:731538/1‑150 (MQ=255)
                                                                                                     aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAcc                                                  <  1:1775274/150‑1 (MQ=255)
                                                                                                           ggAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAAc                                            >  2:951855/1‑150 (MQ=255)
                                                                                                           ggAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAAc                                            >  1:1592595/1‑150 (MQ=255)
                                                                                                                gcacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCgg                                       <  2:1331240/150‑1 (MQ=255)
                                                                                                                 cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                      >  2:1476287/1‑150 (MQ=255)
                                                                                                                  acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                     <  1:2411404/150‑1 (MQ=255)
                                                                                                                  acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                     <  2:1801998/150‑1 (MQ=255)
                                                                                                                     cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATc                                  <  1:2290537/150‑1 (MQ=255)
                                                                                                                     cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATc                                  <  2:585436/150‑1 (MQ=255)
                                                                                                                         tcatGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAAc                              >  1:1323564/1‑150 (MQ=255)
                                                                                                                         tcatGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAAc                              >  1:2198896/1‑150 (MQ=255)
                                                                                                                          catGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACg                             <  2:1012361/150‑1 (MQ=255)
                                                                                                                                 gTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTc                      <  2:1059897/150‑1 (MQ=255)
                                                                                                                                     gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTACCCGACGATGAGTTAATCGCCCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACCGCGGGCATCAAACGTCCAGTCGTTa                  >  2:1371381/1‑150 (MQ=255)
                                                                                                                                     gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                  >  1:1681935/1‑150 (MQ=255)
                                                                                                                                     gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                  >  2:762328/1‑150 (MQ=255)
                                                                                                                                     gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                  >  2:1419234/1‑150 (MQ=255)
                                                                                                                                     gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                  >  1:698348/1‑150 (MQ=255)
                                                                                                                                      cTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTaa                 >  2:1097641/1‑150 (MQ=255)
                                                                                                                                      aaatcgcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCGAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTaa                 >  1:1599511/6‑150 (MQ=255)
                                                                                                                                       ttAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                <  2:65024/150‑1 (MQ=255)
                                                                                                                                       ttAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                >  2:1707293/1‑150 (MQ=255)
                                                                                                                                         tcgcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAt              >  2:1095690/3‑150 (MQ=255)
                                                                                                                                           gcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATgg            >  2:2381331/1‑150 (MQ=255)
                                                                                                                                           gcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATgg            >  1:2484779/1‑150 (MQ=255)
                                                                                                                                             ggCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaa          <  1:352714/150‑1 (MQ=255)
                                                                                                                                              gCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaa         >  1:2363308/1‑150 (MQ=255)
                                                                                                                                              gCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaa         >  1:1754819/1‑150 (MQ=255)
                                                                                                                                                    cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAATCAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg   <  1:1001702/150‑1 (MQ=255)
                                                                                                                                                    cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg   <  1:2016762/150‑1 (MQ=255)
                                                                                                                                                    cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg   >  1:1069974/1‑150 (MQ=255)
                                                                                                                                                    cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAAGGAAAAAAcgg   >  2:923113/1‑150 (MQ=255)
                                                                                                                                                     gCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggc  >  1:1251623/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGC  >  NC_000913/4230855‑4231153

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