Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,657 Δ1 bp 6.6% coding (295/654 nt) rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,194,6570T. 6.6% 337.8 / 15.9 106coding (295/654 nt)rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base T (51/47);  new base . (2/5);  total (53/53)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 106 total aligned reads displayed.

CGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACA  >  NC_000913/1194510‑1194805
                                                                                                                                                    |                                                                                                                                                    
cgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTTTGTGGGAATCCCTTGCACCTGCACATTATAGATTTTACCGGTGCGTTT‑ACCCGGCTGGGTTa                                                                                                                                                    >  2:87373/1‑150 (MQ=255)
  ggTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTaaa                                                                                                                                                  >  1:748727/1‑150 (MQ=255)
  ggTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTaaa                                                                                                                                                  <  2:1554841/150‑1 (MQ=255)
    ttCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcg                                                                                                                                                <  2:582439/150‑1 (MQ=255)
     tCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgc                                                                                                                                               >  1:244704/1‑150 (MQ=255)
       gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACgcgc                                                                                                                                             >  2:1250854/1‑150 (MQ=255)
       gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACgcgc                                                                                                                                             >  1:216357/1‑150 (MQ=255)
        tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAACGCAcc                                                                                                                                            <  2:692161/150‑1 (MQ=255)
        tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                            >  2:1235947/1‑150 (MQ=255)
         cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTACTGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAACGCAccg                                                                                                                                           <  1:1286097/150‑2 (MQ=255)
         cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGG‑TAAACGCAccg                                                                                                                                           <  1:244318/150‑2 (MQ=255)
         cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGCCTGGGTTAAACGCGCCt                                                                                                                                           >  1:667474/1‑150 (MQ=255)
         cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAccg                                                                                                                                           <  2:1426579/150‑2 (MQ=255)
          aaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGGCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAccgg                                                                                                                                          <  1:1187161/150‑3 (MQ=255)
               gCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCCCCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcg                                                                                                                                     >  1:276345/1‑150 (MQ=255)
                cTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCCCCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgc                                                                                                                                    >  1:1489625/1‑150 (MQ=255)
                 tcagcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCGATTGCGCAAGGCTTCAAGTCCGTCTTGTGAGGGAATAACTTCCACCTGCACATAAAAGATTTTTCCGGTGCGTAT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                   >  1:359969/4‑150 (MQ=255)
                 tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                   <  1:805230/150‑1 (MQ=255)
                 tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                   <  1:1428191/150‑1 (MQ=255)
                 tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCAACTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgag                                                                                                                                   >  1:975379/1‑148 (MQ=255)
                 tCCGCGCCGGCGGGCAGGGTAGGGCCAACATTAAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATCCCTTCCACCTGCACAGAATAGATATTTCCGGTGCGATT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                   >  1:226152/1‑150 (MQ=255)
                    gcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCg                                                                                                                                <  2:1187566/150‑1 (MQ=255)
                      gcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGATTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGtt                                                                                                                              >  1:48524/1‑150 (MQ=255)
                       ccGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTa                                                                                                                             >  1:1647887/1‑150 (MQ=255)
                         ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                           >  2:635720/1‑150 (MQ=255)
                           cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                         <  2:574724/150‑1 (MQ=255)
                             ggCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                                       <  2:748727/150‑1 (MQ=255)
                              gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCTCAAGGCTTCATGTGCGTCTAGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                      <  1:96104/150‑1 (MQ=255)
                              gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                      <  2:1098572/150‑1 (MQ=255)
                               cAGGGTAGGGCCCTCATTTAAGGTGAGGCCATTGCGCAAGGCTTCAAGGGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAg                                                                                                                     <  2:266552/150‑1 (MQ=255)
                                aGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGc                                                                                                                    <  2:394882/150‑1 (MQ=255)
                                aGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGc                                                                                                                    <  2:1390501/150‑1 (MQ=255)
                                  ggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAc                                                                                                                  <  2:31631/150‑1 (MQ=255)
                                        gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                            >  2:673429/1‑150 (MQ=255)
                                           atcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAccc                                                                                                         <  1:221154/150‑1 (MQ=255)
                                               tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCCGCACCAGCAACCATTCg                                                                                                     >  2:1244693/1‑150 (MQ=255)
                                               tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                     <  2:1328432/150‑1 (MQ=255)
                                               tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCAACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGACGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                     >  1:65934/1‑150 (MQ=255)
                                               tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                     <  1:1721341/150‑1 (MQ=255)
                                                   aGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTAt                                                                                                 >  1:1620615/1‑150 (MQ=255)
                                                     gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                               >  1:1624676/1‑150 (MQ=255)
                                                     gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                               >  1:1449523/1‑150 (MQ=255)
                                                        aCGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCgg                                                                                            >  2:1667066/1‑150 (MQ=255)
                                                         cGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGt                                                                                           >  1:540633/1‑150 (MQ=255)
                                                         cGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGt                                                                                           >  2:1343686/1‑150 (MQ=255)
                                                          gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTc                                                                                          <  2:526614/150‑1 (MQ=255)
                                                           ccATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCTCTATCGCGGTCa                                                                                         >  2:1182241/1‑150 (MQ=255)
                                                           ccATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCa                                                                                         >  2:1771198/1‑150 (MQ=255)
                                                                gcgcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCg                                                                                    <  1:107272/150‑1 (MQ=255)
                                                                  gcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGAc                                                                                  <  2:1526253/150‑1 (MQ=255)
                                                                    aaGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACct                                                                                <  1:128906/150‑1 (MQ=255)
                                                                       gCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgct                                                                             <  2:882606/150‑1 (MQ=255)
                                                                        cTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgctg                                                                            >  1:1366160/1‑150 (MQ=255)
                                                                        cTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGAAACCCTTCGCTATCGCGGTCAAGGCGACctgctg                                                                            >  1:1549262/1‑150 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           >  2:1436445/1‑150 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           <  1:1171377/150‑1 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           <  2:58738/150‑1 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           <  1:1743181/150‑1 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           >  2:1595795/1‑150 (MQ=255)
                                                                         ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                           <  2:244704/150‑1 (MQ=255)
                                                                          tCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCa                                                                          <  2:1029099/150‑1 (MQ=255)
                                                                          tCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGGCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCa                                                                          <  2:1021908/150‑1 (MQ=255)
                                                                            aaGTGCGTCTTGTGTGGTAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATa                                                                        <  2:1386716/150‑1 (MQ=255)
                                                                                gCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCTCGGTCAAGGCGACCTGCTGCATAAacc                                                                    >  2:1660952/1‑149 (MQ=255)
                                                                                   tCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAccc                                                                 <  2:1461994/150‑1 (MQ=255)
                                                                                   tCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAccc                                                                 <  1:131185/150‑1 (MQ=255)
                                                                                       gtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGaa                                                             >  2:1248854/1‑150 (MQ=255)
                                                                                          tgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCg                                                          <  1:281089/150‑1 (MQ=255)
                                                                                          tgGGAATACCTTCCACCTGCACATAATAGAGTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGCCGACCTGCTGCATAAACACCCTGAACCg                                                          >  1:295324/1‑150 (MQ=255)
                                                                                            ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                        >  2:821859/1‑150 (MQ=255)
                                                                                             gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                                       <  1:723348/150‑1 (MQ=255)
                                                                                                    ttCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTc                                                >  1:452218/1‑150 (MQ=255)
                                                                                                         ccTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                           >  1:1067460/1‑150 (MQ=255)
                                                                                                         ccTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                           >  1:709677/1‑150 (MQ=255)
                                                                                                         ccTGCACATAATAGATTTTTCCGGTGCGGTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                           <  2:1099926/150‑1 (MQ=255)
                                                                                                           tGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCCCCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAATCACCCTGCACCGGGATGAATTCTTTTAAt                                         >  1:1253292/1‑150 (MQ=255)
                                                                                                            gCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATg                                        <  2:1509319/150‑1 (MQ=255)
                                                                                                                ataataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCAGTTATTGGTCAGCACCAGCAACCCTTCGCTATAGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                                    >  1:77164/1‑150 (MQ=255)
                                                                                                                 taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCCGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACACTGAACCGGGATGAATTCTTTTAATGTTTTg                                   <  1:1122550/150‑1 (MQ=255)
                                                                                                                    taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                                >  2:1133783/1‑150 (MQ=255)
                                                                                                                     aGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                               >  2:708528/1‑150 (MQ=255)
                                                                                                                     aGATTTTTCCGGAGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                               <  1:936839/150‑1 (MQ=255)
                                                                                                                      gATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTcc                              >  1:552245/1‑150 (MQ=255)
                                                                                                                        tttttCAGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                            >  1:789113/1‑150 (MQ=255)
                                                                                                                             ccGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGATTTCTTTTATTGTTTTGCGTCCGGCTTCa                       >  1:427995/1‑150 (MQ=255)
                                                                                                                             ccGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCa                       >  1:1612355/1‑150 (MQ=255)
                                                                                                                             caggcgCGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTTTTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAGTCTTTGAATGTTTTGCGTCCGGCTTCa                       >  1:1653449/6‑150 (MQ=255)
                                                                                                                              cGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGCACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCAt                      >  1:1107148/1‑150 (MQ=255)
                                                                                                                                gTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGACAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTCCATCg                    >  1:509049/1‑150 (MQ=255)
                                                                                                                                   cGTTTAACCCAGCCGGG‑TAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGCACCGGGATGAATTCTTTTTATGTTTTGCGTCCGGCTTCATCGGTg                 >  1:724781/1‑150 (MQ=255)
                                                                                                                                   cGTTT‑ACCCGGCTGGGTTTAACGCGCCTGCAACGCGCCGTTATGGGGCTGCACCAGCACCGCTGCGCTATCTCGGTCACGGCGACCTGCTTCATAAACTCCCTTAACCGTGATGAATTCTTTTAATGTTTTGCTTCCGGCTTCATCGGTg                 <  1:1114632/150‑1 (MQ=255)
                                                                                                                                   cGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAATCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg                 <  1:64146/150‑1 (MQ=255)
                                                                                                                                     ttt‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGACCCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGaa               >  1:926488/1‑150 (MQ=255)
                                                                                                                                       t‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt             <  1:1554326/150‑1 (MQ=255)
                                                                                                                                          cccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc           >  2:19290/1‑150 (MQ=255)
                                                                                                                                               cTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGgaag      <  1:1786819/150‑1 (MQ=255)
                                                                                                                                                tGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGTGgaaga     <  2:1449523/150‑1 (MQ=255)
                                                                                                                                                 gggTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGgaagaa    <  2:1366160/150‑1 (MQ=255)
                                                                                                                                                   gTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACa  <  1:706210/150‑1 (MQ=255)
                                                                                                                                                   gTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTAGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACa  >  1:863185/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACA  >  NC_000913/1194510‑1194805

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: ‑