Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,645:1 +A 5.8% coding (307/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,6451.A 5.8% 746.0 / 39.9 171coding (307/654 nt)rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base . (79/82);  new base A (4/6);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 177 total aligned reads displayed.

GGCCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  NC_000913/1194497‑1194794
                                                                                                                                                     |                                                                                                                                                     
    aTAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCg                                                                                                                                                  >  1:1926338/1‑150 (MQ=255)
       aCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCt                                                                                                                                               <  2:1857529/150‑1 (MQ=255)
        ccACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTg                                                                                                                                              <  2:1154673/150‑1 (MQ=255)
            gcgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGtt                                                                                                                                          >  2:1842278/1‑150 (MQ=255)
                   cGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAAcgcg                                                                                                                                   <  2:1114268/150‑1 (MQ=255)
                    gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACgcgc                                                                                                                                  >  1:2748399/1‑150 (MQ=255)
                    gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAAcgcac                                                                                                                                  <  1:1121392/150‑3 (MQ=255)
                    gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAAcgcac                                                                                                                                  <  1:1121584/150‑3 (MQ=255)
                     tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                 >  2:2622692/1‑150 (MQ=255)
                     tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                 >  2:428029/1‑150 (MQ=255)
                     tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                 >  2:1111442/1‑150 (MQ=255)
                     tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGcc                                                                                                                                 >  2:1394887/1‑150 (MQ=255)
                     tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCAGTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGG‑TAAACGCAcc                                                                                                                                 <  2:8730/150‑1 (MQ=255)
                        aCCAGCTCCGCGCCGGCGGGCAGGGTCGGGCCATCATGTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGGGTGGGACTACCTTCCACCTGCACCTCATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                              >  2:304215/1‑150 (MQ=255)
                        aCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGc                                                                                                                              >  2:4895/1‑150 (MQ=255)
                         ccAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCa                                                                                                                             >  2:210514/1‑150 (MQ=255)
                          cAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCaa                                                                                                                            <  2:687586/150‑1 (MQ=255)
                           aGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAc                                                                                                                           <  2:2484738/150‑1 (MQ=255)
                           aGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAc                                                                                                                           >  2:2499773/1‑150 (MQ=255)
                             cTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgc                                                                                                                         >  1:1902331/1‑150 (MQ=255)
                             cTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgc                                                                                                                         <  1:2152524/150‑1 (MQ=255)
                              tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                        >  1:697998/1‑150 (MQ=255)
                              tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGATACGCCATTGCGCAAGGCTTCAAGTGCGTCTAGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                        >  1:2683913/1‑150 (MQ=255)
                               ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACgcgc                                                                                                                       <  2:591413/150‑1 (MQ=255)
                               ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACgcgc                                                                                                                       <  2:2136750/150‑1 (MQ=255)
                               ccGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACgcgc                                                                                                                       <  2:1241260/150‑1 (MQ=255)
                                cgcgCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      <  2:1243793/150‑1 (MQ=255)
                                cgcgCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      >  1:2294602/1‑150 (MQ=255)
                                cgcgCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAAAAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGcc                                                                                                                      >  2:2353318/1‑150 (MQ=255)
                                      ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                <  1:2549105/150‑1 (MQ=255)
                                           gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                           <  1:1164804/150‑1 (MQ=255)
                                           gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                           <  1:1164657/150‑1 (MQ=255)
                                           gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTGTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                           <  2:1328712/150‑1 (MQ=255)
                                             aGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGGTATTGGTCAGc                                                                                                         >  2:1556373/1‑150 (MQ=255)
                                              gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCTAGTGAGTCTTGTGTGGGAATACCTTACACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                        <  2:1952062/150‑1 (MQ=255)
                                              gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                        <  1:2638528/150‑1 (MQ=255)
                                              gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                        <  2:465391/150‑1 (MQ=255)
                                              gggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCa                                                                                                        <  2:2855274/150‑1 (MQ=255)
                                               ggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAc                                                                                                       >  2:1372226/1‑150 (MQ=255)
                                                  aGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAg                                                                                                    <  2:2832117/150‑1 (MQ=255)
                                                     gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                 >  1:1455186/1‑150 (MQ=255)
                                                      ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAc                                                                                                <  2:2371853/150‑1 (MQ=255)
                                                      ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAc                                                                                                <  1:1277514/150‑1 (MQ=255)
                                                      ccATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAc                                                                                                <  1:1058414/150‑1 (MQ=255)
                                                           atTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAGCAGCAACCCTTc                                                                                           <  2:2391746/150‑1 (MQ=255)
                                                            tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                          <  1:2281007/150‑1 (MQ=255)
                                                            tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                          <  2:147133/150‑1 (MQ=255)
                                                            tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                          <  1:921457/150‑1 (MQ=255)
                                                                   ttACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTATTCGCGCCTGCAACGCGCGGTTATTGGTCAGCACCAGCAACCCTTCGCTATcgc                                                                                   >  2:53748/1‑150 (MQ=255)
                                                                    tACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcgcg                                                                                  <  2:24929/150‑1 (MQ=255)
                                                                       gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTc                                                                               >  1:2632639/1‑150 (MQ=255)
                                                                         cATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCaa                                                                             <  2:566370/150‑1 (MQ=255)
                                                                         cATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAAGACCTTCCACCTGCACATAATAGATTGTGCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCaa                                                                             <  2:1372757/150‑1 (MQ=255)
                                                                            tGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGc                                                                          <  1:1419615/150‑1 (MQ=255)
                                                                            tGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGc                                                                          <  1:2471904/150‑1 (MQ=255)
                                                                               gcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGAc                                                                       >  2:1206630/1‑150 (MQ=255)
                                                                                cAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGAcc                                                                      <  2:1280801/150‑1 (MQ=255)
                                                                                 aaGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACct                                                                     >  1:1961799/1‑150 (MQ=255)
                                                                                    gCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgct                                                                  <  2:215637/150‑1 (MQ=255)
                                                                                           gTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATaaa                                                           >  1:1611814/1‑150 (MQ=255)
                                                                                             gCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAaca                                                         >  1:1317971/1‑150 (MQ=255)
                                                                                             gCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAaca                                                         >  1:2901149/1‑150 (MQ=255)
                                                                                                   tgtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGa                                                   <  1:661304/150‑1 (MQ=255)
                                                                                                    gtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGaa                                                  <  1:857859/150‑1 (MQ=255)
                                                                                                       tgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCg                                               <  1:2452944/150‑1 (MQ=255)
                                                                                                        gggAATACCTTCAACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTACGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCgg                                              >  1:1477444/1‑150 (MQ=255)
                                                                                                         ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                             >  1:503496/1‑150 (MQ=255)
                                                                                                          gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGACCCGGGa                                            >  2:2826012/1‑150 (MQ=255)
                                                                                                          gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                            <  1:1897712/150‑1 (MQ=255)
                                                                                                          gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                            >  2:1757896/1‑150 (MQ=255)
                                                                                                          gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                            <  1:2717847/150‑1 (MQ=255)
                                                                                                          gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                            >  1:1398111/1‑150 (MQ=255)
                                                                                                           aaTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGAt                                           <  1:1372226/150‑1 (MQ=255)
                                                                                                             tACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGa                                         <  1:243731/150‑1 (MQ=255)
                                                                                                                    cACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCttt                                  <  1:1080856/150‑1 (MQ=255)
                                                                                                                      ccTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                >  1:2369121/1‑150 (MQ=255)
                                                                                                                       cTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                               <  2:131169/150‑1 (MQ=255)
                                                                                                                       cTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTaa                               <  2:1455186/150‑1 (MQ=255)
                                                                                                                        tGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAAt                              <  1:1733791/150‑1 (MQ=255)
                                                                                                                        cgCACACAATCGATTTTTCCGGCGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAAt                              <  1:1733838/149‑1 (MQ=255)
                                                                                                                         gCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCGGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGAGGAATTCTTTTAATg                             >  2:2428637/1‑150 (MQ=255)
                                                                                                                           acaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtt                           <  2:2894180/150‑1 (MQ=255)
                                                                                                                           acaTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtt                           <  2:1899380/150‑1 (MQ=255)
                                                                                                                            caTAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGttt                          <  2:2581544/150‑1 (MQ=255)
                                                                                                                             ataataGATTTGTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCGGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                         <  1:2082091/150‑1 (MQ=255)
                                                                                                                              taataGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTg                        <  2:544868/150‑1 (MQ=255)
                                                                                                                                 taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                     >  1:2360940/1‑150 (MQ=255)
                                                                                                                                 taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                     >  2:1975149/1‑150 (MQ=255)
                                                                                                                                 taGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGGTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                     >  2:1975164/1‑150 (MQ=255)
                                                                                                                                  aGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                    >  1:1609226/1‑150 (MQ=255)
                                                                                                                                    aTTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                  >  1:2830121/1‑150 (MQ=255)
                                                                                                                                    aTTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCg                  >  2:2054351/1‑150 (MQ=255)
                                                                                                                                      ttttCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGc                >  1:976377/1‑150 (MQ=255)
                                                                                                                                         tCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTc             <  2:2294515/150‑1 (MQ=255)
                                                                                                                                          caggcgCGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTATTGGTCAGCAGCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCa            >  2:1664123/6‑150 (MQ=255)
                                                                                                                                             gcgCGTTTAACCCAGCCGGG‑TAAACGCACCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCg         >  2:142105/3‑150 (MQ=255)
                                                                                                                                               gCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGt       <  1:2201434/150‑1 (MQ=255)
                                                                                                                                                 gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:885829/1‑150 (MQ=255)
                                                                                                                                                 gTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGa     >  1:885569/1‑150 (MQ=255)
                                                                                                                                                    t‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt  >  2:2346669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCATAACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTT‑ACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACT  >  NC_000913/1194497‑1194794

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