Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,231,007 G→A 6.6% L218L (CTG→CTA) ‡ rluF → 23S rRNA pseudouridine(2604) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,231,0070GA 6.6% 266.4 / 6.9 93L218L (CTG→CTA) ‡rluF23S rRNA pseudouridine(2604) synthase
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base A (2/4);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 105 total aligned reads displayed.

GTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGA  >  NC_000913/4230858‑4231156
                                                                                                                                                     |                                                                                                                                                     
gTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTg                                                                                                                                                       <  2:1268972/150‑1 (MQ=255)
gTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTg                                                                                                                                                       <  2:1341347/150‑1 (MQ=255)
  tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCTTGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGATCGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                     >  1:652628/1‑150 (MQ=255)
  tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                     >  1:482817/1‑150 (MQ=255)
  tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGAGCAGGGGCTGAACCGTCTGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggt                                                                                                                                                     <  1:993787/150‑2 (MQ=255)
    aaaaaGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggggg                                                                                                                                                   <  2:533903/150‑1 (MQ=255)
    aaaaaGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCAGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggggg                                                                                                                                                   >  1:462706/1‑150 (MQ=255)
    aaaaaGAAGCGCCGTTTCGCGTCCGCATTACCCTGGGGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTAAGCTTAAGCGGCATTCCGCTggggg                                                                                                                                                   <  2:199242/150‑1 (MQ=255)
    aaaaaGAAGCGCCGTCTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggggg                                                                                                                                                   <  1:819206/150‑1 (MQ=255)
        agaagCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATg                                                                                                                                               <  1:364458/150‑1 (MQ=255)
           agCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCCCCAGCGGAATGccg                                                                                                                                            <  1:1060268/150‑1 (MQ=255)
             cgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                          <  1:928043/150‑2 (MQ=255)
             cgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                          <  1:390490/150‑2 (MQ=255)
              gcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGcgcgt                                                                                                                                         >  2:573065/1‑149 (MQ=255)
              gcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGa                                                                                                                                         >  2:1312188/1‑150 (MQ=255)
              gcCGGTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgctt                                                                                                                                         <  1:498716/150‑3 (MQ=255)
                 gTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAttt                                                                                                                                      >  1:1100382/1‑150 (MQ=255)
                    tGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGATGGGGGAATGGCGCGATTTAAc                                                                                                                                   >  1:1188315/1‑150 (MQ=255)
                     gacTTCCGCATTACCCTGGGGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGAGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCTTGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAcc                                                                                                                                  <  1:890915/148‑1 (MQ=255)
                      gctaccGCATTACCCTGGGGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcg                                                                                                                                 <  2:217838/146‑1 (MQ=255)
                            gCATTACCCTGGTGCAGGGGTTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATg                                                                                                                           >  2:1067969/1‑150 (MQ=255)
                            gCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATg                                                                                                                           <  2:7459/150‑1 (MQ=255)
                                tACCCTGGTGCTGGGGCTGAACATTCTGAACCGGAGCATGCGCGAGCATTGCGGCAATGAAGAGAAAACGCTGGAACGCATGGGCATCATGAACGTTAGCTAAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGtt                                                                                                                       <  2:877682/150‑1 (MQ=255)
                                tACCCTGGGGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTAAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGtt                                                                                                                       <  2:745258/150‑1 (MQ=255)
                                 aCCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTa                                                                                                                      >  2:69556/1‑150 (MQ=255)
                                  cccTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTGGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTACCCGACGATGAGTTaa                                                                                                                     >  2:797026/1‑150 (MQ=255)
                                    cTGGTGCAGGGGCTGAACCGTCCGCTCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAACCTGGAACGCACGCGCATCATGAACGTAAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATc                                                                                                                   <  2:150567/150‑1 (MQ=255)
                                    cTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCCTCATGAACGTTAGCTTACGCGGCATTCCGCTGGGGGAAGGGCGCGATTTAACCGACGATGAGTTAATc                                                                                                                   >  1:196055/1‑150 (MQ=255)
                                            ggggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCtt                                                                                                           >  1:76109/1‑150 (MQ=255)
                                             gggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCttt                                                                                                          <  1:732853/150‑1 (MQ=255)
                                              ggCTGAACCGTCAGATCCGGCGCATGTGCGATCATTTCGGCTATGAAGTGAAAAAGCTGGAACGAACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTa                                                                                                         <  1:809006/150‑1 (MQ=255)
                                                  gAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCt                                                                                                     >  1:311463/1‑150 (MQ=255)
                                                     ccGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAt                                                                                                  >  2:133286/1‑150 (MQ=255)
                                                      cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                                 <  1:554333/150‑1 (MQ=255)
                                                       gTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                >  1:1329487/1‑150 (MQ=255)
                                                       gTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                >  1:1334696/1‑150 (MQ=255)
                                                        tCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGa                                                                                               <  2:1119540/150‑1 (MQ=255)
                                                          aGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaa                                                                                             >  2:303982/1‑150 (MQ=255)
                                                             tCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAtt                                                                                          <  2:1301030/150‑1 (MQ=255)
                                                             tCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACAGACGATGAGTTAATCGACCTCTTTAAGCTCATTGATAAtt                                                                                          <  2:860480/150‑1 (MQ=255)
                                                              ccGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTc                                                                                         <  2:1340485/150‑1 (MQ=255)
                                                              ccGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTc                                                                                         >  2:1106309/1‑150 (MQ=255)
                                                               cGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTcc                                                                                        <  1:307050/150‑1 (MQ=255)
                                                                ggCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCct                                                                                       >  1:968378/1‑150 (MQ=255)
                                                                  cgcATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCACTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCctct                                                                                     >  2:1140743/1‑150 (MQ=255)
                                                                    cATGTGCGAGCATTTCGGCTATGACTTGAAAAAGCGGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTc                                                                                   <  2:222401/150‑1 (MQ=255)
                                                                    aaTGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTc                                                                                   <  2:654488/149‑1 (MQ=255)
                                                                       gtgCGAGCATTTCTGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCTACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGa                                                                                >  1:1139715/1‑150 (MQ=255)
                                                                        tgCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAACGACGATGAGTTAATCGCCCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAg                                                                               >  1:159944/1‑150 (MQ=255)
                                                                         gCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAgg                                                                              <  1:477627/150‑1 (MQ=255)
                                                                         gCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACAGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAgg                                                                              <  1:1005828/150‑1 (MQ=255)
                                                                          cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             >  2:1003296/1‑150 (MQ=255)
                                                                          cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             >  2:801962/1‑150 (MQ=255)
                                                                          cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             >  2:162819/1‑150 (MQ=255)
                                                                          cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             <  1:94533/150‑1 (MQ=255)
                                                                           taGCATTTCGGCTATGAATTGAATAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCCGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTa                                                                            <  1:134105/149‑1 (MQ=255)
                                                                             gCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaaa                                                                          <  1:1285220/150‑1 (MQ=255)
                                                                               aTTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAAc                                                                        <  2:404671/150‑1 (MQ=255)
                                                                                tttCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAAcc                                                                       <  1:69556/150‑1 (MQ=255)
                                                                                 ttCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCt                                                                      <  2:1163230/150‑1 (MQ=255)
                                                                                     gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTATCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAg                                                                  <  1:1176547/150‑1 (MQ=255)
                                                                                     gCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAg                                                                  >  2:1134518/1‑150 (MQ=255)
                                                                                       tATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCg                                                                >  2:400054/1‑150 (MQ=255)
                                                                                        aTGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGa                                                               >  2:391823/1‑150 (MQ=255)
                                                                                         cgAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGaa                                                              >  2:49252/2‑150 (MQ=255)
                                                                                          gAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTACTCGCCCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAg                                                             >  2:52537/1‑150 (MQ=255)
                                                                                          gAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAg                                                             <  2:765029/150‑1 (MQ=255)
                                                                                            aGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGc                                                           >  1:969367/1‑150 (MQ=255)
                                                                                               gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaa                                                        >  1:152213/1‑150 (MQ=255)
                                                                                                  aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAcc                                                     <  1:820785/150‑1 (MQ=255)
                                                                                                   aaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCg                                                    >  2:1116921/1‑150 (MQ=255)
                                                                                                      cTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGaaa                                                 >  1:646715/1‑150 (MQ=255)
                                                                                                       tGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGaaaa                                                >  2:303148/1‑150 (MQ=255)
                                                                                                        ggAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTTATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAAc                                               <  2:1327808/150‑1 (MQ=255)
                                                                                                          aaCGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAACACCTAAAGCGAAGGCCAAACCGAAAACAg                                             >  2:1021031/1‑150 (MQ=255)
                                                                                                            cgcacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCg                                           <  1:1246412/150‑1 (MQ=255)
                                                                                                              cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                         <  2:841753/150‑1 (MQ=255)
                                                                                                               acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                        <  1:695012/150‑1 (MQ=255)
                                                                                                               acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                        <  1:1312188/150‑1 (MQ=255)
                                                                                                               acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCATTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAACCCGAAAACAGCGGGc                                        <  1:198715/150‑1 (MQ=255)
                                                                                                                   gcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCa                                    <  1:1262634/150‑1 (MQ=255)
                                                                                                                   gcATCATGAACGTTAGCTTAAGCGGCATTACGCTGGGGGAATGGCGCGATTTACCCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCa                                    >  1:635225/1‑150 (MQ=255)
                                                                                                                     ttcatGAACGTTAGCTTCAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCaac                                  <  1:854589/149‑2 (MQ=255)
                                                                                                                      tcatGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAGTGAAAATTCCTCTTCCGAGGTAAAACCTAAATCGAAGGCACAACCGAAAACAGCGGGCATccacc                                 >  2:1353511/1‑146 (MQ=255)
                                                                                                                          gAACGTTTGGTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCGGATTGAAAATTCCTCTTCATAGGTAATACCTAAAGCGAAGTCCAAACCGAAAACAGCGGGCATCAAACGTcc                             <  2:534165/150‑1 (MQ=255)
                                                                                                                          gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTcc                             >  1:633771/1‑150 (MQ=255)
                                                                                                                          gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTcc                             >  1:1093646/1‑150 (MQ=255)
                                                                                                                             cGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGt                          >  1:708404/1‑150 (MQ=255)
                                                                                                                              gTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTc                         <  2:1294909/150‑1 (MQ=255)
                                                                                                                               ttAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCg                        <  2:76109/150‑1 (MQ=255)
                                                                                                                                    aatcgcgcCATTCCCCCAGCGGAATGCCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                   >  2:1073096/5‑150 (MQ=255)
                                                                                                                                     tAAGCGGCATTCCGCTGGGGAATGGCGCGATTTAACCGACGATGAGCTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAACCCGAAAACAGCGGGCCTCAACCGTCCAGGCGTTAAGAt                 >  1:1331977/1‑150 (MQ=255)
                                                                                                                                       cgcgcCATTCCCCCAGCGGAATGCCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATg                >  2:701603/2‑150 (MQ=255)
                                                                                                                                         cGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGa              <  1:533485/150‑1 (MQ=255)
                                                                                                                                            cATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaaa           >  2:404041/1‑150 (MQ=255)
                                                                                                                                               tCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAc        >  2:387022/1‑150 (MQ=255)
                                                                                                                                                  gCTGGGGGAATGTCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggc     <  2:240486/150‑1 (MQ=255)
                                                                                                                                                     gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa  <  1:1067969/150‑1 (MQ=255)
                                                                                                                                                     gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa  >  2:3695/1‑150 (MQ=255)
                                                                                                                                                     gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCGTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGAATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa  >  2:1087089/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGA  >  NC_000913/4230858‑4231156

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