Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,225 C→T M44I (ATG→ATA
R132R (CGC→CGT
yaaW ←
mbiA →
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,2250CT92.3% 133.8 / 3.3 53M44I (ATG→ATA
R132R (CGC→CGT
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (18/30);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

AGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTA  >  NC_000913/10981‑11469
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGt                                                                                                                                                                                                                                                >  1:54947/1‑251 (MQ=255)
     ttACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCt                                                                                                                                                                                                                                           >  1:348685/1‑251 (MQ=255)
              aaTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCg                                                                                                                                                                                                                                  <  1:342096/251‑1 (MQ=255)
                 tGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGt                                                                                                                                                                                                                               >  1:381524/1‑251 (MQ=255)
                     ccAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTg                                                                                                                                                                                                                           <  1:447729/251‑1 (MQ=255)
                                  cTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGGTTTGCCTTTTTCATTAt                                                                                                                                                                                                              >  1:587282/1‑251 (MQ=255)
                                         ccAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAg                                                                                                                                                                                                       <  1:423824/250‑1 (MQ=255)
                                          cAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGc                                                                                                                                                                                                      <  1:617259/251‑1 (MQ=255)
                                                  gACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATctggg                                                                                                                                                                                              >  1:618039/1‑247 (MQ=255)
                                                            tGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGCTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGcc                                                                                                                                                                                    >  1:67257/1‑251 (MQ=255)
                                                             gTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGcc                                                                                                                                                                                    >  1:38436/1‑250 (MQ=255)
                                                             gTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACAGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTAGTTGCGTCTCAGTACGCTGGAGAGGAGAGTTTTGCCTTTTTCATTATGGGTGAGAAATCGGGCGAAATTTGcc                                                                                                                                                                                    >  2:406097/1‑250 (MQ=255)
                                                                   tttCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCAAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACATTTCGTTGAGACTCAGTACGCTGGAGTGGCGAGGTTTTCCTTTTTCAATATGGGTGAGGAATCGGGAGAAAttgacacaatttt                                                                                                                                                                             >  2:364139/1‑240 (MQ=255)
                                                                    ttCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGtt                                                                                                                                                                             <  1:29909/250‑1 (MQ=255)
                                                                           tCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACt                                                                                                                                                                      <  1:406097/250‑1 (MQ=255)
                                                                             aaTCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTAc                                                                                                                                                                    >  1:429904/1‑250 (MQ=255)
                                                                                cggtttGAATCATCGAGCAAAATGGTCCGTCACAATTAATCGTGTTCGCGCTGTTTGTTGGCGATACTATCGCCACCAAAATGCTTTAATTCCCCAGCCATCAGCTGCCATTTGCGGCGATGTTGCTCGGGATGCTCTTCCATCGATTTAAACTGTTCGTTGCGTATCAGTACGCTGGAGAGGCGACTTTTGCCTTGTTCATTATGGGTGAGCAAGCGGGCGAAAGTTTCCAACTGTGCCTCACTGCAATg                                                                                                                                                                <  2:677448/248‑1 (MQ=255)
                                                                                           accGAGCAAAAGGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGGTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCATCAATCAGCTGGCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATGTAAACAGTTCTTTGCGTATCAGTACGCTTGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGGTCCGCACGACAGTGCTGAAGAAAAt                                                                                                                                                     <  2:278852/249‑1 (MQ=255)
                                                                                                       ggCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAAt                                                                                                                                         <  1:657977/251‑1 (MQ=255)
                                                                                                              aaCAAGTTACCGTGTTCGCGCTGTTTGTTGGCGATACCATCGCCACCAAAATGCTTTAATTCTCCATCAATCAGCTGCCAGTTGCGGCGATGTTTCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTTTCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAAGTTGCCAACTGTTCCGCACTACAATGCGGAAGAAAATCCAGATCTGAATCATTCAg                                                                                                                                  <  2:323740/250‑1 (MQ=255)
                                                                                                                 aagttACCGTGTCCGCGCAGTTTGTTTGCGATTCTATCGCCACCAAAATGCTTTAATTCTCCAGCAATCAGCTGCCAGGTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGGTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTa                                                                                                                               <  2:611597/248‑1 (MQ=255)
                                                                                                                  aTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTc                                                                                                                                    <  1:87135/245‑1 (MQ=255)
                                                                                                                  aTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTaa                                                                                                                              <  2:73109/251‑1 (MQ=255)
                                                                                                                  aTTTACCGTGTCCACGCAGTTTGTTGGCGCTACTATCGCCACCCACATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTTCCACCTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCaaatcg                                                                                                                              >  2:87135/1‑246 (MQ=255)
                                                                                                                    ttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAAt                                                                                                                             <  1:169097/250‑1 (MQ=255)
                                                                                                                      aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTaa                                                                                                                          <  1:424272/251‑1 (MQ=255)
                                                                                                                      aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTaa                                                                                                                          <  2:638239/251‑1 (MQ=255)
                                                                                                                       ccGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCgg                                                                                                                                                                                                >  1:1031/1‑180 (MQ=255)
                                                                                                                       ccGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCgg                                                                                                                                                                                                <  2:1031/180‑1 (MQ=255)
                                                                                                                                cgcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTAAtttt                                                                                                                >  2:59572/1‑251 (MQ=255)
                                                                                                                                 gcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTTTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCAttttt                                                                                                               <  1:303943/251‑1 (MQ=255)
                                                                                                                                    gTTTGTTGTCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCAGCGATTTAAACAGTTCGTTGCGTATCAGTACGCTTGAGAGTCGAGTTTTGCCttgtt                                                                                                                                                                                                                    <  2:645553/147‑4 (MQ=255)
                                                                                                                                    gTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCAAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCttttt                                                                                                                                                                                                                    >  1:645553/1‑147 (MQ=255)
                                                                                                                                     tttGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAAAATTCATTTTTtgtg                                                                                                           >  1:187195/1‑251 (MQ=255)
                                                                                                                                          tgGCGATAATATCGCCGCCAAAAGGCCGGAATTCGCCATCAATCAGCTGCCAGTTGCGGCGAGGTTGCTCTGGATGCCCCTCCATCGATTTAAACAGTTCGTTGCGTAGCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGAAATCGGGTGAAATTTGCCAATTGTTCCTCACTACAATGCTGTAGAAAATCCAGATCTGAATCAGTCAGGTAAGTAACCTTCATTTTTTGTggttt                                                                                                       <  2:298110/250‑4 (MQ=255)
                                                                                                                                           ggCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAAGCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTc                                                                                                      <  1:91316/250‑1 (MQ=255)
                                                                                                                                                 aCTACCGCCACACAAATGATGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTGCGTTGCGTAGCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTc                                                                                               <  1:297784/251‑1 (MQ=255)
                                                                                                                                                    aTCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTg                                                                                             <  1:522013/250‑1 (MQ=255)
                                                                                                                                                            caAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGt                                                                                                                                                                                                          <  2:345391/133‑1 (MQ=255)
                                                                                                                                                            caAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGt                                                                                                                                                                                                          >  1:345391/1‑133 (MQ=255)
                                                                                                                                                             aaaaatcTGTCATTCTCCAGCAATCAGCTGCCAGGTGCAGCTAGGTTGCTCGGGATGCCCTTCCATCTATTTAAACAGGTCGTGGCGTAGCCGTACGCTGGAGAGGCGTGTTTTGCCTTAGCCATGATGGGTGAGCAAACGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTTAAGAAAATCCAGAGCTGAATCATTCAGGTAAGTAACATTCATTCTTTGTGGTTGCTATAGTCTGGCGTTAgtca                                                                                   <  2:66848/245‑2 (MQ=255)
                                                                                                                                                                  gctggaatgcgccagCAATCAGCTGCCAGTGTCGGAGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGGTCGGGGCGTATCAGTACGCTGGGGAGGCGAGTTTTTCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTACTCACTACAATGCTGAAGAAAAGCCAGATCTGAATCATTCAGGTAATTAACATGCATTTATTGTGGTTTCTATATTCGGGCGTTAGTCGTCGcc                                                                              <  2:428882/237‑1 (MQ=255)
                                                                                                                                                                             ccagCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGACATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAAATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTGGCCGTTAAttttt                                                                    >  2:195586/4‑249 (MQ=255)
                                                                                                                                                                                        cTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGAACTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCat                                                         >  1:72226/1‑250 (MQ=255)
                                                                                                                                                                                           ccAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGCGGCGAGTTTTGCCTTTTTCATTATGGGTGGGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACACTGCTTAAGAAAATCCAGATCTGAATCATTTAGGTAATTAAca                                                                                                                        >  2:341575/1‑184 (MQ=255)
                                                                                                                                                                                           ccAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTACCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCAGTCAGGTAATTAAca                                                                                                                        <  1:341575/184‑1 (MQ=255)
                                                                                                                                                                                                            gCTCGGGATTCCCTTCCATCGCTTTAAACTGTTCTTTGCGTATCACTCCGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGTGTGAGCAATCGGTCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTATCATTCTTTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTAtt                                    <  2:587282/251‑2 (MQ=255)
                                                                                                                                                                                                                    tGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGaa                             >  1:374904/1‑250 (MQ=255)
                                                                                                                                                                                                                      cccTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAg                          <  1:108775/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTg                       >  1:161398/1‑251 (MQ=255)
                                                                                                                                                                                                                           tcATCGATTTAAACAGTTCGTTGCCTACCAGTTCGCTGTAGAGGCGCGTTGTGCCTTTTTCATTATTGGTGAGCAATCTGTCGAAATTTGCCATCTGTTCCTCCCTACAATGCTGCAGAAAATCCAGATCCGACTCATTCAGTTAATTACCATTCATTTTTTGTGTTTTCTATATTCTGGCGGTAGCCGTCGCCGATAATTTTCAGCGTGGCCATATCCTCTGAGTTCTCCGTATGATCGGAAAAGGTGa                      <  2:72226/249‑1 (MQ=255)
                                                                                                                                                                                                                                    tAAACAGTGCGTTGCGCATCATTACGCTGGAGAGGCGAGCTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTTAAGAAAATCCAGATCTTAATCATTCAGGTAATTAACATTCAGTTTTTGTGGCTGCTATATTCTGGCGTTAGTCATCGCCGATAATTTTCAGCGTGTCCAGATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTgaga             <  2:67257/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gTTGCGTATCAGTACGCGGGAGAGGCGAGTTTTGCCTTGTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTa  <  1:61306/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
AGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTA  >  NC_000913/10981‑11469

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