Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 15,073 A→T L302L (CTA→CTT dnaJ → chaperone protein DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,0730AT97.8% 128.9 / NA 45L302L (CTA→CTTdnaJchaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (21/23);  minor base G (0/1);  total (21/24)
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

GGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGC  >  NC_000913/14827‑15311
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ggACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTcc                                                                                                                                                                                                                                            >  1:49177/1‑251 (MQ=255)
 gACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAATCTTTTCCg                                                                                                                                                                                                                                           <  1:185297/251‑1 (MQ=255)
              gCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTaa                                                                                                                                                                                                                               >  2:121119/1‑250 (MQ=255)
                aTCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCCGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTGAAgg                                                                                                                                                                                                                            >  2:106488/1‑251 (MQ=255)
                          cGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCCCCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCt                                                                                                                                                                                                                  <  1:236726/251‑1 (MQ=255)
                                 gggcgAAGCGGGCGAGCATGGCGCACCGGCCGGCGATCTGTACGTTCCGGTTAAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGACGTCCCGATCAACTTCGCTATGGCGGCGATGGGTGGCGAAATCGCAATACCGGCCCTTGATGATCGCGACAAACTGAAAGTGCCGTGCGCAACCCAAAACGGTAAGCTTTTCCCTCTGGGAGGTAAAGGCgaccagtgtggccgcg                                                                                                                                                                                                           >  2:1954/5‑236 (MQ=255)
                                    tGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGGGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGACGTCATGTCTGTCCGCGGTg                                                                                                                                                                                                        >  1:10603/1‑251 (MQ=255)
                                              gAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTGCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAggg                                                                                                                                                                                              <  1:233473/251‑1 (MQ=255)
                                               aGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCTCAGGGt                                                                                                                                                                                             <  2:113923/251‑1 (MQ=255)
                                                            ggCAGGCGATCTGTTCGTTCAGGTTCAGGTGAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGcc                                                                                                                                                                                <  1:229179/251‑1 (MQ=255)
                                                                     tCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTgtcg                                                                                                                                                                       >  1:229243/1‑251 (MQ=255)
                                                                               caggttcaggttAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAacac                                                                                                                                                              <  1:3890/250‑1 (MQ=255)
                                                                                 ggttcaggttAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCgg                                                                                                                                                           <  2:245582/251‑1 (MQ=255)
                                                                                   ttcaggttAAACAGCACCCGATTTTCGAGCGTGCAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAAACCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCAGCGTTGTCGTCGAAACAACGGt                                                                                                                                                          >  1:94361/1‑250 (MQ=255)
                                                                                             aCAGCACCCGAGTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATGTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACg                                                                                                                                               <  1:106488/251‑1 (MQ=255)
                                                                                                 cACCCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGaga                                                                                                                                            >  2:212512/1‑248 (MQ=255)
                                                                                                        ttttCGAGCGTGAATGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATAGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGGGCCTGGCGAAACCCAGAACGGTAAGCTTTTCCGTATGCGCGGTACAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGGTTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGGACgggaagcagga                                                                                                                                     >  2:194332/1‑240 (MQ=255)
                                                                                                             gAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAgct                                                                                                                               >  2:188439/1‑251 (MQ=255)
                                                                                                                cGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTCTGCCGGCGCTGGGTGGCGAAATCCAACTACAGACCCTTGCTGGCCGCGTCCAACTGACAGTGACTGGCGGAACCCCGCCCGGTAAGCTTTTCCGTATGCGCGGTACAGGCGTCGGGTCTGTCCGCGGCGGCGCACCGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAgctgct                                                                                                                            >  1:133428/1‑251 (MQ=255)
                                                                                                                 gTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGCAATACCGACCCTTGATGGTCGCGTCAAACTGAAAGCGCCTGGCGAAACCCAGACCGGGAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTCGGCCTGAACGCGAAGCAGAAACAgctggtg                                                                                                                           >  1:163263/1‑248 (MQ=255)
                                                                                                                 gTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGctgctg                                                                                                                           >  1:103773/1‑251 (MQ=255)
                                                                                                                          acaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctg                                                                                                                  >  1:216412/1‑251 (MQ=255)
                                                                                                                            aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgca                                                                                                                <  1:134156/251‑1 (MQ=255)
                                                                                                                            aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgca                                                                                                                <  1:121119/251‑1 (MQ=255)
                                                                                                                            aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgc                                                                                                                 >  1:63590/1‑250 (MQ=255)
                                                                                                                                 ctATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgcaaga                                                                                                             <  2:165362/248‑1 (MQ=255)
                                                                                                                                        gAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCgg                                                                                                    <  1:145904/251‑1 (MQ=255)
                                                                                                                                         aaGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGt                                                                                                   <  1:247557/251‑1 (MQ=255)
                                                                                                                                               cGATCAGCTTCGCTATGGCGGCGCTCGGTGGCGAAATAGAAGAACCGACCCTTGATGGTCGCGCCCACCTGAAAGTGCCTGGCGAAACCCCGACCGGGAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGGATGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGACGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAAACGCTGCTGCAAAAGCTGCCAGAAAGcgttgggggcccc                                                                                             >  2:38220/1‑239 (MQ=255)
                                                                                                                                                          gCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGGGGCCCAACCGGGGAGc                                                                                   >  2:69687/1‑250 (MQ=255)
                                                                                                                                                             aTGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGGATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCTGTGAGcac                                                                                 <  2:125152/249‑1 (MQ=255)
                                                                                                                                                             aTGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGGCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGcaca                                                                                <  2:71884/250‑1 (MQ=255)
                                                                                                                                                                   gcgcTGGGTGGCGAAAACGAAGTACCGACCCTTGATTGTCGCGGCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGGCCGCGGTGGCGCACAGGGTGAGTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGccc                                                                         <  2:214015/251‑1 (MQ=255)
                                                                                                                                                                             gCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCGTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGTTTCGGTGGTCCAACCGGTGAGCACAACAGCCCGCGCTCAAag                                                               <  2:124613/251‑1 (MQ=255)
                                                                                                                                                                                      aaGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAACGCTTCGGTGGCCCAACCGGTGAGCACACCAGCCCGCGCTCAAAGAGCTTCttt                                                      >  1:55629/1‑251 (MQ=255)
                                                                                                                                                                                           ccGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATgg                                                 <  1:188439/251‑1 (MQ=255)
                                                                                                                                                                                                   tGATGGTCGCGTCAAACTGAAAGTGCCTGGCTAAACCCAGACCGGTAAGCTGTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACGCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaa                                           <  2:216412/249‑1 (MQ=255)
                                                                                                                                                                                                    gATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAATCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaagaa                                        <  1:202968/251‑1 (MQ=255)
                                                                                                                                                                                                           gcgTCAAACTGAAAGTGGCTGGGGAAACCCATACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCATGTCTGTCCGCGGTGGCGCACAGGGTGAGGTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGGTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCGAACCGGTGAGCACAACAGCCTGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGtttttt                                 <  2:16560/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaCTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTGCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGGCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTgacgac                           <  2:49177/251‑1 (MQ=255)
                                                                                                                                                                                                                     gAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCATCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTg                        >  1:56761/1‑250 (MQ=255)
                                                                                                                                                                                                                        aGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGAccc                    <  1:212512/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGCGATTTGCTGTGCCGCGTTGTCGTCGGAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTcc          >  2:230625/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCaaa     >  1:98942/1‑249 (MQ=255)
                                                                                                                                                                                                                                          gACCGGTAAGCTTTTCCTTATTCGCGGTAGAGGCTTCAAGTCTGTCCGCGGTGGCGCACACGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACGCCGGTAGGTCTGAACGAGAAGCAGAAACAGCTGCTGCAGGTGCTGCGAGAAAGCTTCGGTGGCCTAACCGGTGATCACAACAGCCCGCTCGCAGAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGc  <  1:120085/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
GGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGC  >  NC_000913/14827‑15311

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