Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,331 A→G T397T (ACT→ACC thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,3310AG95.5% 203.7 / 0.5 66T397T (ACT→ACCthiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (52/11);  total (54/12)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

CGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTG  >  NC_000913/73087‑73573
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtgttg                                                                                                                                                                                                                                                 >  1:628762/1‑251 (MQ=255)
cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtgttg                                                                                                                                                                                                                                                 >  1:547762/1‑251 (MQ=255)
 gCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtgttgt                                                                                                                                                                                                                                                >  1:239667/1‑251 (MQ=255)
   cAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGttgttg                                                                                                                                                                                                                                              >  1:257383/1‑251 (MQ=255)
        tCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGc                                                                                                                                                                                                                                           >  1:489217/1‑249 (MQ=255)
        ccACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAg                                                                                                                                                                                                                                         <  2:702875/250‑1 (MQ=255)
           ccAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTaa                                                                                                                                                                                                                                      <  2:595239/251‑1 (MQ=255)
            cAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTaaa                                                                                                                                                                                                                                     <  1:574164/251‑1 (MQ=255)
                 gACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCAGCAGGCCACTGGGGTTGTTGAGCAGTAAAaaga                                                                                                                                                                                                                                 >  2:358489/1‑250 (MQ=255)
                    aGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGcc                                                                                                                                                                                                                            >  1:238228/1‑251 (MQ=255)
                           aggCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGcc                                                                                                                                                                                                                     >  1:497060/3‑251 (MQ=255)
                                     cAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATcc                                                                                                                                                                                                          >  1:582967/1‑251 (MQ=255)
                                               gcCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATTGCGGCAGGCCAATGGGGGTTATGAGCAGTAAAAAGAAGCAGGTACCCAGCACAATCCCCGGAATGGc                                                                                                                                                                                                >  2:479102/1‑251 (MQ=255)
                                                 cAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCg                                                                                                                                                                                               >  1:235300/1‑250 (MQ=255)
                                                  aGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCCTGGCg                                                                                                                                                                                               >  1:444623/1‑249 (MQ=255)
                                                     ggACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGAACCATCCCCGGCATGGCGAggtt                                                                                                                                                                                          >  2:409091/1‑249 (MQ=255)
                                                            tCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATg                                                                                                                                                                                   >  1:648279/1‑251 (MQ=255)
                                                                           cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACCTACTGTAGCGGGCGGTGATCTCGCGCCTCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGAATTTGTGAAAATCACAATGACGTCATCAGATTGCGGCAGGCCAATGGTGTTGTTGCGAAGTAAAAAGCAGCCGGTAGCCAGCACAATCCCCGGCATAGCTAGGATCAACGTGCCgatcctttcaacc                                                                                                                                                                    >  2:124544/3‑239 (MQ=255)
                                                                           cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGc                                                                                                                                                                    <  1:532623/249‑1 (MQ=255)
                                                                            gaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTCCTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACCATGCCGTCAGCAGATTGCGGCAGGCCGATGGTGTTGATGAGCAGGAAAAAGAAGCCGGTAGGCAGCACAATCCCCGGCATGGCGAGGATACACGTGCCGCTCATCTccggcg                                                                                                                                                                   >  2:30381/2‑247 (MQ=255)
                                                                                        gcgcgACCAGCCTTCAATGCCCAGCGACTGACACAACATACGGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCa                                                                                                                                                       <  1:409091/247‑1 (MQ=255)
                                                                                         cgcgACCAGCCTTTAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATAGCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAg                                                                                                                                                      <  1:487243/248‑1 (MQ=255)
                                                                                                aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                               >  1:701057/1‑251 (MQ=255)
                                                                                                aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                               >  1:15299/1‑251 (MQ=255)
                                                                                                    ttCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACAAGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                           >  2:675116/1‑251 (MQ=255)
                                                                                                       aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcc                                                                                                                                        >  1:180463/5‑250 (MQ=255)
                                                                                                            ccAGCGACTGACACAACATACTGTAGCGGGCGGTGATCTCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGCAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAACAGCCGGCAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCACATCTCAAGCGCCTGACCCGCCCGCATTTTCTGcccaccccc                                                                                                                                    >  2:655913/1‑243 (MQ=255)
                                                                                                                 gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAg                                                                                                                                                                                                                                            >  2:522224/1‑143 (MQ=255)
                                                                                                                 gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAg                                                                                                                                                                                                                                            <  1:522224/143‑1 (MQ=255)
                                                                                                                 gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                              >  1:628589/1‑251 (MQ=255)
                                                                                                                  aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGAAGTTcg                                                                                                                             >  2:322371/1‑251 (MQ=255)
                                                                                                                          aaCATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCCGCAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCCCCTGACCCCACAGCCTTTTCTGCCGCGccctaaattctcgactgcc                                                                                                                     >  2:142582/1‑235 (MQ=255)
                                                                                                                             atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTcc                                                                                                                   >  1:130781/4‑250 (MQ=255)
                                                                                                                              tacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATAAACCTGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCCCGCCCGAAGTTCGCGACTGCTccct                                                                                                                 >  2:71379/3‑249 (MQ=255)
                                                                                                                                  gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGGTCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCACTGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGACCAATCCCCGGCATGGCGAGGTTCAACATGCCGCTCATCTCAAGCGCCTGACCCGCCAGCATTTTCAGCCGCGCCCGAAGTTCGCGACGGCTACAtcac                                                                                                              >  2:27616/1‑247 (MQ=255)
                                                                                                                                       gggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGACCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATGGCGGCAGGCCAATGGTGTTGTTGAGCATAAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                                        <  1:251620/251‑1 (MQ=255)
                                                                                                                                              gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAAGGATAGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAAGGGTGTTGGTGAGCAGTAAAAAGAAGCCGGTCGCCAGGACAATCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCCGACCAGCCAGCATTTTCTACAGAGCCAGGAGTACGCGACGGCTCCATAGCAGAATCATGGTc                                                                                                  >  2:368314/1‑250 (MQ=255)
                                                                                                                                                 atCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                               >  1:108139/1‑250 (MQ=255)
                                                                                                                                                    gcgcATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCATTAACGCATTGGTGAAAATCCCAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGGTTGTTGAGCAGTAAAACGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGCGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCGCGCCCGCACTTCGCGCCCGCTCCCTAGCCGCCTCCTGGTCAGccac                                                                                            >  2:77243/1‑247 (MQ=255)
                                                                                                                                                     cgcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAc                                                                                         >  1:168520/1‑251 (MQ=255)
                                                                                                                                                          cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtgtt                                                                                                                                                                                                                                                  <  1:180645/96‑1 (MQ=255)
                                                                                                                                                          cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGAGGCAGGCCAATGGtgtt                                                                                                                                                                                                                                                  >  2:180645/1‑96 (MQ=255)
                                                                                                                                                          cGGGGTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGGTTGTTGAGCAGAAAAAAGAAGCCGGGAGCCAGCACAATCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACCACGcac                                                                                     >  1:124368/1‑250 (MQ=255)
                                                                                                                                                             gTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCCGCACTACGCACAAt                                                                                  >  1:492711/1‑250 (MQ=255)
                                                                                                                                                               tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCCGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGATCCCTAGCAGCATCATGGTCAGCACTACGCCCAATac                                                                                >  1:710645/1‑250 (MQ=255)
                                                                                                                                                               ttcccagttcttacggcgcATAAGGGATCGCATATAACGTATGGGTAAAAATCACAGTGCCTGCAGCAGATGGCGGCAGCCCAATGGTGTTGTTGAGCATAAAAAAGAGGCCGGTAGCCAGACAACTCCCCGGCATGGCGAGTATCAACATGCCTCCCCTCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATac                                                                                <  2:560012/235‑1 (MQ=255)
                                                                                                                                                                   cagtaCTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCAGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTATGCCGCGCCCGAAGTTCGGGACTGCTCCATAGCAGCATCATGGACAGCACTACGCACAATACACCtt                                                                           >  2:498160/5‑250 (MQ=255)
                                                                                                                                                                      taCTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccc                                                                        >  1:170675/2‑250 (MQ=255)
                                                                                                                                                                            cAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGCGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAgcc                                                                  >  2:110044/1‑250 (MQ=255)
                                                                                                                                                                               cgaaTAAGGGATCCCCATTACTCCATTGGTAAAAATCACAATGCCGCCAGCAGCTTGCGGCGGGCAATGGGGGTGGTTGAACAGAAAAAGAAAGCGGGTGCCCACCACAACCCCCGGCATGGCGAGACCCAACATGCCGCTCATCTCCAGCGCTGGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACGACGCACAATACACCTGCCGCCAGCGCaa                                                               <  2:31981/248‑1 (MQ=255)
                                                                                                                                                                                gcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTTGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCCCAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAAt                                                              >  1:315422/1‑251 (MQ=255)
                                                                                                                                                                                gcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAAt                                                              >  1:150998/1‑251 (MQ=255)
                                                                                                                                                                                  aTAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATAc                                                            >  1:534879/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:378357/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:361226/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:172415/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgc                                           >  1:91573/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATTGGTGAAAATAACTATGCCGTCAGCAGATTGCGGCAGGCCCATGGTGTTGTTGCGCAGTAAACAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCACTCATCTCCACCGCCTGACCCGCCAGCAATTTATGCCGCGCCCGCAGCTCGCGCATGCTCCCTAACCGCATCATGGTCACCAACAACACCAACACCCCTGACGCCAGCGCCATCCCCAACGAGGCCcacacccc                                         >  2:618329/1‑247 (MQ=255)
                                                                                                                                                                                                        aTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                                       >  1:78648/1‑251 (MQ=255)
                                                                                                                                                                                                           ggTGAAAATCACAATGCCGTCAGCAGATTGCGGGAGGCCAATGGTGTTGTAGAGCAGAAAAAAGAAGCGGGTAGCCAGCACAACCCCCGGCAGGGCGGGGATCAACAGGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGcc                                    <  1:289090/251‑1 (MQ=255)
                                                                                                                                                                                                            gagAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCCTGGGCAGCACTACGCACAATACACCTCACGCCCGCGCCACACGCCACGCGGTCCACAgcccccgcc                                    >  2:28280/3‑243 (MQ=255)
                                                                                                                                                                                                              gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCACCGAGGGCCACAGCGCCTGCcacc                                 >  2:18106/1‑250 (MQ=255)
                                                                                                                                                                                                              gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTg                                                                           >  1:292238/1‑209 (MQ=255)
                                                                                                                                                                                                               aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGGCCAAAGCGCCTGCcac                                  >  1:39739/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCACTACGCAACGAGGTCCACAGCGCCTGCCACAGCAACGGATGTGCCCGCACTTc           >  2:555823/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATCCACATGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAAcgg  >  2:156659/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
CGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTG  >  NC_000913/73087‑73573

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