Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,327,216 A→G 7.2% D122G (GAT→GGT)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2160AG 7.2% 457.3 / 16.5 152D122G (GAT→GGT) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (63/78);  new base G (7/4);  total (70/82)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCA  >  NC_000913/1327067‑1327365
                                                                                                                                                     |                                                                                                                                                     
ttACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                                                                       >  1:450270/1‑150 (MQ=255)
 tACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGAt                                                                                                                                                      >  1:402606/1‑150 (MQ=255)
   cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                    >  1:901164/1‑150 (MQ=255)
   cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                    >  1:1300867/1‑150 (MQ=255)
   cAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATgg                                                                                                                                                    >  1:2033091/1‑150 (MQ=255)
    aaaCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGt                                                                                                                                                   <  2:2282018/150‑1 (MQ=255)
     aaCCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTg                                                                                                                                                  >  1:1136151/1‑150 (MQ=255)
           aaGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagc                                                                                                                                            <  1:905838/150‑4 (MQ=255)
           aaGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTgg                                                                                                                                            <  2:795309/150‑1 (MQ=255)
            aGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                                                                           <  1:582794/150‑1 (MQ=255)
            aGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGc                                                                                                                                           <  2:143815/150‑1 (MQ=255)
            aGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTGTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGGGAACTGGc                                                                                                                                           <  1:1031209/150‑1 (MQ=255)
             ggTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCa                                                                                                                                          >  1:656459/1‑150 (MQ=255)
              gTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCATCGGTGGTGAAcagcagc                                                                                                                                         <  1:384066/150‑7 (MQ=255)
                gAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAAc                                                                                                                                       <  1:2259790/150‑1 (MQ=255)
                  gTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCg                                                                                                                                     <  2:652053/150‑1 (MQ=255)
                   ttGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGt                                                                                                                                    >  1:1903162/1‑150 (MQ=255)
                     gtgtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGttt                                                                                                                                  <  1:664839/150‑1 (MQ=255)
                      tgtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTa                                                                                                                                 >  1:1752939/1‑150 (MQ=255)
                       gtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaa                                                                                                                                <  1:1434423/150‑1 (MQ=255)
                       gtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaa                                                                                                                                >  1:1763368/1‑150 (MQ=255)
                       gtACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGAGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTaa                                                                                                                                >  2:817874/1‑150 (MQ=255)
                         aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                                                                              <  2:760118/150‑1 (MQ=255)
                         aCCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATg                                                                                                                              <  1:1236787/150‑1 (MQ=255)
                           ccGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCa                                                                                                                            <  2:872321/150‑1 (MQ=255)
                            cGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGa                                                                                                                                                       >  1:652240/1‑122 (MQ=255)
                            cGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAc                                                                                                                           >  1:680417/1‑150 (MQ=255)
                              tAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAccc                                                                                                                         <  1:2152873/150‑1 (MQ=255)
                                   aTCCGGAAGGACGTCAAACGGTGTTTGACAGTCTGCCAAAACTGCGTGCAGCACGCTGGATTTCCGTGGGTCCCTAGGACGTTAATACCTGTGCTCTGCTGCTGTTCACAACCGATCGTGAACTGGCAAACCGTTTAATGCACCCAAgac                                                                                                                    >  2:1363097/1‑148 (MQ=255)
                                    tCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCg                                                                                                                   <  1:276925/150‑1 (MQ=255)
                                       ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                                <  2:367198/150‑1 (MQ=255)
                                       ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                                >  1:2083750/1‑150 (MQ=255)
                                       ggAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGa                                                                                                                <  1:636982/150‑1 (MQ=255)
                                         aaGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAg                                                                                                              <  1:1826174/150‑1 (MQ=255)
                                          aGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGt                                                                                                             <  2:1827836/150‑1 (MQ=255)
                                                   aaCGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGa                                                                                                    >  2:162602/1‑150 (MQ=255)
                                                      ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                                                                                 <  2:147817/150‑1 (MQ=255)
                                                      ggTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGGCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAata                                                                                                 <  2:2383569/150‑1 (MQ=255)
                                                         gtTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGc                                                                                              >  1:1829836/1‑150 (MQ=255)
                                                            tGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgt                                                                                           >  2:1068327/1‑150 (MQ=255)
                                                             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                          >  2:1004416/1‑150 (MQ=255)
                                                             gACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtg                                                                                          >  1:480968/1‑150 (MQ=255)
                                                               ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                        <  2:1134698/150‑1 (MQ=255)
                                                               ccGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCcgtgcg                                                                                        >  2:2233199/1‑150 (MQ=255)
                                                                  tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCgtgt                                                                                     <  2:1656127/150‑1 (MQ=255)
                                                                   cTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCTTGAAGTTGAACGTGAATATGCCGTGCGTGTa                                                                                    <  2:1071824/150‑1 (MQ=255)
                                                                    tGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAt                                                                                   >  1:327590/1‑150 (MQ=255)
                                                                     gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTAtt                                                                                  >  2:33977/1‑150 (MQ=255)
                                                                       cAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTg                                                                                <  1:237262/150‑1 (MQ=255)
                                                                        aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                                                                               <  2:1214905/150‑1 (MQ=255)
                                                                        aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                                                                               <  2:906672/150‑1 (MQ=255)
                                                                        aaaaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTgg                                                                               <  1:725331/150‑1 (MQ=255)
                                                                          aaCTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTc                                                                             <  2:1721041/150‑1 (MQ=255)
                                                                                gTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTg                                                                       <  1:2035631/150‑1 (MQ=255)
                                                                                  ggCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                                                                     >  1:1163010/1‑150 (MQ=255)
                                                                                     gcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgacgac                                                                  <  1:592975/150‑1 (MQ=255)
                                                                                                tGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                       <  1:778022/150‑1 (MQ=255)
                                                                                                 gCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt                                                      >  1:2282654/1‑150 (MQ=255)
                                                                                                    tcgggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                   >  1:635133/3‑150 (MQ=255)
                                                                                                    gTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGAt                                                   <  1:549509/150‑1 (MQ=255)
                                                                                                       ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                                <  2:1903162/150‑1 (MQ=255)
                                                                                                       ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                                >  1:1607685/1‑150 (MQ=255)
                                                                                                       ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                                <  2:2081213/150‑1 (MQ=255)
                                                                                                       ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTg                                                >  2:1498488/1‑150 (MQ=255)
                                                                                                        gTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGa                                               >  1:188297/1‑150 (MQ=255)
                                                                                                            cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTgac                                                                     <  1:2221615/124‑1 (MQ=255)
                                                                                                              tGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTg                                         >  1:2290749/1‑150 (MQ=255)
                                                                                                               ggACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTgg                                        <  1:1036018/150‑1 (MQ=255)
                                                                                                                 aCGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCg                                      <  2:1085678/150‑1 (MQ=255)
                                                                                                                 aCGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCg                                      <  2:347511/150‑1 (MQ=255)
                                                                                                                   gTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                    >  2:178473/1‑150 (MQ=255)
                                                                                                                    ttAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCAAAGCCGTTAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGAGAAACTGCGTGATTTGAGTCGTGGCGTGc                                   <  2:2033091/150‑1 (MQ=255)
                                                                                                                     tAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGGGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCa                                  <  1:220277/150‑1 (MQ=255)
                                                                                                                           gtgtgGTCTGCTGCTGGTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCAGCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGa                            <  2:1397840/149‑1 (MQ=255)
                                                                                                                           cTGTGGTCTGCTGCTGTTCACAACCGATGGTGAACTGGCAAACCGTTTAATGCCCCCAAGCCGTGAAGTTGAACGTGAATATGCCGGGCGTGTATTTGGTCAGGTTGACGCCGCGAAACAGCGTGATTTGAGTCGTGGCGTGCAGTTGGa                            >  2:2476694/1‑150 (MQ=255)
                                                                                                                            tgtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGaa                           >  2:143991/1‑150 (MQ=255)
                                                                                                                                gTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATTCACCCAAGCCGTGAAGTTCAACGAGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATg                       >  2:1361051/1‑150 (MQ=255)
                                                                                                                                gTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATg                       <  2:2148486/150‑1 (MQ=255)
                                                                                                                                     tgccagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCg                  >  2:1177637/6‑150 (MQ=255)
                                                                                                                                     ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCg                  >  1:576347/1‑150 (MQ=255)
                                                                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGGAGATGGTCCGGc                >  2:2304382/1‑150 (MQ=255)
                                                                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGc                >  2:446570/1‑150 (MQ=255)
                                                                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGc                >  2:11363/1‑150 (MQ=255)
                                                                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGc                >  2:1678194/1‑150 (MQ=255)
                                                                                                                                       gctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGc                >  2:322907/1‑150 (MQ=255)
                                                                                                                                        cagTTCACCATCGGTGGTGAACTGGCAGACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGGATTTGGTCAGGTTGCCGACGCGAACCTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCa               >  1:2486669/3‑150 (MQ=255)
                                                                                                                                        cagTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCa               >  1:676869/3‑150 (MQ=255)
                                                                                                                                        cagTTAACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATTCCGTGCGTGTATTTGGTCAGGTTGAAGACGCGAAACTGCGTGATTTGAGTCGTGGAGTGCAGTTGGAAGATGGTCCGGCa               >  2:2156769/3‑150 (MQ=255)
                                                                                                                                         agTTCACCATCGGTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAg              >  1:1525447/2‑150 (MQ=255)
                                                                                                                                          gTTCACGACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGCCGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGGTGGAAGATGGTCCGGCAGc             >  1:2126808/1‑150 (MQ=255)
                                                                                                                                            tCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCtt           >  1:222469/1‑150 (MQ=255)
                                                                                                                                              accaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCtttt         <  2:1188142/150‑1 (MQ=255)
                                                                                                                                               ccaccGATGGTTAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa        <  2:267480/150‑1 (MQ=255)
                                                                                                                                               ccaccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTa        <  2:276038/150‑1 (MQ=255)
                                                                                                                                                  ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaaa     <  2:676368/150‑1 (MQ=255)
                                                                                                                                                   cGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAAc    <  1:817874/150‑1 (MQ=255)
                                                                                                                                                   cGATGGTGAACTGGCAAACCGTTTAAAGCACCCAAGCCGTGAAGTTGAACTTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACAGCGTGATTTTAGTCGTCGCGTTCAGTTGGAACATGGTCAGGCAGCTTTTAAAAc    >  2:56616/1‑150 (MQ=255)
                                                                                                                                                     aTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCa  <  2:1752939/150‑1 (MQ=255)
                                                                                                                                                     aTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCa  <  1:2233199/150‑1 (MQ=255)
                                                                                                                                                     aTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGGGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCa  <  1:184337/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTACAAACCGGAAGGTGAGTTGTGTACCCGTAACGATCCGGAAGGACGTCCAACGGTGTTTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCA  >  NC_000913/1327067‑1327365

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