Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,796,197 G→T A73A (GCG→GCT waaC → ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,1970GT98.4% 211.7 / NA 62A73A (GCG→GCTwaaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/38);  total (23/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 = 9.98e-01

TGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGC  >  NC_000913/3795956‑3796437
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCaa                                                                                                                                                                                                                                          <  2:104708/250‑1 (MQ=255)
 gTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCaaa                                                                                                                                                                                                                                         <  1:331404/250‑1 (MQ=255)
   tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGc                                                                                                                                                                                                                                       <  2:377580/250‑1 (MQ=255)
        ggggAAGCCTGACTGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCg                                                                                                                                                                                                                                 <  2:345391/249‑1 (MQ=255)
                    cGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACa                                                                                                                                                                                                                     <  1:117701/251‑1 (MQ=255)
                      gATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAg                                                                                                                                                                                                                   <  1:209574/251‑1 (MQ=255)
                          cGGGTTTTGATCGTTAAAACATCGTCGATGCGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGc                                                                                                                                                                                                                  >  1:326800/1‑248 (MQ=255)
                             gTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAgaga                                                                                                                                                                                                             >  1:190962/1‑250 (MQ=255)
                                ttGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACt                                                                                                                                                                                                          >  1:376643/1‑250 (MQ=255)
                                   aTCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGa                                                                                                                                                                                                      <  1:216756/251‑1 (MQ=255)
                                           aaCATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGtt                                                                                                                                                                                               >  2:205151/1‑250 (MQ=255)
                                              aTCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCATTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCg                                                                                                                                                                                           >  1:22249/1‑251 (MQ=255)
                                                cgtcgATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGat                                                                                                                                                                                         <  2:315204/251‑2 (MQ=255)
                                                            aTGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTg                                                                                                                                                                             >  2:363267/1‑251 (MQ=255)
                                                              gTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGt                                                                                                                                                                           >  1:311138/1‑251 (MQ=255)
                                                                   ccATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTaaaaa                                                                                                                                                                      <  1:217794/251‑1 (MQ=255)
                                                                       aCGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAgcg                                                                                                                                                                   <  2:123215/250‑1 (MQ=255)
                                                                         gTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCgcgg                                                                                                                                                                <  2:313254/251‑1 (MQ=255)
                                                                           tGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCgcggc                                                                                                                                                               <  1:201402/250‑1 (MQ=255)
                                                                                 cactcactGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTg                                                                                                                                                        >  1:157048/1‑251 (MQ=255)
                                                                                  actcactGATGCCCCGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGATCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTg                                                                                                                                                        >  1:69861/1‑250 (MQ=255)
                                                                                     cactGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTGTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGAcg                                                                                                                                                     <  1:66701/250‑1 (MQ=255)
                                                                                        tGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTc                                                                                                                                                 >  1:416528/1‑251 (MQ=255)
                                                                                        ggATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGGTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTc                                                                                                                                                 <  2:311138/250‑1 (MQ=255)
                                                                                         gATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCt                                                                                                                                                <  1:348449/251‑1 (MQ=255)
                                                                                          aTGCCCAGCAGGCAAGCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTGGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTg                                                                                                                                               <  1:390603/251‑1 (MQ=255)
                                                                                                      cAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGCTTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCAAGGCGt                                                                                                                                    >  2:388341/1‑250 (MQ=255)
                                                                                                          cccAGGGATTAAGTTTGACTGGTTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATGCGTCGCTTTCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAGCTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgc                                                                                                                               <  1:416469/251‑1 (MQ=255)
                                                                                                                  ttAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGAACAGATTCCTTGCTGGCACGCTGCCGTTGAGAGAGTTATTCCTGTGGCAATAAGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGACGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAGgcatggattg                                                                                                                       >  1:418621/1‑247 (MQ=255)
                                                                                                                   tAAGTTTGATTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattg                                                                                                                       <  2:78550/250‑5 (MQ=255)
                                                                                                                   tAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattg                                                                                                                       <  1:90498/250‑5 (MQ=255)
                                                                                                                    aaGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGCGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGACCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATAGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggattgg                                                                                                                      >  1:427394/1‑245 (MQ=255)
                                                                                                                         tGACTGGGTGGTGTAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAATCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCGCAGGGGCTGTTAAAAAGCGCGGCGCGGGTGACGCGTCTGGCGCAGGGCGTAAAGCATTGATTGGACTgg                                                                                                                 <  2:434070/250‑1 (MQ=255)
                                                                                                                              gggTGGTGGAAGAAGGGTTCGCACAGATTCCGTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAAcc                                                                                                           <  2:190962/251‑1 (MQ=255)
                                                                                                                                gtggtgGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGc                                                                                                         >  1:128292/1‑251 (MQ=255)
                                                                                                                                  ggtggAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTc                                                                                                       <  2:326800/251‑1 (MQ=255)
                                                                                                                                            gggTGCGCACAGAGTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCGGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCTGCCCCCATAAAAGCTGAACGCATAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGTTGACGCGTCGGGCGCATGGCGTAAAGCATGGATTGGACTGGCAATCCGCTCGCGAAccgtt                                                                                             <  2:37955/251‑4 (MQ=255)
                                                                                                                                                  gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGCTGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagccac                                                                                       >  2:330714/1‑250 (MQ=255)
                                                                                                                                                  gCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagcc                                                                                         >  1:66211/1‑249 (MQ=255)
                                                                                                                                                   cacaGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTagccagc                                                                                      >  2:258296/1‑251 (MQ=255)
                                                                                                                                                     caGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAATCCGCTCGCTAACCTTTagccagc                                                                                      <  2:398355/249‑1 (MQ=255)
                                                                                                                                                      aGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTg                                                                                   <  1:207815/251‑1 (MQ=255)
                                                                                                                                                        aTTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGtt                                                                                 >  1:222317/1‑251 (MQ=255)
                                                                                                                                                         ttccttccTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGttt                                                                                >  1:93760/1‑251 (MQ=255)
                                                                                                                                                              tcctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACa                                                                            <  2:66211/250‑1 (MQ=255)
                                                                                                                                                               cctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACaa                                                                           <  1:258296/250‑1 (MQ=255)
                                                                                                                                                                 tGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatgg                                                                                                                           <  1:201291/200‑1 (MQ=255)
                                                                                                                                                                 tGGCACGCTGCCGTTGAGAGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatgg                                                                                                                           >  2:201291/1‑200 (MQ=255)
                                                                                                                                                                  ggCACGCTGCCGTTGATCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTTGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTGCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTTGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGTCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                       <  1:173339/251‑1 (MQ=255)
                                                                                                                                                                  ggCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGt                                                                       >  1:90213/1‑251 (MQ=255)
                                                                                                                                                                               ggAGCGAGTTATTCCTGTGGCAATACGTCGCTGGTGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAATCATCATATTg                                                          <  1:102405/250‑1 (MQ=255)
                                                                                                                                                                                  gtgAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCaa                                                       <  2:112926/249‑1 (MQ=255)
                                                                                                                                                                                      gTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca                                                   <  2:222317/251‑1 (MQ=255)
                                                                                                                                                                                                     aTACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGa                                    <  2:128292/251‑1 (MQ=255)
                                                                                                                                                                                                       aCGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAAc                                  <  2:22249/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAcc                             <  1:363267/251‑1 (MQ=255)
                                                                                                                                                                                                               ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTTGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAGTGCAAAACAGCAGCACGCCGTAGAGCGCAccc                            <  2:254669/249‑1 (MQ=255)
                                                                                                                                                                                                               ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCcg                           <  1:205151/250‑1 (MQ=255)
                                                                                                                                                                                                               ggCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACGATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAccc                            <  1:43721/249‑1 (MQ=255)
                                                                                                                                                                                                                  gTAAAGCGTGGTTCTCGGCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCATTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGGCTGTCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCTCTCTCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGaa                       <  1:267040/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCCCCATAAAAGCTGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGtct      >  1:415183/1‑249 (MQ=255)
                                                                                                                                                                                                                                       ccATAAACGCTGAACGCAACGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGc  >  1:318111/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
TGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGC  >  NC_000913/3795956‑3796437

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