Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,445 C→A L359L (CTG→CTT)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,4450CA93.8% 172.2 / NA 66L359L (CTG→CTT) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (54/7);  total (59/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAG  >  NC_000913/73200‑73691
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                                                                                                                                                     >  1:628589/1‑251 (MQ=255)
 aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGAAGTTcg                                                                                                                                                                                                                                                    >  2:322371/1‑251 (MQ=255)
            atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTcc                                                                                                                                                                                                                                          >  1:130781/4‑250 (MQ=255)
             tacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATAAACCTGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCCCGCCCGAAGTTCGCGACTGCTccct                                                                                                                                                                                                                                        >  2:71379/3‑249 (MQ=255)
                 gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGGTCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCACTGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGACCAA‑TCCCCGGCATGGCGAGGTTCAACATGCCGCTCATCTCAAGCGCCTGACCCGCCAGCATTTTCAGCCGCGCCCGAAGTTCGCGACGGCTACAtcac                                                                                                                                                                                                                                     >  2:27616/1‑247 (MQ=255)
                      gggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGACCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATGGCGGCAGGCCAATGGTG‑TTGTTGAGCATAAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                                                                                                                                                               <  1:251620/251‑1 (MQ=255)
                             gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAAGGATAGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAAGGGTG‑TTGGTGAGCAGTAAAAAGAA‑GCCGGTCGCC‑AGGACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCCGACCAGCCAGCATTTTCTACAGAGCCAGGAGTACGCGACGGCTCCATAGCAGAATCATGGTc                                                                                                                                                                                                                         >  2:368314/1‑250 (MQ=255)
                                atCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                                                                                                                                                      >  1:108139/1‑250 (MQ=255)
                                   gcgcATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCCCAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAACGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGCGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCGCGCCCGCACTTCGCGCCCGCTCCCTAGCCGCCTCCTGGTCAGccac                                                                                                                                                                                                                   >  2:77243/1‑247 (MQ=255)
                                    cgcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTAc                                                                                                                                                                                                                >  1:168520/1‑251 (MQ=255)
                                         cGGGGTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGAAAAAAGAA‑GCCGGGAGCC‑AGCACAA‑TCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CCACGcac                                                                                                                                                                                                            >  1:124368/1‑250 (MQ=255)
                                            gTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCCGCA‑CTACGCACAAt                                                                                                                                                                                                         >  1:492711/1‑250 (MQ=255)
                                              tttCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCCGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGATCCCTAGCAGCATCATGGTCAGCA‑CTACGCCCAATac                                                                                                                                                                                                       >  1:710645/1‑250 (MQ=255)
                                              ttcccagttcttacggcgcATAAGGGATCGC‑ATATAACGTATGGGTAAAAATCACAGTGCCTGCAGCAGATGGCGGCAGCCCAATGGTG‑TTGTTGAGCATAAAAAAGAG‑GCCGGTAGCC‑AG‑ACAACTCCCCGGCATGGCGAGTATCAACATGCCTCCCCTCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATac                                                                                                                                                                                                       <  2:560012/235‑1 (MQ=255)
                                                  cagtaCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCAGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTATGCCGCGCCCGAAGTTCGGGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCtt                                                                                                                                                                                                  >  2:498160/5‑250 (MQ=255)
                                                     taCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTgccc                                                                                                                                                                                               >  1:170675/2‑250 (MQ=255)
                                                           cAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGCG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAgcc                                                                                                                                                                                         >  2:110044/1‑250 (MQ=255)
                                                              cgaaTAAGGGATCCCCAT‑TACTCCATTGGTAAAAATCACAATGCCGCCAGCAGCTTGCGGC‑GGGCAATGGGGGTGGTTGAACAG‑AAAAAGAAAGCGGGTGCCC‑ACCACAA‑CCCCCGGCATGGCGAGACCCAACATGCCGCTCATCTCCAGCGCTGGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CGACGCACAATACACCTGCCGCCAGCGCaa                                                                                                                                                                                      <  2:31981/248‑1 (MQ=255)
                                                               gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTTGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCCCAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                                                                                                                                                     >  1:315422/1‑251 (MQ=255)
                                                               gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                                                                                                                                                     >  1:150998/1‑251 (MQ=255)
                                                                 aTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAc                                                                                                                                                                                   >  1:534879/1‑251 (MQ=255)
                                                                       gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                                                            >  1:172415/1‑251 (MQ=255)
                                                                       gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                                                            >  1:361226/1‑251 (MQ=255)
                                                                       gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                                                            >  1:378357/1‑251 (MQ=255)
                                                                                   aCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAgc                                                                                                                                                                 >  1:91573/1‑251 (MQ=255)
                                                                                     gCATTGGTGAAAATAACTATGCCGTCAGCAGATTGCGGCAGGCCCATGGTG‑TTGTTGCGCAGTAAACAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCACTCATCTCCACCGCCTGACCCGCCAGCAATTTATGCCGCGCCCGCAGCTCGCGCATGCTCCCTAACCGCATCATGGTCACCAACAACAC‑CAACACCCCTGACGCCAGCGCCATCC‑CCAACGAGGCCcacacccc                                                                                                                                                               >  2:618329/1‑247 (MQ=255)
                                                                                       aTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCt                                                                                                                                                             >  1:78648/1‑251 (MQ=255)
                                                                                          ggTGAAAATCACAATGCCGTCAGCAGATTGCGGGAGGCCAATGGTG‑TTGTAGAGCAGAAAAAAGAA‑GCGGGTAGCC‑AGCACAA‑CCCCCGGCAGGGCGGGGATCAACAGGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGcc                                                                                                                                                          <  1:289090/251‑1 (MQ=255)
                                                                                           gagAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCCTGGGCAGCA‑CTACGCACAATACACCTCACGCCCGCGCCACAC‑GCCACGCGGTCCACAgcccccgcc                                                                                                                                                          >  2:28280/3‑243 (MQ=255)
                                                                                             gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCACCGAGGGCCACAGCGCCTGCcacc                                                                                                                                                       >  2:18106/1‑250 (MQ=255)
                                                                                             gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTg                                                                                                                                                                                                  >  1:292238/1‑209 (MQ=255)
                                                                                              aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGGCCAAAGCGCCTGCcac                                                                                                                                                        >  1:39739/1‑249 (MQ=255)
                                                                                                                   gATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCACTAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCAACGGATGTGCCCGCACTTc                                                                                                                                 >  2:555823/1‑251 (MQ=255)
                                                                                                                         gggcaggccagtggtgttttggAGCAGTAAAAAGAA‑GCCGG‑AGCCCAGCACAA‑CCCCCGGCATGGGGAGGAAAAACAGGCCCCTCTTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCa                                                                                                                            <  1:555823/230‑1 (MQ=255)
                                                                                                                            aGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATCCACATGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAAcgg                                                                                                                        >  2:156659/1‑249 (MQ=255)
                                                                                                                               ccaatggtg‑tTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTACCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCCTGGTCAGCA‑CTACGCACCATACACCTGCCGCCAGCGCCATAC‑GCAACGAGGTCCACAGCGCATTCCACACCCCCGGTTGAGCCAGCACTTCCGACAACCGGCg                                                                                                                     >  2:32313/7‑251 (MQ=255)
                                                                                                                                                  aGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGACGATTTAGCCCATCGACGAt                                                                                                   >  1:90608/1‑251 (MQ=255)
                                                                                                                                                  aGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGAATAGCGAGGAAAAACATGCCGCTCATCTCCAGCGCCTGACAAGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCCGCGCAATAC‑GCAACGAGGTCAACAGCGCATGCCACAGCACCGGCTGTGCCAGCCATTCCGGCAACTGGCGATTTATCCCATCGACGAt                                                                                                   >  1:308805/1‑251 (MQ=255)
                                                                                                                                                    tAAAAAGCA‑GCCGAAAGCC‑AGCACAA‑TCCACGTCAGGGCGAGGATAAAAATGCCGCTCATCTCCAACGCCTGACCCGCCAGCATTTTAAGCCGCGCCCGACGATCGCGACTGCTCCTTAGCAGCATCCTGGTCAGCA‑CTACGCACAATACACCTGCCGCCCGCGCAA‑ACAGAAACGAGGTCCACACCGCCGGCAACAGCAACGGTTGAGCCAGCACCTGCGGCAACTGGCCACTTAGCCCATccacgtacc                                                                                                 >  2:639821/1‑247 (MQ=255)
                                                                                                                                                     aaaaaGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCAc                                                                                                >  1:253242/1‑251 (MQ=255)
                                                                                                                                                        aagaa‑gCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGc                                                                                             >  1:429579/1‑251 (MQ=255)
                                                                                                                                                                   tAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATAAACATGCCGCACATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCTCCCGACGCTCGCGACTGCTCCCTAGCAGAATCATGGTCAGCA‑CTACGCACAATCCACCTGCCGCCCGCGCAATCC‑GCAACGAGCTCCACAGCGCCTGCCACCGCACCGGTTGTGCCAGCACTTCCTGCAAAGGGCGATTTAGCCCACCGACGATCCCCCCCAGCAACGGt                                                                                   >  2:433258/1‑251 (MQ=255)
                                                                                                                                                                       c‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGgc                                                                                >  1:442278/1‑250 (MQ=255)
                                                                                                                                                                       c‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATAGACGATCACCGCCAGCAACGGTGgc                                                                                >  2:300624/1‑250 (MQ=255)
                                                                                                                                                                           cacaA‑TCACCGGCCAGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACACGCCAGCATTTTCTGCCGCGCCAGAAGTTCGCGACTGCACCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACAACTGCCGCCAGCGCAACAC‑GCAACGAGGACCAAAGCGCCGGCCAACGCACCGGTTGTGCCAGCACTTCCGGCCACCGGCGCTTTAGCCCATCAACAATCACCGCCAgacacaggtgaaga                                                                             >  2:340289/1‑237 (MQ=255)
                                                                                                                                                                            acaA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGACAGCAACGGTGgcagca                                                                            >  2:72938/1‑250 (MQ=255)
                                                                                                                                                                               a‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAAc                                                                          >  1:680889/1‑249 (MQ=255)
                                                                                                                                                                                       gCATGGAGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCAGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGAACTTCCGGCAACGGGCGATTTAGCCCATCGAAGATCACCGCCAGCAACGGAGGCAGCGACagcagcagc                                                                 >  2:525529/1‑251 (MQ=255)
                                                                                                                                                                                                 gATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca                                                        >  1:534243/1‑250 (MQ=255)
                                                                                                                                                                                                         tGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGAGGCCGAAATAGCCGCAGCGCCAGCACAATTAACAcc                                               >  2:485717/1‑251 (MQ=255)
                                                                                                                                                                                                          gCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCg                                              >  1:220389/1‑251 (MQ=255)
                                                                                                                                                                                                            cGCTCATCTCCAGCGCCTGACCCGACAGCATTTACTGCCGCGCCCGACGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCCACGGTGGAAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtg                                            >  1:242181/1‑251 (MQ=255)
                                                                                                                                                                                                             gCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCCACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCACCAGCAGCAGCGCCAGCACAATTaaacccgtgc                                           >  2:591415/1‑243 (MQ=255)
                                                                                                                                                                                                              cTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCAc                                                                                                                                    >  1:287573/1‑161 (MQ=255)
                                                                                                                                                                                                               tCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCCTGGTCAGCA‑CTACGCACACTACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGTTTGTGCCAGCACTTCCGGCAACTGGCGATTTAACCCATCGACGATCACCGCCAGCACCGGTGGCAGCAACCGCAGCAGCGCAAGcgaaactagacccgtgtcc                                         >  2:496099/1‑233 (MQ=255)
                                                                                                                                                                                                                 aTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCACCGCCGATCACCGCCAGCAACGGTGGCAGCCACAGCAGCAGCGCCAGCACACTTAACACCGAGTcccc                                       >  2:647953/1‑248 (MQ=255)
                                                                                                                                                                                                                  tctcCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGCTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCGTCGACGCTCACCGCCAGCAACGGAGGCAGCACCAGCAGCAGCGCCCGCCCAATTAACACCGTGTCGCa                                       >  2:97435/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccccagcgcgtcaCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCAAAGCAGCATCATGGTCAGCA‑CTACGCAAAATACCCCTGCCGCCAGCGCAATAC‑GCAACGAGTTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTAAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCaaa                                     <  1:479102/239‑1 (MQ=255)
                                                                                                                                                                                                                      cAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc                                  >  1:505423/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATCGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCCGCACTTCCGGCAACTGGCGATTTAGCCCACCGCCGCTCACCGCCAGCAACGGTCGCAGCAACAGCAGCAGCGCCCGCACAATTAACACCGGGACGCACATGCggat                              >  2:230898/1‑248 (MQ=255)
                                                                                                                                                                                                                           ccTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACCCCGTTGCGCAAATGCGGCt                              >  2:356275/1‑250 (MQ=255)
                                                                                                                                                                                                                              gACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGc                          >  1:431524/1‑251 (MQ=255)
                                                                                                                                                                                                                                cccGCCAGCATTTTCTGCGCCGCCCGAAGTTCGCGACTGCCCCATAGCAGCATCAGGGTCAGCA‑CTCCGCACAATCCCCCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAg                        <  1:18106/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTACAACCGTGTCGCAAATGCGGCTATGCCAACGATc                  >  2:197318/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGt          >  1:537503/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cccGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACCCCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGCCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCCGCAACCTCAGCTGCGCCCGCACACTTACCACCGTCTCGCAACTGCGGCGATTCAGCAGATCGTccagttccccccac  >  2:80427/1‑239 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
GACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAG  >  NC_000913/73200‑73691

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