Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,457,392 C→T E158E (GAG→GAA yjgA ← putative ribosome biogenesis factor YjgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,3920CT96.2% 165.6 / NA 54E158E (GAG→GAAyjgAputative ribosome biogenesis factor YjgA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (33/18);  minor base G (0/2);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCG  >  NC_000913/4457149‑4457628
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
cgcgCAGTAAGTCGGCCTTATCCAGGTATAGTAGGCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCATTTGAGCATAGTGGCTGAGGCGACTGTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACACACATCAACAACGAAAAGTACCCTGCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTctttgtctttct                                                                                                                                                                                                                                       <  2:397090/251‑1 (MQ=255)
        aaTTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTg                                                                                                                                                                                                                               <  1:17151/251‑1 (MQ=255)
             ggCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGAt                                                                                                                                                                                                                          >  1:472592/1‑251 (MQ=255)
                cTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCggcac                                                                                                                                                                                                                         >  2:114086/1‑246 (MQ=255)
                   aGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGt                                                                                                                                                                                                                    <  1:441794/251‑1 (MQ=255)
                             aGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGc                                                                                                                                                                                                          >  1:492827/1‑251 (MQ=255)
                              gTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCt                                                                                                                                                                                                         <  1:340243/251‑1 (MQ=255)
                                      aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCg                                                                                                                                                                                                 <  1:422750/251‑1 (MQ=255)
                                         ccGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGATTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGcc                                                                                                                                                                                              >  1:327767/1‑250 (MQ=255)
                                          cGGATAAGGCGTTCACACCACAACAAGCATTTGATCATAACGCCTGGTACGCCACTTACTGCGTCCTATCAGACCTCCAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCCTCGTTCACCGCTAAATCGCACAGATACTGGACAATCTGGCGAGCGGCTTTCGGCGGCATATTCACTTCTTTTTATTTCTTCGCGATGCGGATCAGAGTAAGCAGATGCTGACTATCGGCAt                                                                                                                                                                                             >  2:320557/1‑251 (MQ=255)
                                               acggcggtcaagacgcatcggtgagttGAGCAGAGGGCCTGATGCGACTCTTGATGCGGCTTATCAGGCCTACAGGTGCACCGTATCCGGCAACCCATACTCAGACCAACGACGAAAATTACCCTTCGTGCTCCGCTAACTCGCGCAGATGCTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCGTCTTTGTCTCTCTGCTCGAAGTGGATCAGTGTACGCAGTTTCTGACGGTCGCCATAggg                                                                                                                                                                                         <  2:70631/225‑1 (MQ=255)
                                                   cGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCcac                                                                                                                                                                                     >  1:220852/1‑250 (MQ=255)
                                                    gTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAg                                                                                                                                                                                   <  2:482466/251‑1 (MQ=255)
                                                       cACGCCGCATCCGGCATTTGAGCATAGTGCATGATGCGACGCTTGATGCGTCTTATCAGGCCTACCGGTGAACCGTATCCGGCAAACCATACTCACATCCACAACGAAAATTACCCTTCGATCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTACCTTCATTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGAtt                                                                                                                                                                                >  1:376875/1‑251 (MQ=255)
                                                             gcATCCGGCATGTGAGCATAGTGACTGATGCGACTCTGGATGCGTCTGAGCAGGCCTGCAGGTGCACCGTAGCCGGCAAACCATACTCACATCAACAACGAAAATTCTCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTGGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAAttc                                                                                                                                                                          <  2:46592/251‑3 (MQ=255)
                                                                         tGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGc                                                                                                                                                              >  1:242009/1‑251 (MQ=255)
                                                                            gCATAGTGCCTGATGCGACGCTTGATGCGACTTATCAGGCCTACAGGTGCACCGTCTCCGGAAAACCATACTCACATAAACAACGAAAATTACCCTTCGTTCTCCGCTACCTAGCGCAGATACTGGAACATCTGGCGTGCGGTTTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAATTGCTGACGATCGGCATCGGGCAACAGATTTAATACCTCAGCGATGGCGtc                                                                                                                                                           >  2:320881/1‑251 (MQ=255)
                                                                                    ccTGATACGACGCTTGATGAGTCTTAACAGGCCTACAGGTGCACCGTCTCCGGCCAACCATACTCACCTCACCAAAGAAAATTGCCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGACAATCTGGCGTGCGGATTTAGGAGGCTTATTCACTTCTTTTTCTTTCTTCGAGTTGCGGATCAGAGTACGCAGTTGCTGCAGATCGGCCTCGGGCCACAGTTTTAATACCACAGCGATGGCGTCATCAccct                                                                                                                                                   >  2:69827/1‑249 (MQ=255)
                                                                                           cGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTgatcgat                                                                                                                                            <  1:553786/251‑1 (MQ=255)
                                                                                             atgCTTGATGAGTCTTATCAGGTCTACAGGTTCACATTATCCGGCAAACCATACCCCCATCAACAACGAAAACTACCCTTCTTTCTCCGCTAACGCGTGCCGTTACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCGTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGGTGCTTACGATCGGCAGCGGGCCACAGATTTAATACCTCAGCGATGGCGTCTTCACCTTGATAGATCa                                                                                                                                          <  2:421717/249‑1 (MQ=255)
                                                                                                          cTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGCAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCgcgc                                                                                                                               >  2:285462/1‑249 (MQ=255)
                                                                                                                 ggCCTACAGGTGCACCGTATCCGGCAAACCATACTCACCTCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTcc                                                                                                                      >  1:219227/1‑251 (MQ=255)
                                                                                                                           tGCACCGTATCCGGCAAACCATACTCACATCAACACCGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATgg                                                                                                             >  1:119830/1‑250 (MQ=255)
                                                                                                                           ggcaacttatcctgaaaactttACTCACAACAGAAACGAAAATTACCCGTCGGTCTCCGCTAACTCGCGCAGAGACTGGAAAATCTGGCTTGCGGATTGCGGCGGCTTATTACCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGTCCTCAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGTGCATGTTTTCCTGTTTATGGa                                                                                                            <  2:557324/230‑1 (MQ=255)
                                                                                                                                 gTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGca                                                                                                      <  1:285462/251‑1 (MQ=255)
                                                                                                                                 gTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCAGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGcc                                                                                                      >  1:240505/1‑250 (MQ=255)
                                                                                                                                   aTCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGGACAGcacc                                                                                                    >  2:166786/1‑251 (MQ=255)
                                                                                                                                              cATACCCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGtt                                                                                         <  1:43312/251‑1 (MQ=255)
                                                                                                                                                 aCTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCCCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCTGCATCGGGCCACAGCTTTAATACCTCAGCGATGTCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACCGCACCACCTGCTGGtttgt                                                                                      >  2:105464/1‑248 (MQ=255)
                                                                                                                                                      caTCAACAACGAAAATTTCCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCTGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACTATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGt                                                                                 <  1:69827/251‑1 (MQ=255)
                                                                                                                                                       aTCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGtt                                                                                <  1:97700/251‑1 (MQ=255)
                                                                                                                                                           acaacGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGCTCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTCACGttttttt                                                                            >  2:138898/1‑250 (MQ=255)
                                                                                                                                                                   aaTTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGCTCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGAGACGGGTTTTCAGCTTATc                                                                    >  1:341970/1‑251 (MQ=255)
                                                                                                                                                                               ttCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACg                                                        >  1:461927/1‑251 (MQ=255)
                                                                                                                                                                                tctcCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACg                                                        <  1:545697/250‑1 (MQ=255)
                                                                                                                                                                                   ccGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGCTCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTACTACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCTCGCAGGTTTTCCAGTTTATGGGCCAGCACCACCTGCTGGTTGTTACGGTTTTTCAgatttcccaccgccggccgact                                                     >  2:254708/1‑229 (MQ=255)
                                                                                                                                                                                     gCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAgg                                                  >  2:263959/1‑251 (MQ=255)
                                                                                                                                                                                      cTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCCGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGc                                                 >  1:272905/1‑251 (MQ=255)
                                                                                                                                                                                       tAACTCGCGCAGATACTGGAAAATCTGGCGTGCAGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCGTCGGGCCACAGATTTAATACCTCAGCCAAGGCGTCATCACCTTGATCGATCAGACGATTGCGCAGGTTTTCCAGTTTATGGAACAGACCCCCCTGCTGGTAGTGACGGTTTTTCAGCTTATCCCGAGCCTGACGACTAGGct                                                >  2:562645/1‑251 (MQ=255)
                                                                                                                                                                                            cgcgCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGCTTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCCGCGCCTGACGAATAGGCTTTAcc                                           >  2:144098/1‑250 (MQ=255)
                                                                                                                                                                                            cgcgCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTa                                             <  1:505935/249‑1 (MQ=255)
                                                                                                                                                                                              cgccgatgcTGGAAAATCTGGCGTGCGGATTTCGGCCGCTTACTCCCGTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGGACGCAGTTGCTGACTATCGGCATCGGGCCACAGATTTAATTCCTCAGCGATGGCGTCATCACCTGGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATTGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGTCTCTACGTc                                         <  2:242009/243‑1 (MQ=255)
                                                                                                                                                                                                taGATACAGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCATTCTTCGCGTTGCGGGTCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGCTGGCGTCATCAACTTGATCGATCCGACGATCGCGCAGGTTTTCCAGTTTATGGAACCGAACCCACTGCTGGTTGTGACGGTTTTTCCACTTATCACGCGCCTGACAAATAGGCTCTACGTcgc                                       >  2:429740/2‑251 (MQ=255)
                                                                                                                                                                                                 aGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGATCAGCACCACCTGCTGATTGTGACGGTTTTTCAGCTTATCCAGCGCCGGACGAAAAGGCTCTACGTcacc                                      >  2:140442/1‑248 (MQ=255)
                                                                                                                                                                                                 aGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTcgcg                                      <  1:430493/251‑1 (MQ=255)
                                                                                                                                                                                                   aTACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCt                                    >  1:211573/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGTGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGCCGAATAGGCTCTACGTCGCGCTGgcgc                               >  2:220463/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGgcg                                <  1:114086/250‑1 (MQ=255)
                                                                                                                                                                                                          aaaaTCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGgcgc                               >  1:110701/1‑249 (MQ=255)
                                                                                                                                                                                                                   cGTGCGGAGTCCGGCGGCTTATTCCCTTGTTTGTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCTCCTTGAACGATCAGACTATCGCTCAGTTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATTGGCTCTACGTCGCGCTGGCGTCGCCTCTTAcc                    <  2:215416/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCg                   >  1:462824/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGc             >  1:324711/1‑251 (MQ=255)
                                                                                                                                                                                                                                gcggcTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGCTTATGGAACAGCACCGCCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAAc       >  1:69336/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggcTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTg     >  1:102149/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ttATTCCCTTCTTTTTCTTTCTTAGAGTTGCGGATCCGAGTACGCAGTTACTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCAACTTTATCGATCACACGCTCGCGCAGGTGTTACAGTTTATGGAACAGCACAACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGCCGAATAGGCTCTAAGTCGCGCTCGCTCAGCATCTTACCGAAGAGTTGCacatggcg  >  2:491483/1‑244 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
CGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCG  >  NC_000913/4457149‑4457628

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