Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 628,713 G→A D55N (GAC→AAC)  entA → 2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913628,7130GA95.8% 141.9 / NA 48D55N (GAC→AAC) entA2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/27);  minor base C (2/0);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCA  >  NC_000913/628469‑628954
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
aTGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGCAAGGTACAGGGTTTGATCAAGCGCTCACTCAGGAGCAATATCCCTTTTCGACCGAAGTGATGGCTGTTGCCAACGCTg                                                                                                                                                                                                                                              >  1:135787/1‑251 (MQ=255)
     gATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGGTTGTTGAGGCTGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCa                                                                                                                                                                                                                                          <  2:205704/250‑1 (MQ=255)
          cGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGGGAAAGTTACCAGATTTGATTCAGCGTTTACTCAAGAGCAATATCCCTTTGCGACCGAAATGATGGTTGTTGGCCACGCTGCGCAgtgcgc                                                                                                                                                                                                                                    >  2:125075/1‑246 (MQ=255)
                      gCTGGCGAAAAACCCGACCATCGACGCCTGGTGGCAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGCAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAAAATTCCTTTGCGACCGAAGTGATGGGTTTTGCCAACGATGCGCAGGTCGCGCAAGTGTGTcc                                                                                                                                                                                                                        >  1:164717/1‑250 (MQ=255)
                      gCTGGCGAAAAACCCGACCATCGACGCCTGGTGGACGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGGCGCGCAAGTGTGTc                                                                                                                                                                                                                         >  1:123332/1‑250 (MQ=255)
                                     gACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGGGGTGAAATAATGGATGTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATGGGCTACGGCACGTCGGTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAAGATCCCTTTGCGACCGGAGTGATGGATGTTGCCAACTCTGCTCAGGTCGCGCAAGTGTGTCAGCTACTGTTAGCTGa                                                                                                                                                                                                         <  2:269089/251‑1 (MQ=255)
                                       ccATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGaa                                                                                                                                                                                                        >  1:25050/1‑250 (MQ=255)
                                            gACGCCTGGTGGAAGCTTCTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGGAAAGGTATCTGCTGCGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTTATCAAGCGTTCACTCAGGAGCTATATCCCTTGGCGACCGAAGTGATGGATGTTGCCAACGTGGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAg                                                                                                                                                                                                  <  1:292341/251‑1 (MQ=255)
                                              cGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGGGCGAAAGGTACAGGGTTTGATTAAGGGGTCCCTCAGGAGCAATATTCCTTTGCGACCGAAGTGATGGATGTTGCCCACGCCGCGCAGGGCGCGCAAGAGTGTCAGCGGCTGGTTGGTGAAACggcggg                                                                                                                                                                                                >  1:298853/1‑247 (MQ=255)
                                                 cTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGGGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGGCCGAAGTGATGGATGTTGCCAACGCTGCGCAGGGCGCGCAAGGGTGTCAGCGACTGTTAGCTGAAACGGAGCGACt                                                                                                                                                                                             >  1:353282/1‑251 (MQ=255)
                                                   ggtggAAGCTACTATCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGGGGAAGTGAAAAGGTTTGAATAAGCGATAACTCAAGAGCAATAACCCCTTTGGAACGAAATGGTGGATGTTGCCAACGCTGAGCAAGACGGGCAAGTGTGGCCGCGGATGTTAGCTGAAAACGAGCGACttg                                                                                                                                                                                           >  2:103791/1‑249 (MQ=255)
                                                      ggAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGCAAGTACAGGGTTTGATCTAGCGTTCACACCGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGGCAAAGCTGCGCAGGGCGCGCAAGAGTGTCAGCGACTGTTAGCTGAAACCGAgacaccggg                                                                                                                                                                                          >  1:405544/1‑241 (MQ=255)
                                                        aaGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGAcgcg                                                                                                                                                                                      <  2:82743/251‑1 (MQ=255)
                                                           cTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGAACGAAGTGATGGATGTTGCCAAAGCCGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGGTGGAACGGAGAGACTGGGAGCGCTg                                                                                                                                                                                   >  2:344301/1‑251 (MQ=255)
                                                             gctcacccGAGAGGTGGAAGAATGGATGGCAGCGGTAAAAATGTCTGGGCAACCGGCGCAGGTAAAAGTATCGGCTACGCCACGGCGCTGGCGCTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTTATCAAGCGTGCACTCAGGAGGAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTg                                                                                                                                                                                   <  2:7417/244‑1 (MQ=255)
                                                                        aGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTATGTCAGCGACTGTTAGCTGAAAAGGGGCGACTGGACGCGCTGGGCAATGcggcgg                                                                                                                                                                      >  1:443457/1‑251 (MQ=255)
                                                                            gAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACTCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCTGCGGGAAt                                                                                                                                                                  <  2:298853/251‑1 (MQ=255)
                                                                              aataatGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCCCGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAAttt                                                                                                                                                                <  2:427573/251‑1 (MQ=255)
                                                                                    ggATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGGACAAGCGTTCACTCAGGAGCAATATCCCTTTTGGGCCGAAGTGGTGGATGTTGCCAACGCTGCGCCGGGTGCGCAAGGGTGTCAGCGACTTTTTGCTGAAACGGAGCGACTTGGCGCGCCGGGCAATGCGGCggggggttttcgac                                                                                                                                                          >  2:310758/1‑240 (MQ=255)
                                                                                        ttCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATgg                                                                                                                                                       <  2:305542/250‑1 (MQ=255)
                                                                                        ttCAGCGGTAAAAATGGCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTCCACTCAGGAGCAATATCCCGTTGCGACCGAAGTGAGGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGTGAATTTTACGCATgg                                                                                                                                                       <  1:17633/250‑1 (MQ=255)
                                                                                          cAGCGGTAAAAATGTCTGGGTAACCGGAGCAGGGAAATGTATCGGCTACTCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCTCAAGTGTGTCAGCGACTGTTAGCTGAAACCGAGCGACTGGACGCGCTGGTCAAGGCGACGTGAATTTTGCGCATTGgc                                                                                                                                                     <  2:437772/250‑1 (MQ=255)
                                                                                                     aTGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAg                                                                                                                                         <  1:358811/251‑1 (MQ=255)
                                                                                                     aTGTCTGGGTAACCGGCGAAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTTCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGGCGCGCAAATGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGGATTTTACGCaatggggggaccggtc                                                                                                                                           >  1:362848/1‑234 (MQ=255)
                                                                                                       gTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCt                                                                                                                                       <  2:362848/251‑1 (MQ=255)
                                                                                                           gggTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCa                                                                                                                                     <  2:273140/249‑1 (MQ=255)
                                                                                                                                 gaTCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGGGTTCACTCATGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACtt                                                                                                             <  2:163656/250‑1 (MQ=255)
                                                                                                                                         aCGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGtt                                                                                                     <  1:121132/251‑1 (MQ=255)
                                                                                                                                                       gCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGc                                                                                       <  1:152825/251‑1 (MQ=255)
                                                                                                                                                            tgttgAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCTGTTAACGTCGGCGGTGCGttt                                                                                   <  2:16578/250‑1 (MQ=255)
                                                                                                                                                             gttgAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCGCGCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGGGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTa                                                                                  <  1:94185/250‑1 (MQ=255)
                                                                                                                                                                              aaGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTATGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAAcc                                                                <  2:443457/251‑1 (MQ=255)
                                                                                                                                                                                     aGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCa                                                               <  2:78148/245‑1 (MQ=255)
                                                                                                                                                                                     aGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTAAGCATGGGCGCGACCGATCAGCCCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGGTAAAGTCGGCGGGGCGGTTAAACTGGTCCAGCAAAcacagaacg                                                         >  1:78148/1‑243 (MQ=255)
                                                                                                                                                                                              tCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCCGTTCCAGCAAACCATGAACCAGTTTCGcc                                                >  1:88838/1‑251 (MQ=255)
                                                                                                                                                                                                aaGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACtt                                                                                                             <  1:132896/189‑1 (MQ=255)
                                                                                                                                                                                                    gTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAgc                                          <  2:76231/251‑1 (MQ=255)
                                                                                                                                                                                                            aGGAGCAATATCCCTTTGCGACCGGAGTGATGGATGTTGCCAACCCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTTAAGAGGACTGGCCGCAGACTTTTGCGGTTAACCTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATCAACCACATTCGTCGTCCgcggggcggg                                  >  1:387522/1‑251 (MQ=255)
                                                                                                                                                                                                               aGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCCTGAACCAGGTTCGCCGTCAACGgggcggggcg                               >  1:356533/1‑251 (MQ=255)
                                                                                                                                                                                                               aGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGgggcggggcg                               <  1:230812/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaTATCCCTTTGCGACCGAAGTGATGGATGTTGCCAAAGCTGCGCAGGTCGAGCAAGTGTGTCAGCGACTGGTAGCTGAAACGGAGCGACTGGAGGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGAtt                            <  1:433503/251‑1 (MQ=255)
                                                                                                                                                                                                                   atatCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTg                                                                            <  1:317116/202‑1 (MQ=255)
                                                                                                                                                                                                                   aganCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCATCTCAGTAAAGAGGACTGGCAGGAGAGTGTTGCGGTTAACGTCGGCGGTGCGTTTAACGTGTTCCAGCAAACCATGAACCAGTTTCTCCGTCAGCGTGGCGGGGCGATTg                           <  2:413361/249‑1 (MQ=255)
                                                                                                                                                                                                                     atCCCTTTGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACATTTGCGGTTAACGTCGGCGGTGGGTTTAAACTGTTCCAGCAAAACATGAAACAGTTTCGCCGTCAGcgggggggggggggtgtt                         >  1:355613/1‑243 (MQ=255)
                                                                                                                                                                                                                     atCCCTTTGCGACCGAAGTGATGGATGGTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGGATGTTACGCACGGGCGCGACCGATCAGCTCAGTAAAGAGGACAGGCAGCAGACTTTTGCGGTGAACGGAGTCGGGGCGGGTTACCTAGTTCAGCAAACCCTGAACCAATTTCGCCGTCAGCAgggcggggcggtatgg                         >  2:232125/1‑245 (MQ=255)
                                                                                                                                                                                                                           tttGCGACCGAAGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGGTAACGTCGGCGGGGCGTTTAACCTGTTTCAGCAAAACATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTTACtct                    >  1:328765/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aGTGATGGATGTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGAcg         <  1:145674/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gatgTTGCCAACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGTTAACGTCGGCGGTGCGTTTAACCTGTTTCAGCAAACCATGAACCAGGTTCGCCGTCCGCGGGGCGGGGGGATTGGCACTGTGGCGTCCGACGCCccgcc  >  1:216234/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ATGGTGATATCGACTTTGTCATGCTGGCGAAAAACCCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGTTCCAGCAAACCATGAACCAGTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCA  >  NC_000913/628469‑628954

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