Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,532,537 G→C 7.1% A216G (GCC→GGC)  yjhP ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,5370GC 7.1% 287.2 / 11.9 85A216G (GCC→GGC) yjhPputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/36);  new base C (3/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-02

CTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGG  >  NC_000913/4532317‑4532797
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
cTCCTTATGGTGACATCGTGATTAATAAACATAGTGCATATGTATTGAGCCGGGGATACAGCAGGTCAAAGCAAAGGTAAAGCATCGCTATTTGCAGCCAGAAATACGGGGAACCTGTTGTTACCGGGGGATGAGTGCGAACACGCCCCAGCCAAGGCAGTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                                                                                                                                                <  4:268182/275‑1 (MQ=255)
            aCATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATAGTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                                                                                                                                                                     <  2:308848/274‑1 (MQ=255)
             cATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCCTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATa                                                                                                                                                                                                      >  2:70665/1‑272 (MQ=255)
                    aTTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCTCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCCAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCCCATGGCCAGCCATTTTTCGTCTTCATACCTGCCCCa                                                                                                                                                                                             >  4:172529/1‑274 (MQ=255)
                             cagaCTGAATAGGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTgg                                                                                                                                                                                   <  2:108572/272‑1 (MQ=255)
                                      tatGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGAGATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt                                                                                                                                                                          <  2:172531/275‑1 (MQ=255)
                                      gaggTAGTGTGCAGGTGAGACCAGAGATCAACGCCAAGGTAAAGCCGGGCAATTGCCACCCATGTATCACGGGAACCTGTGGTTGCCGTGCGAGGAGTGCGAACACGCCCCAGCGAACGCGTGCCCGCGCGTCTGTCACGTAGCGTCTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCCGGATTCGCTTCCAGCCAGCGGCGCATGGGCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt                                                                                                                                                                          <  4:227453/272‑1 (MQ=255)
                                              ttgCAGGTGAAACATCAGATCAAAGCAAAGGTAAAGCCTGGATATTTCCAGCCAGAAAGACAGGAAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTc                                                                                                                                                                  <  2:497879/274‑1 (MQ=255)
                                              gtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCa                                                                                                                                                                                                                                                              >  2:235407/1‑183 (MQ=255)
                                              gtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGACTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATGTTGAGGCTTACTACCGGTCCCAGCCTTCCTGGTCTGCCAGTCCCCttt                                                                                                                                                                   >  2:302616/1‑274 (MQ=255)
                                                 cAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAg                                                                                                                                                                                                                                                       >  4:251053/1‑187 (MQ=255)
                                                  aGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                                                                                                                                          <  4:235403/183‑1 (MQ=255)
                                                  aGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATGTCCAGCCAGAAATACCGGGACCCTGTTGTTACCGGGCGGTGAGTGCGAACACGCCCAAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTgg                                                                                                                                                                                   <  2:167220/254‑1 (MQ=255)
                                                  aGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGACATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTcc                                                                                                                                                                                               <  2:173376/242‑1 (MQ=255)
                                                   ggTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                                                                                                                                                                                      <  1:14839/202‑1 (MQ=255)
                                                    gTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                                                                                                                                                                                            >  1:121966/1‑243 (MQ=255)
                                                    gTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCATTCCCCTTTccacca                                                                                                                                                             >  2:147457/1‑274 (MQ=255)
                                                    gTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTCTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTCCCATTTccacct                                                                                                                                                             >  2:227457/1‑273 (MQ=255)
                                                     tGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAtt                                                                                                                                                                                                                                                   <  2:251057/187‑1 (MQ=255)
                                                         aCATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTCGGTCACGTAGCGTTTCGGCGCGATGTTGAGTTCGGCCCTGACTTCGGCGGCGAGGTCATCGTACGGATTCGCTTCCAGCCAGCGGCGCATGTTCAGCCCTTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCCGTACCATTTCCACCACGTcct                                                                                                                                                       >  1:91002/1‑273 (MQ=255)
                                                              aGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTTCCAGCCTTCCTGTTCTTCCAGTACCATTCCCACCACGTCATAGcc                                                                                                                                                   >  2:54106/1‑274 (MQ=255)
                                                                aTCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCga                                                                                                                                                 >  2:263655/1‑274 (MQ=255)
                                                                        aaagGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcggcgg                                                                                                                                                                                                                                                                      >  8:156158/1‑149 (MQ=255)
                                                                             tAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAg                                                                                                                                                                                                                                                                 >  6:402238/1‑149 (MQ=255)
                                                                             tAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAg                                                                                                                                                                                                                                                                 >  8:39755/1‑149 (MQ=255)
                                                                              aaaGCCTGGATATTTCCAGCCAGAAATACCGGGAACCGGTTGTTACCGGGCGTAGAGGGCGACCAAGCCCCAGCCAAAGAATCCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTCACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCt                                                                                                                                   <  2:289529/274‑1 (MQ=255)
                                                                               aaGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGCACACGCCCCACACAAAGCATTCCCGCGCGTAGGTCACGTAGCGTCTCTGCGCCATGTTTAGTTCGGCCCTGACTTCGACTGACAAGTCCTCCTCAGACTTCGCTTCCAGCCAGCGGCTCATGCTCAGCCATTTTGCCGCTTCATCACGGTCCCGGCCTTCCTGGTCTACCAGTGTCATTTCCACTACGTCATAGCCGAGAtcgacaagcgccg                                                                                                                                  >  4:492510/1‑264 (MQ=255)
                                                                                 gCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                                                                                                                                                                                             <  5:437912/149‑2 (MQ=255)
                                                                                 gCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGCCGCCGAAGTCag                                                                                                                                                                                                                                                             <  7:39755/149‑2 (MQ=255)
                                                                                   tgcGATATTTCTGGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCAAGCCAAAGCATTCCCGGGGGTAGGTCACGTAGCGTTTCGGCGCGATGTTGAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCt                                                                                                                                   <  2:492515/267‑1 (MQ=255)
                                                                                   cTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                                                                                                                                                                                           >  5:372034/1‑149 (MQ=255)
                                                                                    tCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                                                                                                                                                                          >  7:19217/1‑149 (MQ=255)
                                                                                      gATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                                                                                                                                                                        >  7:163429/1‑149 (MQ=255)
                                                                                      gATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAAGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCCGTACCATTTCCACCACGTCATAGCCGAGAGCGTCGCACGCTCCGACAAGt                                                                                                                          >  2:199895/1‑275 (MQ=255)
                                                                                                        tACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                                                                                                                                                                                      >  6:332975/1‑149 (MQ=255)
                                                                                                                  cTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                                                                                                                                                                            >  5:114866/1‑149 (MQ=255)
                                                                                                                          aCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                                                                                                                                                                                              <  6:240367/139‑1 (MQ=255)
                                                                                                                          aCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAt                                                                                                                                                                                                                              >  5:240367/1‑139 (MQ=255)
                                                                                                                                         cGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAgg                                                                                                           >  2:27799/1‑239 (MQ=255)
                                                                                                                                             cacGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGt                                                                                                                                                                                                 <  8:163429/149‑1 (MQ=255)
                                                                                                                                             cacGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGt                                                                                                                                                                                                 >  7:431255/1‑149 (MQ=255)
                                                                                                                                               cGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCCCGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCTAGCTGACCCCGCAGGCCTGGCCTTTttct                                                                 >  2:242973/1‑275 (MQ=255)
                                                                                                                                                  cccAGCCAAAGCATTCCCGCGCGTCGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCTAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTTGCAGCGTCAGGAAATCGCTTGTCGAGCTGACACCGAAGGCCGTGGCTATttcgtgt                                                              >  2:101283/1‑271 (MQ=255)
                                                                                                                                                      gCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCtt                                                                                                                                                                                        <  6:126472/149‑1 (MQ=255)
                                                                                                                                                         aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTGGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                                                                                                                                                                     >  7:404135/1‑133 (MQ=255)
                                                                                                                                                         aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                                                                                                                                                                     <  8:404135/133‑1 (MQ=255)
                                                                                                                                                         aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTc                                                                                                                                                                                       >  5:307938/1‑147 (MQ=255)
                                                                                                                                                         aaaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTc                                                                                                                                                                                       <  6:307938/147‑1 (MQ=255)
                                                                                                                                                          aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGACCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                                                                                                                                            >  7:209031/1‑77 (MQ=255)
                                                                                                                                                          aaGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGACCCTGACTTCGGCGGCGAAGTCATc                                                                                                                                                                                                                                                            <  8:209031/77‑1 (MQ=255)
                                                                                                                                                            gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGa                                                                                                                                                                                                                                                     >  4:11521/1‑82 (MQ=255)
                                                                                                                                                            gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                                                                                                                                                                  <  8:19217/149‑1 (MQ=255)
                                                                                                                                                            gCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                                                                                                                                                                  >  6:66171/1‑149 (MQ=255)
                                                                                                                                                             cATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTc                                                                                                                                                                                 <  7:295539/149‑1 (MQ=255)
                                                                                                                                                              aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                                                                                                                                                <  5:332975/149‑1 (MQ=255)
                                                                                                                                                              aTTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                                                                                                                                                                                >  5:87227/1‑149 (MQ=255)
                                                                                                                                                                tCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCg                                                                                                                                                                                                                                                 <  2:11521/82‑1 (MQ=255)
                                                                                                                                                                tCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGc                                                                                                                                                                              <  5:402238/149‑1 (MQ=255)
                                                                                                                                                                 cccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc                                                                                                                                                                             <  5:128388/149‑1 (MQ=255)
                                                                                                                                                                  ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                                                                                                                                               <  5:477477/114‑1 (MQ=255)
                                                                                                                                                                  ccGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                                                                                                                                                               >  6:477477/1‑114 (MQ=255)
                                                                                                                                                                              cACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTcca                                                                                                                                                                <  6:31848/149‑1 (MQ=255)
                                                                                                                                                                              accgtatgttttCGTGGCGATTTGTGGTGNNGCCCTGACTTGGGCGGGCAAGTCACTCTCGNGATTCGCATTCAGCCAGCGTCGGTTGGTCGGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCGGGTCTGCCNGTACCATTTCCCCCACGTCAGAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGtt                                  <  3:62098/266‑1 (MQ=255)
                                                                                                                                                                               aCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccac                                                                                                                                                               <  8:431255/149‑1 (MQ=255)
                                                                                                                                                                                cGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccacca                                                                                                                                                             >  2:442878/1‑150 (MQ=255)
                                                                                                                                                                                 gTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTc                                                                                                                                                                                                >  2:545904/1‑114 (MQ=255)
                                                                                                                                                                                    gCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTc                                                                                                                                                         <  4:442873/150‑1 (MQ=255)
                                                                                                                                                                                     cGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                                                                                                                                                                                            <  4:545899/114‑1 (MQ=255)
                                                                                                                                                                                       tttCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCTAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAt                                                                                                                                                       <  6:372034/149‑1 (MQ=255)
                                                                                                                                                                                         tCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCGGCCGCCCGTAg                                      >  1:62098/1‑260 (MQ=255)
                                                                                                                                                                                            tcgcgATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGa                                                                                                                                                                                                                                                     >  7:270536/2‑50 (MQ=255)
                                                                                                                                                                                            gcgcGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGa                                                                                                                                                                                                                                                     <  8:270536/50‑1 (MQ=255)
                                                                                                                                                                                                   gTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgt                                                                                                                                           >  6:417304/1‑149 (MQ=255)
                                                                                                                                                                                                    tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtc                                                                                                                                          >  7:117004/1‑149 (MQ=255)
                                                                                                                                                                                                    tttAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGAtcgtc                                                                                                                                          <  6:114866/149‑1 (MQ=255)
                                                                                                                                                                                                           tgacgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCt                                                                                                                                   >  6:125010/7‑149 (MQ=255)
                                                                                                                                                                                                            ggCCCTGACTTCGGCGGCGAAGTCATCGTCCGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAACTCCTAGCCGCTTCAGGCAATCCCTTGTCCAGCCGACCCCGCCGGCCTGCGCCATTTCTTTCGTTCCTGGTAGCTGACGCCCGTATGCATCGCCGCTAACCCTGATCCTCCCCGGCTTc     >  4:309256/1‑274 (MQ=255)
                                                                                                                                                                                                             acgatGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTcc                                                                                                                                 <  5:125010/145‑1 (MQ=255)
                                                                                                                                                                                                              cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                                                                                                                                                                                                    >  2:100943/1‑65 (MQ=255)
                                                                                                                                                                                                              cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTTCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCCACAAGTCCTTGCAGGGTCAGGCAATCGCCTTTCGCACTTACGCCGCCGTCCTGGGCTATCTCTTTTGTTTCTGGTAGCTCACTCCAGCAGGCTTCCCCGACAGCCATGCTccaccccatcttcag   >  4:118983/1‑261 (MQ=255)
                                                                                                                                                                                                              cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGAAGCGTCAGGAAATCGCTTGTAGAGCTGACGCCGCAGGCCTGGGCTCTTTCTTCTGTTGGTGTTAGCTGACGCCAGTATGGGTGGCCGAGAAGCAAGATcacgcaggggttgcgg  >  1:54406/1‑264 (MQ=255)
                                                                                                                                                                                                              cccTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTCCGGCTTCATACCGGTCCCAGCCTTCCTTGTCTGCCAGTACCATTTCCACAACGTCATAGCCGAGATCCTCTCACGCTCCGACAAGTCCTGGCAGCGTAAGGAAATCTCTTGTCCAGCTGACGCCGCCGGCCTGGGCTATATCTTCTGTTGCTGGTAGCTGACCCCAGTAAGGTTCGCCGTTCAGCAGGCCACCGCCCCGttcccg   >  4:193032/1‑270 (MQ=255)
                                                                                                                                                                                                                atGACTTCGCCGCCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGAc                                                                                                                              >  8:165077/2‑149 (MQ=255)
                                                                                                                                                                                                                 tGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACa                                                                                                                             >  7:362330/1‑149 (MQ=255)
                                                                                                                                                                                                                  gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTg                                                                                                                                                                                                                <  4:100942/65‑1 (MQ=255)
                                                                                                                                                                                                                  gACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACaa                                                                                                                            <  8:107958/149‑1 (MQ=255)
                                                                                                                                                                                                                          cggcgAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGc                                                                                                                    >  5:44744/1‑149 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                    
CTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAGATCGTCGAACGCTCCGACAAGTCCTGGCAGCGTCAGGAAATCGCTTGTCGAGCTGACGCCGCAGGCCTGGGCTATTTCTTCTGTTGCTGGTAGCTGACGCCAGTAGGGTTCGCCGATAAGCATGATCCCGCCCGGTTTCAGG  >  NC_000913/4532317‑4532797

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: