Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,156 G→A A67A (GCC→GCT
P109P (CCG→CCA
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,1560GA95.4% 183.0 / ‑0.3 65A67A (GCC→GCT
P109P (CCG→CCA
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (23/39);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCG  >  NC_000913/10911‑11388
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
tCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCa                                                                                                                                                                                                                                       <  1:60940/249‑1 (MQ=255)
  cAGCTCATTCCCCCTGGCACCGCCCGAGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTAGCGCCACCAAAATGCTGTAATGCTCCAGCAATCa                                                                                                                                                                                                                                   <  2:654241/251‑1 (MQ=255)
       cATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTg                                                                                                                                                                                                                               <  1:217794/250‑1 (MQ=255)
             cccTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACACTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCACCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTg                                                                                                                                                                                                                        >  2:424272/1‑251 (MQ=255)
                    aTCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGCAACATCGAGCAAAATGGCCCGATACACTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTTCGGCGAt                                                                                                                                                                                                                 >  2:547100/1‑251 (MQ=255)
                       gACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGtt                                                                                                                                                                                                              <  1:547100/251‑1 (MQ=255)
                        aCCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATg                                                                                                                                                                                                                <  2:314550/248‑1 (MQ=255)
                           gcgTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGGTGCTc                                                                                                                                                                                                          <  1:108117/251‑1 (MQ=255)
                               gCAGAAACACCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGGTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCTCCAAAATGCTGTAATTCTCCAGCAATCATCTGCCAGTTGCGGAGATGTTGCTCggg                                                                                                                                                                                                       <  2:340001/250‑1 (MQ=255)
                                cAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGAt                                                                                                                                                                                                     <  1:215341/251‑1 (MQ=255)
                                 aGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATg                                                                                                                                                                                                    <  1:9616/251‑1 (MQ=255)
                                 aGAAACCCCTGATTATGTTCCTCGTCCATATTCTTCCAGGGGGTACGCATAAATTGTTCCATTAACTGTTGCTCAATCCCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAATCATCGAGCAAAATGTCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATg                                                                                                                                                                                                    <  2:562704/251‑1 (MQ=255)
                                    aaCTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGccc                                                                                                                                                                                                 >  1:638239/1‑251 (MQ=255)
                                     aCTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCt                                                                                                                                                                                                <  1:679355/251‑1 (MQ=255)
                                       tCCTGATTATGTTCCTCGTCCATATTCTTACAGGTGCTACGCAGACATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCAACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGTTGCCCTTc                                                                                                                                                                                              >  2:388952/1‑251 (MQ=255)
                                                    ccTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaaa                                                                                                                                                                                  <  1:594745/250‑1 (MQ=255)
                                                          ccATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCg                                                                                                                                                                           <  1:662177/251‑1 (MQ=255)
                                                           catatacttcCAGGCGTTACGCAGAAATTGTTCCAGTAACTGTTGATCAATCTCCAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATCCAATTTACCGTGTCCGCGAAGTTTGTTGGCGATACTATCGCCACAAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGTTGCCCTTCACTCGCTTTAAACATTTCGt                                                                                                                                                                          >  2:433074/7‑251 (MQ=255)
                                                                      aGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGt                                                                                                                                                               >  1:54947/1‑251 (MQ=255)
                                                                           ttACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCt                                                                                                                                                          >  1:348685/1‑251 (MQ=255)
                                                                                   aaaTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                                                                                                                                            >  2:287997/1‑225 (MQ=255)
                                                                                   aaaTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                                                                                                                                            <  1:287997/225‑1 (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)
                                                                                                         tCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCgg                                                                                                                                                                                                        >  1:169430/1‑175 (MQ=255)
                                                                                                           aatctCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAAt                                                                                                                                                                                                                                     >  1:537458/5‑144 (MQ=255)
                                                                                                            atctCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTgc                                                                                                                                                                       >  1:657322/4‑205 (MQ=255)
                                                                                                            aactCAAACGTATACATCTCGTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCAAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAATCAGTTCGTTgc                                                                                                                                                                       <  2:657322/202‑1 (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                                         <  2:638239/251‑1 (MQ=255)
                                                                                                                                                                                            aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTaa                                         <  1:424272/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                                                                                                                         >  1:345391/1‑133 (MQ=255)
                                                                                                                                                                                                                                  caAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGt                                                                                                                         <  2:345391/133‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaaatcTGTCATTCTCCAGCAATCAGCTGCCAGGTGCAGCTAGGTTGCTCGGGATGCCCTTCCATCTATTTAAACAGGTCGTGGCGTAGCCGTACGCTGGAGAGGCGTGTTTTGCCTTAGCCATGATGGGTGAGCAAACGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTTAAGAAAATCCAGAGCTGAATCATTCAGGTAAGTAACATTCATTCTTTGTGGTTGCTATAGTCTGGCGTTAgtca  <  2:66848/245‑2 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
TCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCG  >  NC_000913/10911‑11388

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