Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,653 T→C 6.5% Q100R (CAG→CGG)  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,6530TC 6.5% 317.2 / 9.7 108Q100R (CAG→CGG) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base T (55/46);  new base C (2/5);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-01

Only 100 of 110 total aligned reads displayed.

ACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  NC_000913/1194504‑1194802
                                                                                                                                                      |                                                                                                                                                     
aCCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCt                                                                                                                                                       >  1:1813450/1‑150 (MQ=255)
 ccACGCGGGTTCGTGAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGACATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTg                                                                                                                                                      <  2:255526/150‑1 (MQ=255)
  cACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTgg                                                                                                                                                     >  2:131185/1‑150 (MQ=255)
  cACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTgg                                                                                                                                                     >  2:64146/1‑150 (MQ=255)
  cACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTgg                                                                                                                                                     >  1:449009/1‑150 (MQ=255)
  cACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTgg                                                                                                                                                     >  1:314075/1‑150 (MQ=255)
  cACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTgg                                                                                                                                                     >  2:1554326/1‑150 (MQ=255)
      cgGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTTTGTGGGAATCCCTTGCACCTGCACATTATAGATTTTACCGGTGCGTTTACCCGGCTGGGTTa                                                                                                                                                 >  2:87373/1‑150 (MQ=255)
        ggTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTaaa                                                                                                                                               >  1:748727/1‑150 (MQ=255)
        ggTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTaaa                                                                                                                                               <  2:1554841/150‑1 (MQ=255)
           tCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAAcgc                                                                                                                                            >  1:244704/1‑150 (MQ=255)
             gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACgcgc                                                                                                                                          >  1:216357/1‑150 (MQ=255)
             gTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACgcgc                                                                                                                                          >  2:1250854/1‑150 (MQ=255)
              tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAcc                                                                                                                                         <  2:692161/150‑1 (MQ=255)
              tCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGcc                                                                                                                                         >  2:1235947/1‑150 (MQ=255)
              gcAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGcc                                                                                                                                         <  1:386517/149‑1 (MQ=255)
               cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTACTGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                        <  1:1286097/150‑2 (MQ=255)
               cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                        <  1:244318/150‑2 (MQ=255)
               cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGCCTGGGTTAAACGCGCCt                                                                                                                                        >  1:667474/1‑150 (MQ=255)
               cAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccg                                                                                                                                        <  2:1426579/150‑2 (MQ=255)
                aaCCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGGCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGCGCGTTTAACCCAGCCGGGTAAACGCAccgg                                                                                                                                       <  1:1187161/150‑3 (MQ=255)
                      cTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCCCCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAAcgc                                                                                                                                 >  1:1489625/1‑150 (MQ=255)
                       tcagcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCGATTGCGCAAGGCTTCAAGTCCGTCTTGTGAGGGAATAACTTCCACCTGCACATAAAAGATTTTTCCGGTGCGTATACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                >  1:359969/4‑150 (MQ=255)
                       tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                <  1:1428191/150‑1 (MQ=255)
                       tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                <  1:805230/150‑1 (MQ=255)
                       tCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCAACTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAAcgag                                                                                                                                >  1:975379/1‑148 (MQ=255)
                       tCCGCGCCGGCGGGCAGGGTAGGGCCAACATTAAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATCCCTTCCACCTGCACAGAATAGATATTTCCGGTGCGATTACCCGGCTGGGTTAAACGCGCCTGCAAcgcg                                                                                                                                >  1:226152/1‑150 (MQ=255)
                          gcgcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCg                                                                                                                             <  2:1187566/150‑1 (MQ=255)
                            gcCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGATTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGtt                                                                                                                           >  1:48524/1‑150 (MQ=255)
                               ggcggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTAtt                                                                                                                        >  2:635720/1‑150 (MQ=255)
                                 cggGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTgg                                                                                                                      <  2:574724/150‑1 (MQ=255)
                                   ggCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTc                                                                                                                    <  2:748727/150‑1 (MQ=255)
                                    gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCTCAAGGCTTCATGTGCGTCTAGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                   <  1:96104/150‑1 (MQ=255)
                                    gCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCa                                                                                                                   <  2:1098572/150‑1 (MQ=255)
                                     cAGGGTAGGGCCCTCATTTAAGGTGAGGCCATTGCGCAAGGCTTCAAGGGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAg                                                                                                                  <  2:266552/150‑1 (MQ=255)
                                      aGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGc                                                                                                                 <  2:394882/150‑1 (MQ=255)
                                        ggTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCAc                                                                                                               <  2:31631/150‑1 (MQ=255)
                                              gCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCaa                                                                                                         >  2:673429/1‑150 (MQ=255)
                                                 atcatTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAAccc                                                                                                      <  1:221154/150‑1 (MQ=255)
                                                     tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCACCGCGCCGTTATTGGTCCGCACCAGCAACCATTCg                                                                                                  >  2:1244693/1‑150 (MQ=255)
                                                     tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                  <  2:1328432/150‑1 (MQ=255)
                                                     tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCAACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGACGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                  >  1:65934/1‑150 (MQ=255)
                                                     tttAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGGCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCg                                                                                                  <  1:1721341/150‑1 (MQ=255)
                                                         aGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTAt                                                                                              >  1:1620615/1‑150 (MQ=255)
                                                           gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            >  1:1624676/1‑150 (MQ=255)
                                                           gTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATcg                                                                                            >  1:1449523/1‑150 (MQ=255)
                                                               cGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGt                                                                                        >  2:1343686/1‑150 (MQ=255)
                                                                gCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTc                                                                                       <  2:526614/150‑1 (MQ=255)
                                                                 ccATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCa                                                                                      >  2:1771198/1‑150 (MQ=255)
                                                                      gcgcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCg                                                                                 <  1:107272/150‑1 (MQ=255)
                                                                        gcAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGAc                                                                               <  2:1526253/150‑1 (MQ=255)
                                                                          aaGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACct                                                                             <  1:128906/150‑1 (MQ=255)
                                                                             gCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgct                                                                          <  2:882606/150‑1 (MQ=255)
                                                                              cTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACctgctg                                                                         >  1:1366160/1‑150 (MQ=255)
                                                                              cTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGAAACCCTTCGCTATCGCGGTCAAGGCGACctgctg                                                                         >  1:1549262/1‑150 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        >  2:1436445/1‑150 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        <  2:58738/150‑1 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        <  1:1743181/150‑1 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        >  2:1595795/1‑150 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        <  2:244704/150‑1 (MQ=255)
                                                                               ttCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCtgctgc                                                                        <  1:1171377/150‑1 (MQ=255)
                                                                                tCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCa                                                                       <  2:1029099/150‑1 (MQ=255)
                                                                                tCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGGCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCa                                                                       <  2:1021908/150‑1 (MQ=255)
                                                                                  aaGTGCGTCTTGTGTGGTAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATa                                                                     <  2:1386716/150‑1 (MQ=255)
                                                                                      gCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCTCGGTCAAGGCGACCTGCTGCATAAacc                                                                 >  2:1660952/1‑149 (MQ=255)
                                                                                         tCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACAccc                                                              <  2:1461994/150‑1 (MQ=255)
                                                                                             gtgtgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGaa                                                          >  2:1248854/1‑150 (MQ=255)
                                                                                                tgGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGGTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCg                                                       <  1:281089/150‑1 (MQ=255)
                                                                                                tgGGAATACCTTCCACCTGCACATAATAGAGTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGCCGACCTGCTGCATAAACACCCTGAACCg                                                       >  1:295324/1‑150 (MQ=255)
                                                                                                  ggAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCggg                                                     >  2:821859/1‑150 (MQ=255)
                                                                                                   gAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGa                                                    <  1:723348/150‑1 (MQ=255)
                                                                                                     aTACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATg                                                  <  2:216357/150‑1 (MQ=255)
                                                                                                          ttCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTc                                             >  1:452218/1‑150 (MQ=255)
                                                                                                               ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  1:709677/1‑150 (MQ=255)
                                                                                                               ccTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        >  1:1067460/1‑150 (MQ=255)
                                                                                                               ccTGCACATAATAGATTTTTCCGGTGCGGTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTa                                        <  2:1099926/150‑1 (MQ=255)
                                                                                                                 tGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCCCCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAATCACCCTGCACCGGGATGAATTCTTTTAAt                                      >  1:1253292/1‑150 (MQ=255)
                                                                                                                  gCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATg                                     <  2:1509319/150‑1 (MQ=255)
                                                                                                                      ataataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCAGTTATTGGTCAGCACCAGCAACCCTTCGCTATAGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGtttt                                 >  1:77164/1‑150 (MQ=255)
                                                                                                                       taataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCCGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACACTGAACCGGGATGAATTCTTTTAATGTTTTg                                <  1:1122550/150‑1 (MQ=255)
                                                                                                                         ataGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCg                              >  1:1590783/1‑150 (MQ=255)
                                                                                                                          taGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCACCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGt                             >  2:1133783/1‑150 (MQ=255)
                                                                                                                           aGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            >  2:708528/1‑150 (MQ=255)
                                                                                                                           aGATTTTTCCGGAGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTc                            <  1:936839/150‑1 (MQ=255)
                                                                                                                            gATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTcc                           >  1:552245/1‑150 (MQ=255)
                                                                                                                              tttttCAGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCgg                         >  1:789113/1‑150 (MQ=255)
                                                                                                                                   ccGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGATTTCTTTTATTGTTTTGCGTCCGGCTTCa                    >  1:427995/1‑150 (MQ=255)
                                                                                                                                   ccGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCa                    >  1:1612355/1‑150 (MQ=255)
                                                                                                                                   caggcgCGTTTAACCCAGCCGGGTAAACGCACCTGCAACGCGCCGTTTTTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAGTCTTTGAATGTTTTGCGTCCGGCTTCa                    >  1:1653449/6‑150 (MQ=255)
                                                                                                                                    cGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGCACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCAt                   >  1:1107148/1‑150 (MQ=255)
                                                                                                                                      gTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGACAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTCCATCg                 >  1:509049/1‑150 (MQ=255)
                                                                                                                                         cGTTTAACCCAGCCGGGTAAACGCGCCTGCACCGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGCACCGGGATGAATTCTTTTTATGTTTTGCGTCCGGCTTCATCGGTg              >  1:724781/1‑150 (MQ=255)
                                                                                                                                         cGTTTACCCGGCTGGGTTTAACGCGCCTGCAACGCGCCGTTATGGGGCTGCACCAGCACCGCTGCGCTATCTCGGTCACGGCGACCTGCTTCATAAACTCCCTTAACCGTGATGAATTCTTTTAATGTTTTGCTTCCGGCTTCATCGGTg              <  1:1114632/150‑1 (MQ=255)
                                                                                                                                         cGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAAATCTTTTAATGTTTTGCGTCCGGCTTCATCGGTg              <  1:64146/150‑1 (MQ=255)
                                                                                                                                           tttACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGACCCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGaa            >  1:926488/1‑150 (MQ=255)
                                                                                                                                             tACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACt          <  1:1554326/150‑1 (MQ=255)
                                                                                                                                                cccGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGc        >  2:19290/1‑150 (MQ=255)
                                                                                                                                                  cGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCgg      <  2:1269136/150‑1 (MQ=255)
                                                                                                                                                     cTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGgaag   <  1:1786819/150‑1 (MQ=255)
                                                                                                                                                      tGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGTGgaaga  <  2:1449523/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ACCACGCGGGTTCGTCAACCAGCTCCGCGCCGGCGGGCAGGGTAGGGCCATCATTTAAGGTTACGCCATTGCGCAAGGCTTCAAGTGCGTCTTGTGTGGGAATACCTTCCACCTGCACATAATAGATTTTTCCGGTGCGTTTACCCGGCTGGGTTAAACGCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGA  >  NC_000913/1194504‑1194802

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: