Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,328 G→A I398I (ATC→ATT)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,3280GA93.8% 191.3 / 0.4 66I398I (ATC→ATT) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (50/11);  total (54/12)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTG  >  NC_000913/73087‑73573
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttg                                                                                                                                                                                                                                                 >  1:628762/1‑251 (MQ=255)
cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttg                                                                                                                                                                                                                                                 >  1:547762/1‑251 (MQ=255)
 gCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttgt                                                                                                                                                                                                                                                >  1:239667/1‑251 (MQ=255)
   cAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑ttgttg                                                                                                                                                                                                                                              >  1:257383/1‑251 (MQ=255)
        tCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGc                                                                                                                                                                                                                                           >  1:489217/1‑249 (MQ=255)
        ccACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAg                                                                                                                                                                                                                                         <  2:702875/250‑1 (MQ=255)
           ccAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTaa                                                                                                                                                                                                                                      <  2:595239/251‑1 (MQ=255)
            cAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTaaa                                                                                                                                                                                                                                     <  1:574164/251‑1 (MQ=255)
                 gACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCAGCAGGCCACTGGGG‑TTGTTGAGCAGTAAAaaga                                                                                                                                                                                                                                 >  2:358489/1‑250 (MQ=255)
                    aGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑Gcc                                                                                                                                                                                                                            >  1:238228/1‑251 (MQ=255)
                           aggCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGcc                                                                                                                                                                                                                     >  1:497060/3‑251 (MQ=255)
                                     cAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑Tcc                                                                                                                                                                                                          >  1:582967/1‑251 (MQ=255)
                                               gcCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATTGCGGCAGGCCAATGGGG‑GTTATGAGCAGTAAAAAGAA‑GCAGGTACCCAGCACAA‑TCCCCGGAATGGc                                                                                                                                                                                                >  2:479102/1‑251 (MQ=255)
                                                 cAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCg                                                                                                                                                                                               >  1:235300/1‑250 (MQ=255)
                                                  aGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCCTGGCg                                                                                                                                                                                               >  1:444623/1‑249 (MQ=255)
                                                     ggACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGAACCA‑TCCCCGGCATGGCGAggtt                                                                                                                                                                                          >  2:409091/1‑249 (MQ=255)
                                                            tCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATg                                                                                                                                                                                   >  1:648279/1‑251 (MQ=255)
                                                                           cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACCTACTGTAGCGGGCGGTGATCTCGCGCCTCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGAATTTGTGAAAATCACAATGACGTCATCAGATTGCGGCAGGCCAATGGTG‑TTGTTGCGAAGTAAAAAGCA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATAGCTAGGATCAACGTGCCgatcctttcaacc                                                                                                                                                                    >  2:124544/3‑239 (MQ=255)
                                                                           cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGc                                                                                                                                                                    <  1:532623/249‑1 (MQ=255)
                                                                            gaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTCCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACCATGCCGTCAGCAGATTGCGGCAGGCCGATGGTG‑TTGATGAGCAGGAAAAAGAA‑GCCGGTAGGCAGCACAA‑TCCCCGGCATGGCGAGGATACACGTGCCGCTCATCTccggcg                                                                                                                                                                   >  2:30381/2‑247 (MQ=255)
                                                                                        gcgcgACCAGCCTTCAATGCCCAGCGACTGACACAACATACGGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCa                                                                                                                                                       <  1:409091/247‑1 (MQ=255)
                                                                                         cgcgACCAGCCTTTAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATAGCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAg                                                                                                                                                      <  1:487243/248‑1 (MQ=255)
                                                                                                aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                               >  1:701057/1‑251 (MQ=255)
                                                                                                aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                               >  1:15299/1‑251 (MQ=255)
                                                                                                    ttCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACAAGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                           >  2:675116/1‑251 (MQ=255)
                                                                                                       aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcc                                                                                                                                        >  1:180463/5‑250 (MQ=255)
                                                                                                            ccAGCGACTGACACAACATACTGTAGCGGGCGGTGATCTCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGCAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAACA‑GCCGGCAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCACATCTCAAGCGCCTGACCCGCCCGCATTTTCTGcccaccccc                                                                                                                                    >  2:655913/1‑243 (MQ=255)
                                                                                                                 gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAg                                                                                                                                                                                                                                            >  2:522224/1‑143 (MQ=255)
                                                                                                                 gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAg                                                                                                                                                                                                                                            <  1:522224/143‑1 (MQ=255)
                                                                                                                 gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                              >  1:628589/1‑251 (MQ=255)
                                                                                                                  aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGAAGTTcg                                                                                                                             >  2:322371/1‑251 (MQ=255)
                                                                                                                          aaCATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCCGCAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCCCCTGACCCCACAGCCTTTTCTGCCGCGccctaaattctcgactgcc                                                                                                                     >  2:142582/1‑235 (MQ=255)
                                                                                                                             atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTcc                                                                                                                   >  1:130781/4‑250 (MQ=255)
                                                                                                                              tacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATAAACCTGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCCCGCCCGAAGTTCGCGACTGCTccct                                                                                                                 >  2:71379/3‑249 (MQ=255)
                                                                                                                                  gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGGTCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCACTGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGACCAA‑TCCCCGGCATGGCGAGGTTCAACATGCCGCTCATCTCAAGCGCCTGACCCGCCAGCATTTTCAGCCGCGCCCGAAGTTCGCGACGGCTACAtcac                                                                                                              >  2:27616/1‑247 (MQ=255)
                                                                                                                                       gggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGACCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATGGCGGCAGGCCAATGGTG‑TTGTTGAGCATAAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                                        <  1:251620/251‑1 (MQ=255)
                                                                                                                                              gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAAGGATAGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAAGGGTG‑TTGGTGAGCAGTAAAAAGAA‑GCCGGTCGCCAGGACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCCGACCAGCCAGCATTTTCTACAGAGCCAGGAGTACGCGACGGCTCCATAGCAGAATCATGGTc                                                                                                  >  2:368314/1‑250 (MQ=255)
                                                                                                                                                 atCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                               >  1:108139/1‑250 (MQ=255)
                                                                                                                                                    gcgcATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCCCAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAACGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGCGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCGCGCCCGCACTTCGCGCCCGCTCCCTAGCCGCCTCCTGGTCAGccac                                                                                            >  2:77243/1‑247 (MQ=255)
                                                                                                                                                     cgcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTAc                                                                                         >  1:168520/1‑251 (MQ=255)
                                                                                                                                                          cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑tt                                                                                                                                                                                                                                                  <  1:180645/96‑1 (MQ=255)
                                                                                                                                                          cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGAGGCAGGCCAATGGtg‑tt                                                                                                                                                                                                                                                  >  2:180645/1‑96 (MQ=255)
                                                                                                                                                          cGGGGTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGAAAAAAGAA‑GCCGGGAGCCAGCACAA‑TCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CCACGcac                                                                                     >  1:124368/1‑250 (MQ=255)
                                                                                                                                                             gTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCCGCA‑CTACGCACAAt                                                                                  >  1:492711/1‑250 (MQ=255)
                                                                                                                                                               tttCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCCGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGATCCCTAGCAGCATCATGGTCAGCA‑CTACGCCCAATac                                                                                >  1:710645/1‑250 (MQ=255)
                                                                                                                                                               ttcccagttcttacggcgcATAAGGGATCGC‑ATATAACGTATGGGTAAAAATCACAGTGCCTGCAGCAGATGGCGGCAGCCCAATGGTG‑TTGTTGAGCATAAAAAAGAG‑GCCGGTAGCCAG‑ACAACTCCCCGGCATGGCGAGTATCAACATGCCTCCCCTCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATac                                                                                <  2:560012/235‑1 (MQ=255)
                                                                                                                                                                   cagtaCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCAGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTATGCCGCGCCCGAAGTTCGGGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCtt                                                                           >  2:498160/5‑250 (MQ=255)
                                                                                                                                                                      taCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTgccc                                                                        >  1:170675/2‑250 (MQ=255)
                                                                                                                                                                            cAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGCG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAgcc                                                                  >  2:110044/1‑250 (MQ=255)
                                                                                                                                                                               cgaaTAAGGGATCCCCAT‑TACTCCATTGGTAAAAATCACAATGCCGCCAGCAGCTTGCGGC‑GGGCAATGGGGGTGGTTGAACAG‑AAAAAGAAAGCGGGTGCCCACCACAA‑CCCCCGGCATGGCGAGACCCAACATGCCGCTCATCTCCAGCGCTGGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CGACGCACAATACACCTGCCGCCAGCGCaa                                                               <  2:31981/248‑1 (MQ=255)
                                                                                                                                                                                gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTTGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCCCAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                              >  1:315422/1‑251 (MQ=255)
                                                                                                                                                                                gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                              >  1:150998/1‑251 (MQ=255)
                                                                                                                                                                                  aTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAc                                                            >  1:534879/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:378357/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:361226/1‑251 (MQ=255)
                                                                                                                                                                                        gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                                      >  1:172415/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgc                                           >  1:91573/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATTGGTGAAAATAACTATGCCGTCAGCAGATTGCGGCAGGCCCATGGTG‑TTGTTGCGCAGTAAACAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCACTCATCTCCACCGCCTGACCCGCCAGCAATTTATGCCGCGCCCGCAGCTCGCGCATGCTCCCTAACCGCATCATGGTCACCAACAACAC‑CAACACCCCTGACGCCAGCGCCATCCCCAACGAGGCCcacacccc                                         >  2:618329/1‑247 (MQ=255)
                                                                                                                                                                                                        aTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                                       >  1:78648/1‑251 (MQ=255)
                                                                                                                                                                                                           ggTGAAAATCACAATGCCGTCAGCAGATTGCGGGAGGCCAATGGTG‑TTGTAGAGCAGAAAAAAGAA‑GCGGGTAGCCAGCACAA‑CCCCCGGCAGGGCGGGGATCAACAGGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGcc                                    <  1:289090/251‑1 (MQ=255)
                                                                                                                                                                                                            gagAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCCTGGGCAGCA‑CTACGCACAATACACCTCACGCCCGCGCCACACGCCACGCGGTCCACAgcccccgcc                                    >  2:28280/3‑243 (MQ=255)
                                                                                                                                                                                                              gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCACCGAGGGCCACAGCGCCTGCcacc                                 >  2:18106/1‑250 (MQ=255)
                                                                                                                                                                                                              gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTg                                                                           >  1:292238/1‑209 (MQ=255)
                                                                                                                                                                                                               aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGGCCAAAGCGCCTGCcac                                  >  1:39739/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCACTACGCAACGAGGTCCACAGCGCCTGCCACAGCAACGGATGTGCCCGCACTTc           >  2:555823/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATCCACATGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAAcgg  >  2:156659/1‑249 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
CGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTG  >  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: ‑