Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 27,367 A→C A25A (GCA→GCC rihC → ribonucleoside hydrolase RihC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091327,3670AC98.6% 225.4 / NA 72A25A (GCA→GCCrihCribonucleoside hydrolase RihC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/37);  total (34/38)
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.00e+00

GGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCACCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGC  >  NC_000913/27132‑27609
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
ggCGGTGTGGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGGTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCAGTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTg                                                                                                                                                                                                                                      <  2:671525/251‑1 (MQ=255)
    gtggtgAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAAc                                                                                                                                                                                                                                  <  1:522603/251‑1 (MQ=255)
      ggagAAGCCATTCCGCTGGAAGGCCGTGAAGAAACCATTGTTTTCGCCGTGCCGAAAGAGCTGCGTGTCGATATTAGTGAAGACGTTTAAGTCATTAGCCGCCTACGTTATGCGAAAATGCCGTTCTTGTTACTGGCATTTTTTATGAAGAAAACATGCGTTTCCCTATCTTCATCGATACTGACCTCGGCATTGCAGATGCCGCCGCCATTGCCGACGCGATTTTTGCCACCGATCTCGACATGCAACtg                                                                                                                                                                                                                                >  2:355629/4‑251 (MQ=255)
          aaGCCATTCCTCTTGAAGGTCGTGTAGAAAACAGTGTTTTAGAAGGGCCGAAAGAGCTGCGGGTCGATATTCTGGAAGGCGATTAAGTCATTAGCAGGTTGAGTTATGCGAAAACGCTGGTCTTGTTACCGGCATTTTTTATGGAGAATACATGCGTTTACCGATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCAGTGCCGCCGCTATTTTTGCCCCCGAACGCGAGCTGCAACtgatga                                                                                                                                                                                                                            <  2:648313/251‑1 (MQ=255)
                 tCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCg                                                                                                                                                                                                                      <  1:626639/250‑1 (MQ=255)
                  ccGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTc                                                                                                                                                                                                                    <  1:510819/251‑1 (MQ=255)
                  ccGCTGGAAGGCCGTGAAGAAAACATGGTGTTCGAAGTTCCGAAAGAGTTGCGTGTCGTTATTCGTGAAGTCTATTAAGGCATTAGCAGACTAAGGTATGCGAAAATGCCGGTCTGGTTACCGGCATTTTTTATGGATAAAACAGGCGTTTACCTATCTTCCTCGATACTGACCCCTGCTTTGACGATGCCGTCGCCATTGCCGCCGCGATTTTGGCCCCCGAACCCGACCTGCAACTGATGACCTCCGTc                                                                                                                                                                                                                    <  2:386927/251‑1 (MQ=255)
                     cTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGCTTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTcgcg                                                                                                                                                                                                                 >  2:512504/1‑251 (MQ=255)
                      tGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCgg                                                                                                                                                                                                                >  1:138571/1‑251 (MQ=255)
                               gTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGtct                                                                                                                                                                                                       <  1:180762/251‑1 (MQ=255)
                                ggaagaaAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGtctc                                                                                                                                                                                                      <  1:369203/250‑1 (MQ=255)
                                 gaagaaAACATTGTTTTCGGAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGACGACCACCGTCGCGGGTAATGCCTCg                                                                                                                                                                                                     >  2:208789/1‑251 (MQ=255)
                                           ttGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTATTGTCTCGGTTGCGaaa                                                                                                                                                                                            >  2:325348/1‑250 (MQ=255)
                                                  cGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTAccc                                                                                                                                                                                     <  1:110912/250‑1 (MQ=255)
                                                    aaGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCAGGTAATGTCTCGGTTGAGAAAACTACCCGc                                                                                                                                                                                   >  1:393145/1‑250 (MQ=255)
                                                      gTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTAGCGGGTAATGCCTCGGTTGAGAACACTACCCgaagt                                                                                                                                                                                >  2:108540/1‑247 (MQ=255)
                                                         ccGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGCTTAAGTCATTAGCAGCCAAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGcc                                                                                                                                                                             >  1:340731/1‑251 (MQ=255)
                                                          cGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGccc                                                                                                                                                                            >  1:153034/1‑251 (MQ=255)
                                                            aaaGAGCTGCGTGTCGATATTCGGGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGGCCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTg                                                                                                                                                                          <  2:138571/251‑1 (MQ=255)
                                                                agCTGCGTTGCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAATTTATGCGAAAATGCCGGTCTTTTTACCGGCATTTTTTATGGAGAAAACATGCTTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCTTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAAc                                                                                                                                                                      <  1:429225/251‑1 (MQ=255)
                                                                agCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAAATACCCGCAATGCCCTGCAAc                                                                                                                                                                      >  1:409694/1‑251 (MQ=255)
                                                                      gtgtCGATATTCGTGAAGTCGATTAAGTCATTAGCAGTCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACtgctgc                                                                                                                                                                <  1:653032/251‑1 (MQ=255)
                                                                           gATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTc                                                                                                                                                           >  1:153603/1‑251 (MQ=255)
                                                                            atatTCGAGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCAGACGCGATTTTTGCCCCCGAACTCGACCAGCACCTGGTGACCACCGTCGCGGGTACTGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCt                                                                                                                                                          >  2:4504/1‑251 (MQ=255)
                                                                                     gatgcGATTAAGGCATTAGCAGCCTAATTTATGCTAAAATGCCTGTTGTGTTACCGGCTTTTTTGATGGATAAAACATGCTTTTACCTATCCTCCTCGATTCTGACCCCGTCATTGACGATGCCGTCGCCATTGCCGCCGCGGTTTGTGCCCCCGAACTCGACCGGCAACGGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAAGGCCCTGCAGCTGCTGCATTTCTGGAATGCgg                                                                                                                                                 <  2:409694/247‑1 (MQ=255)
                                                                                     aaGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCgg                                                                                                                                                 <  1:440505/251‑1 (MQ=255)
                                                                                           aTTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTc                                                                                                                                           <  1:42170/251‑1 (MQ=255)
                                                                                           aTTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCCGCATTTCTGGAATGCGGAGATTc                                                                                                                                           >  1:508058/1‑251 (MQ=255)
                                                                                              aaGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCg                                                                                                                                         <  1:16135/250‑1 (MQ=255)
                                                                                                            ctgatctATGCGAAAATGCCGGTCAAGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGGCCCCGACATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCACGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCAGCAATGCCCTGCAACTGCGGCCTTTCTGGAATGAGGAGAATCAGCTCGCCCaacgagcc                                                                                                                          >  2:456914/7‑245 (MQ=255)
                                                                                                            cTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGcc                                                                                                                          >  1:623472/1‑251 (MQ=255)
                                                                                                               aGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCt                                                                                                                       <  1:334679/251‑1 (MQ=255)
                                                                                                                   aTGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCtgtg                                                                                                                   <  1:200873/250‑1 (MQ=255)
                                                                                                                     gCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCa                                                                                                                <  1:621649/251‑1 (MQ=255)
                                                                                                                       gAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCAGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGAGCCACt                                                                                                              >  2:156321/1‑251 (MQ=255)
                                                                                                                          aaTGCCGGCATTGTTACCGGCATTTTTTATGGAGCAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGACGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGAAACTACTGCATTTCTGGAATGCGGAGATTCCGCCCGCCCCAGGGGCCGCTGTGcaaatggt                                                                                                           >  2:640792/1‑244 (MQ=255)
                                                                                                                                 gTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACgcgc                                                                                                     >  1:344892/1‑250 (MQ=255)
                                                                                                                                  tCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCAc                                                                                                   >  1:478/1‑251 (MQ=255)
                                                                                                                                         aCCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTg                                                                                             >  1:84730/1‑250 (MQ=255)
                                                                                                                                          ccGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTCCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGAt                                                                                           <  1:126214/251‑1 (MQ=255)
                                                                                                                                                   tttATGGACAAAACATGCGTTTACCTATCTTCCTAGATACTGACCACGGCATTAACGATGCAGTCGCCATTGCCGCCGCGCTTTTTGCCCCCGCACTCGTCCTGTAACTGATGACCACCGTCGCGGGTAATGTCTCGGCTGGAAAAACTAACCGCAATGCCCTGCCACTCATGCTTTTCTGGAACGCGGAGATTCCGTGCGCCCAAGCGCCCGCTGAGCCCCTGGTAcgacacccgcgtggtgcgccatct                                                                                  >  2:228574/1‑229 (MQ=255)
                                                                                                                                                             aaaCATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCg                                                                        >  1:544219/1‑251 (MQ=255)
                                                                                                                                                                 aTGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATCGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGc                                                                         <  1:47282/246‑1 (MQ=255)
                                                                                                                                                                 aTGCGTTTACCTATCTTACTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACGACCGTCGCGGGTAATGTCTCGGTTGAGCAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGGATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCCCTGGTACGCGCACCGCGTGAGGCGGCATCtgggaccggccgat                                                                     >  2:47282/1‑238 (MQ=255)
                                                                                                                                                                  ggtgtttcccattctcccccGATACGGCCCCCGGCATGGCCGAGGCCTTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCCCCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCGGCATTTCTGGAATGCTGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCg                                                                   <  2:83010/232‑1 (MQ=255)
                                                                                                                                                                      tttACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGa                                                                >  1:565770/1‑250 (MQ=255)
                                                                                                                                                                          ccTATCTTCCCCGATTCTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACGCGACCTGAAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGATTCGTGAAtggc                                                           <  2:429558/251‑1 (MQ=255)
                                                                                                                                                                                  ccTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCTGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACg                                                   >  1:236340/1‑251 (MQ=255)
                                                                                                                                                                                  ccTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACg                                                   >  1:167334/1‑251 (MQ=255)
                                                                                                                                                                                     cGATACTGCCCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCTGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGTAGATTCCGCTCGCCCAAGTTGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGa                                                  <  1:421532/249‑1 (MQ=255)
                                                                                                                                                                                         aCTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGga                                                                                                                                                <  1:111374/152‑1 (MQ=255)
                                                                                                                                                                                         aCTGACCACGGCATTGACTATGCCGCCGTCAGTGCGGCCGCGTTGTTTGCACCCGAACTCGCCCTGCAACTGATGACCTCCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCCGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCTACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTttg                                             <  2:508058/250‑1 (MQ=255)
                                                                                                                                                                                         aCTGACACCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGga                                                                                                                                                >  2:111374/1‑152 (MQ=255)
                                                                                                                                                                                           tGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACCttggtg                                          >  2:535567/1‑248 (MQ=255)
                                                                                                                                                                                              ccccGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGc                                       >  1:320111/1‑251 (MQ=255)
                                                                                                                                                                                                    cATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGGAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACGGGTACGCGCACCGCGTGATGCGGCGTCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAAcc                                 <  1:640792/251‑1 (MQ=255)
                                                                                                                                                                                                          cGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGa                                                                >  2:7885/1‑214 (MQ=255)
                                                                                                                                                                                                          cGATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGa                                                                <  1:7885/214‑1 (MQ=255)
                                                                                                                                                                                                           gATGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGcc                          <  1:19496/251‑1 (MQ=255)
                                                                                                                                                                                                            aTGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCg                         >  1:47496/1‑251 (MQ=255)
                                                                                                                                                                                                            aTGCCGTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCg                         <  1:85147/251‑1 (MQ=255)
                                                                                                                                                                                                              gCCGTCGCCATTGCCTCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACGACCCGCAATGCCCTGCAACTGCTGCATTTCTTGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACTGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACATCCGAAAGCCGCt                       <  1:504842/251‑1 (MQ=255)
                                                                                                                                                                                                                cGTCGCCATTGCCGCCGCGATTTTCGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCg                     <  1:236943/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTCGCCATTGCCGCCGCGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCgg                    >  1:509849/1‑251 (MQ=255)
                                                                                                                                                                                                                    cctgttGCGGCCGCGATTTTCGCCCCCGAACGCTACATGTAACTTATGACCGCCGTCGAGGGTGATGTCTCGGTTGAGAAAACTACCCGCAATGCCCGGCAAAGGCTGCATGTCTGGAATGCGGAGATGCCGCTCGCCCAAGGGGGCGCTGGTGCGCTGGTACGCGCACCGCTTGATGCGGCATCTGGGCACGGCTAATCGGGAATGTCTGGGTACGACTTTGTTGATCTCTACCTAAAGCCGCtttggat                 <  2:617231/247‑7 (MQ=255)
                                                                                                                                                                                                                       aTTGCCGCCGCGATTTTTGCCCCCGAACGCGACCTGCAACTGATGACCACCGTCGCGTGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTGCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATTCGGCATCTGTGCACGGCGATTCGGGAATGGCTGGTTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGAtccc              <  1:672222/251‑4 (MQ=255)
                                                                                                                                                                                                                             gccgcGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGtt        >  1:234133/1‑251 (MQ=255)
                                                                                                                                                                                                                              ccgcGATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGttt       >  1:108418/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cgATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTg    <  2:478/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTgg   <  1:133497/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTTTTGCCCCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGc  <  1:341547/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTTTTGCCCCCGAACTCGACCGGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTgg   <  2:565770/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
GGCGGTGGTGAAGCCATTCCGCTGGAAGGCCGTGAAGAAAACATTGTTTTCGAAGTGCCGAAAGAGCTGCGTGTCGATATTCGTGAAGTCGATTAAGTCATTAGCAGCCTAAGTTATGCGAAAATGCCGGTCTTGTTACCGGCATTTTTTATGGAGAAAACATGCGTTTACCTATCTTCCTCGATACTGACCCCGGCATTGACGATGCCGTCGCCATTGCCGCCGCGATTTTTGCACCCGAACTCGACCTGCAACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGC  >  NC_000913/27132‑27609

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