Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,231,009 G→C 6.5% G219A (GGG→GCG)  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,0090GC 6.5% 271.2 / 11.1 93G219A (GGG→GCG) rluF23S rRNA pseudouridine(2604) synthase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (2/4);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 102 total aligned reads displayed.

TAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGA  >  NC_000913/4230860‑4231156
                                                                                                                                                     |                                                                                                                                                   
tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCTTGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGATCGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                     >  1:652628/1‑150 (MQ=255)
tAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTggg                                                                                                                                                     >  1:482817/1‑150 (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:1334696/1‑150 (MQ=255)
                                                     gTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                >  1:1329487/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)
                                                            ccGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAtta                                                                                         >  2:394312/1‑149 (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)
                                                              ggCGCATGTGCGAGCATTTCGACTATGAAGTGAAAAAGCTGGAACGCCCGCGCATCATGAACGTTAGCTTAAGCTTCATTCAGCTGGGGGAATTGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCct                                                                                       >  1:813808/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)
                                                                      tgCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAg                                                                               >  2:1248229/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                                                                             <  1:94533/150‑1 (MQ=255)
                                                                        cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             >  2:801962/1‑150 (MQ=255)
                                                                        cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                             >  2:162819/1‑150 (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)
                                                                                                                        gAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTcc                             <  1:23534/150‑1 (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  >  2:3695/1‑150 (MQ=255)
                                                                                                                                                   gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa  <  1:1067969/150‑1 (MQ=255)
                                                                                                                                                   gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCGTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGAATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa  >  2:1087089/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGA  >  NC_000913/4230860‑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: