Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,974 C→T L516L (CTG→CTA)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,9740CT97.0% 207.6 / ‑1.6 67L516L (CTG→CTA) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (33/31);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-02

AGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACAT  >  NC_000913/72733‑73213
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
cgtggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGCTCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcagc                                                                                                                                                                                                                                             >  2:366371/5‑248 (MQ=255)
    gtgtgATCAACGCCGTCGATAGTCAGCGCACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                                                                                                                                                                                       >  2:220108/1‑250 (MQ=255)
    gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACCGAGTAGCAGCAGAAtt                                                                                                                                                                                                                                       >  2:140136/1‑250 (MQ=255)
    gtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCCGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAtt                                                                                                                                                                                                                                       >  1:653924/1‑250 (MQ=255)
     tgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATAGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTa                                                                                                                                                                                                                                      >  1:469040/1‑250 (MQ=255)
      gtgATCAACGACGTCGATAGTCAGCGAACCGATGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCGCGGTAGTGTTTCAATCACGGTAAACAGTAGAAAACAGAGTAGCAGCAGAATTAAc                                                                                                                                                                                                                                    <  1:237695/251‑1 (MQ=255)
           caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                                                                                                                                                                               <  1:513750/249‑1 (MQ=255)
           caaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTGACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGt                                                                                                                                                                                                                               <  2:83218/249‑1 (MQ=255)
            aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                                                                                                                                                                                              >  1:702875/2‑251 (MQ=255)
             acgccggcGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGa                                                                                                                                                                                                                             <  1:563089/244‑1 (MQ=255)
                      aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGt                                                                                                                                                                                                                    >  1:30549/1‑251 (MQ=255)
                         gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCt                                                                                                                                                                                                                 <  1:486554/251‑1 (MQ=255)
                             gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt                                                                                                                                                                                                             <  2:653924/251‑1 (MQ=255)
                                   cGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGAt                                                                                                                                                                                                        >  1:67130/1‑250 (MQ=255)
                                      gggctggcGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGTCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCTGTCCTGGCTGCGATAGGa                                                                                                                                                                                                    <  1:156338/250‑1 (MQ=255)
                                       ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                                                                                                                                                                                                   <  1:73713/251‑1 (MQ=255)
                                               tCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAtt                                                                                                                                                                                             <  1:309548/249‑1 (MQ=255)
                                                  gAAAACCGGCTATCAAGTTCAGCAGGGTACTTG‑ACCCGGCGCCGCTTGGCCCGAGGATCGCTACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                                                                                                        <  2:26406/251‑1 (MQ=255)
                                                      aCCGGCGAGCAAGTTACGCAGGGTACTTTTCCCTG‑CGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTag                                                                                                                                                                                    <  2:93357/251‑1 (MQ=255)
                                                            gATCAAGTTCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCAACTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACAGCGCCGTACTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATa                                                                                                                                                                               >  1:219478/1‑250 (MQ=255)
                                                                 agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                                                         <  2:365887/217‑1 (MQ=255)
                                                                 agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                                                         >  1:365887/3‑219 (MQ=255)
                                                                  gttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACACGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGCGATAAAACGGc                                                                                                                                                                        >  1:232299/2‑251 (MQ=255)
                                                                   ttCAGCGGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                                                                        >  2:600197/1‑250 (MQ=255)
                                                                    tCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                                                                        >  1:628535/1‑249 (MQ=255)
                                                                       gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                                                                   <  1:425865/251‑1 (MQ=255)
                                                                       gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                                                                   <  1:274904/251‑1 (MQ=255)
                                                                          gggTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                <  2:30549/251‑1 (MQ=255)
                                                                                ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                                                                      <  1:344441/223‑1 (MQ=255)
                                                                                ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                                                                      >  2:344441/1‑223 (MQ=255)
                                                                                ttttACCTG‑CGCCGCTTGGCACGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGAGCGGAAatc                                                                                                                                                          >  1:564100/1‑251 (MQ=255)
                                                                                          cgccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTAcc                                                                                                                                                 >  1:142364/1‑251 (MQ=255)
                                                                                            ccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCg                                                                                                                                                >  1:429390/1‑250 (MQ=255)
                                                                                             cgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTACGCTAAAACGCATCGGAAAATGTTGGAAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGAGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATCAAACGGCAGGGTGCGGAAATCATCGTTACCGa                                                                                                                                               >  1:609750/1‑250 (MQ=255)
                                                                                                                        cgccgcGTTCTACCGGTAAGCTAAAACGCAACGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTGTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTGGCTGGTAGAGATAAAACGGCAGGGTGCGGAAAGCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑Ca                                                                                                                   <  2:322304/249‑1 (MQ=255)
                                                                                                                         gccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTCGCAGCAGAATTAACGCGGTGACCGCGACGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGAGGAAATCATCGTTACCGAAC‑AACGCCACCCCACCAAAATCAC‑AAATc                                                                                                                >  1:453709/1‑251 (MQ=255)
                                                                                                                             cgTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTAGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAAGCAC‑CAATCGAc                                                                                                             >  1:413406/1‑250 (MQ=255)
                                                                                                                              gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCACCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGTAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCTGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                                                           <  1:42736/251‑1 (MQ=255)
                                                                                                                              gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                                                           <  2:447914/251‑1 (MQ=255)
                                                                                                                                   aCCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCa                                                                                                         <  1:621366/248‑1 (MQ=255)
                                                                                                                                       ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCCCCACGCCAAAATCAC‑CAATCCACAGCACGCaggc                                                                                                  >  2:295334/1‑248 (MQ=255)
                                                                                                                                       ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCaggc                                                                                                  >  1:301157/1‑248 (MQ=255)
                                                                                                                                            cTAAAACGCATCGGCCAATGTTGGTAAATCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAACAGAATTAACGCGGTGACCGCGCCGTCCTGGCCGCGAAAGGAGCCAATTTGCTGGTAGAGCTAAAACGGAAGGGTGCGGCAATCATCGTTACCTACCCAAC‑CCACCCCGCCAAAAT‑ACACAAACGACAACACGCA‑GCaaaaa                                                                                             >  2:356494/1‑248 (MQ=255)
                                                                                                                                                    cATCGGCAAATGTTGGTACAGCCAGGTGATATCAGTCCGTTTTAACCTCTCGACCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCCAGCAAAAACCAGCGcc                                                                                     >  2:324342/1‑251 (MQ=255)
                                                                                                                                                       cGGCAAATGTTTGTAAATCCAGGTGAAATCAGGCAGTTGTAACATCGCGCCCCGGTAGTTTTTCAAGCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACAACGCCAAAATCCC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTgc                                                                                  <  2:653675/251‑1 (MQ=255)
                                                                                                                                                       cGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTg                                                                                   <  1:602166/250‑1 (MQ=255)
                                                                                                                                                                     aaaGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGttt                                                                    <  1:220108/251‑1 (MQ=255)
                                                                                                                                                                        gCCAGGTGATATCAGTCAGTTGTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAg                                                                 <  1:65405/251‑1 (MQ=255)
                                                                                                                                                                          cAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAACCAGCGCCTGCGCCAGTGGACGTTTCCggg                                                               >  2:577812/1‑251 (MQ=255)
                                                                                                                                                                                       tcagtTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTcgac                                                  >  1:560012/1‑248 (MQ=255)
                                                                                                                                                                                              tAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGCCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGCAATCATCGTTACCGAAC‑CACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTACAGGGCGCGCAGCTCGACCACtttt                                           >  2:251620/1‑251 (MQ=255)
                                                                                                                                                                                                  atctCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCCCCCGTGGACGTTTCCGGGCGCGCAACTCGACCACTTTTacgc                                       >  2:334310/4‑248 (MQ=255)
                                                                                                                                                                                                      cGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCgcg                                    <  1:366371/250‑3 (MQ=255)
                                                                                                                                                                                                          ccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                <  2:177041/124‑1 (MQ=255)
                                                                                                                                                                                                          ccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                >  1:177041/1‑124 (MQ=255)
                                                                                                                                                                                                          ccGGTAGTTTTGAAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCTGTGACCTCGCCTTCCTGGCTTCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGTGCGCGCAGCTCGACCACTTTTAAGCGCGACCAg                               <  2:699840/251‑1 (MQ=255)
                                                                                                                                                                                                           cGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGcc                                                                               <  2:611658/202‑1 (MQ=255)
                                                                                                                                                                                                           cGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGAGACCGCGCCGTCCTGGCTGCGATAGGAGCCACTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGcc                                                                               >  1:611658/1‑202 (MQ=255)
                                                                                                                                                                                                               aGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCtt                           <  1:467399/250‑1 (MQ=255)
                                                                                                                                                                                                                    ttCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAAtgc                     >  1:617722/1‑249 (MQ=255)
                                                                                                                                                                                                                      cAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAAtgcc                    <  1:577812/250‑4 (MQ=255)
                                                                                                                                                                                                                          cACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCg               >  1:595239/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGGTAAACAGCAGAAAACAGAGGAGCAGCAGAATGAACGCGTTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCTCGACCAGCCTTCAATGCCCAGCGAc             <  2:564100/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACt            <  1:140136/251‑1 (MQ=255)
                                                                                                                                                                                                                              gTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTg           >  1:578797/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gAAAACAGAGTAGCAGCAGAATTAACGCGATGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAACTCAC‑CAATCGGCAGAACGCAAGCAAAAGCCAGCGCATGCGCCAGCGGAAGATGCAGGACGCGCAGCTCGACCACTTTTAAACGCGACCAGCCTTCAATGCCCAGCGACTAACACAACAt  >  1:564957/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aaaCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACt                                              >  2:187371/1‑204 (MQ=255)
                                                                                                                                                                                                                                          aaaCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACt                                              <  1:187371/204‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
AGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACAT  >  NC_000913/72733‑73213

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