Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,338 T→A 14.2% intergenic (+140/‑74) pgm → / → ybfP phosphoglucomutase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,3380TA14.2% 164.3 / 3.6 68intergenic (+140/‑74)pgm/ybfPphosphoglucomutase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (42/15);  new base A (2/9);  total (44/24)
Fisher's exact test for biased strand distribution p-value = 8.65e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  NC_000913/715128‑715562
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
gttaagaaCATCGGAAGCAGATTGAAAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACCCATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                                                                                <  1:96546/226‑1 (MQ=255)
 tGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAAGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGACTGATGCGTTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGGATTTAAGCGATTTTGTAAGCCGGATAAGGCGTTAACGCCGCATc                                                                                                                                                                                                               >  1:448155/1‑229 (MQ=255)
      aaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAACGATTTTGTAGCCCGGATAAGGCGTCCACGCAGCATCCGGcc                                                                                                                                                                                                          >  1:432867/1‑228 (MQ=255)
         aTCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTGAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaaa                                                                                                                                                                                                       <  1:571026/229‑1 (MQ=255)
           cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAaac                                                                                                                                                                                                     >  1:651338/1‑229 (MQ=255)
             cAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCGGATGGTGCAATATATTGACTTTAAGCGTTTTTGTAGGCCGGATAAGGCGTTCACGCGCCATCCGGCCAAaacaa                                                                                                                                                                                                   >  1:473433/1‑229 (MQ=255)
                    aTTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGGTAAGGCGTTCACGCCGAATCCGGCAAAAAAAACGAaaac                                                                                                                                                                                            >  2:64285/1‑229 (MQ=255)
                              aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGCTTTTGTAGGCCGGTAAAGGCGTTCAGGCCGCATCCGGCAAAAAAAACGAAAAATTTTAcaac                                                                                                                                                                                   >  2:286820/1‑228 (MQ=38)
                               agCTGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACAATTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  1:375314/229‑1 (MQ=255)
                                agggTTTAGATTTGTTAGCGAAGTTTGTAAAACCGGGTAAACTCATTTAAAAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGAAAAACACAGAATTGCCTGATGGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  2:52135/227‑1 (MQ=38)
                                      gagaTTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGAAAAGGCGTTCACGCCGCATCCGGCAAAAGCAACGAACACTTTGTCAACAAAccaaga                                                                                                                                                                          >  2:445466/1‑224 (MQ=255)
                                              ttAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAACTCaa                                                                                                                                                                   >  1:256011/1‑229 (MQ=255)
                                                 ccGAATTTCTCAAACACGCTTAAAACATTTTAAAAAAAAAAGGGCCGGCGTAAGATCCACATTTTTTACTTAGGACAAACACAGAAATGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGTAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGa                                                                                                                                                                <  2:568844/228‑1 (MQ=21)
                                                        tCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCCGGCCTACGAGGATGGTGCAATATATTGAATTTAACCGATTTTGTAGTCCGGCTAAGGCGTTCACGCCCCATCCTGCAAAAACAACAAACACTTTGTCACCAACCTGAgaggccacgggaaatccaa                                                                                                                                                         >  2:611225/1‑211 (MQ=255)
                                                           gAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACTCTTAGCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATg                                                                                                                                                      <  1:624072/229‑1 (MQ=255)
                                                               caaggaggtgtttaagagaacaaaaaaaGGGGGGGCGTAAGATCGCCCTATAGTACGTAGGACAAATGTAGAATTGCCTGATGCGCTACGCTTAGCAGGCCGACGAGGATGGTGCAATATATACAATTTAAGCGATTATGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAAGCCCAATGAAg                                                                                                                                                   <  2:8900/208‑1 (MQ=255)
                                                                 cgcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCAAAAACACCGAACACTTTGTACAAAAACTGAGTAGATCAACGAAAACCAAATaaaaaac                                                                                                                                                >  2:214086/1‑225 (MQ=37)
                                                                 cgcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGa                                                                                                                                                  >  1:346693/1‑227 (MQ=255)
                                                                 cgcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAATAGCCCAAGGAAAACCCAATGAAAAca                                                                                                                                                >  2:280876/1‑228 (MQ=255)
                                                                 cgcgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTAAGTAGCCCAAGGAAATCCCAAAGAAAAcc                                                                                                                                                >  2:422247/1‑229 (MQ=255)
                                                                   cgTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCaa                                                                                                                                              >  1:499886/1‑229 (MQ=255)
                                                                     ttaaCACATTGAAAAAAAAAAGGGCTGTTCCAAGATCGCTCTTTTTTCCGGAGTCCAAACATATAATTGCTTGATTGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACa                                                                                                                                            <  1:744314/227‑1 (MQ=37)
                                                                               ttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTCCGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGCGGAGCTCAAGGAAAACCAAAGAAAGACCAACAagagaagtg                                                                                                                                   >  2:756026/1‑224 (MQ=37)
                                                                                     aaaaaaGGGCGGTCGCAAGGTCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCa                                                                                                                             <  1:411911/229‑1 (MQ=255)
                                                                                       aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCACCAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGACTGGCCGTCCTa                                                                                                                           >  1:334582/1‑229 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       >  1:625598/1‑211 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGGATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       <  2:625597/211‑1 (MQ=255)
                                                                                              cgttCGCAAGAGCGCCTTTTTTTACGTATGACAAACACAGATTTGGCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGAAAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCAc                                                                                                                    <  2:334582/226‑1 (MQ=255)
                                                                                               ggTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGCCAAAAGGTCACt                                                                                                                   >  2:340060/1‑229 (MQ=255)
                                                                                                gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAg                                                                                                                                                   >  1:8900/1‑196 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCAACAGGAGTCTGGGCGTCATAGTGTCACTaa                                                                                                                 >  1:311813/1‑229 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACCAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATATGGTCACTac                                                                                                                 >  1:338225/1‑228 (MQ=255)
                                                                                                         tccgctttttttACATTTGACAAACAAAGAATTGCTGGATGCGCTGCGCTTATGAGGCATACGAGGAGGGGGAAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAAAAAGCAACGAACACTTTGCCAACAAACTGAGTAGCCCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGc                                                                                                          <  1:445466/222‑1 (MQ=37)
                                                                                                             ccTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTTGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACCGAGTAGCCCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGGCATAgtggccccaataaaagccaaca                                                                                                     >  2:563831/1‑210 (MQ=25)
                                                                                                                 ttttACGTATGACAAACACAGAATTGCATGATGCGCTACGCTTATCAGGCCTACGAGGTTGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGCTAAGGCGTTCAGGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGAGTCGCTCAAGGAAATCCCAATGACAACCACCAGGAGTATGGTCGTCATAGTGCAAGTACTTACCGcaccatgat                                                                                                  >  1:315891/1‑220 (MQ=25)
                                                                                                                           gACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCACCAGGAGTCTGGTCGTCATAGTGTAACAAATTACCGCCACATTACTGCTGACttc                                                                                        >  1:205877/1‑226 (MQ=255)
                                                                                                                                            cTGATGGGCTACGTTTACCGGGCCTAAGGGGATGGTGCAATTAATGGAATTAAAGCTATTTTGAAGGCCGGAAAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      <  1:64285/229‑1 (MQ=255)
                                                                                                                                               aTGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGACGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      >  1:543375/1‑226 (MQ=255)
                                                                                                                                                tGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGTCATACTGTCACTAATTACCGCCCCATTACTGCTGACTGCATGCGCCCAACCTGAACAatc                                                                   >  1:116748/1‑228 (MQ=255)
                                                                                                                                                   ccTACGTTTAACAGGCCTACGGGGATGGTGCAATATATGAAATTTAAGCGATTTTGTAGGCCGTAAAAGGCTTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAAtcatca                                                               <  1:422247/228‑1 (MQ=255)
                                                                                                                                                    attacctgacgaggcactcGAGGGTGGTGCAATATAGTAAATTTAAGCAATTTTTGAGGCCGAATAAGGCGTTCCAGCCGGATCCGGCAAAAACAACGAACACTTTGCCAACAAACTGAGTAGCTAAATGAAATCCCACTGAAGCCCAACAGGATGCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATcatcat                                                              <  1:286820/211‑1 (MQ=255)
                                                                                                                                                            aTCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGCCTGCATGCGCCCAACCTGAAAAATCACCATTAGCTg                                                        >  2:481984/1‑227 (MQ=255)
                                                                                                                                                                ggCCTACGGGGGGGGCGAAATATATTAATTTAAAGCGATTTGGTAGCCCGGAAAAGGTGTACACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  <  1:340060/229‑1 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  >  1:133682/1‑228 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGACTGGTCGTCATAGTGTCACTAATTACCGCCACATAACTTCTGACTGCATCCGCCCAACCTGAACAATCATCATAAACTGGCGAcg                                                  >  2:655613/1‑227 (MQ=255)
                                                                                                                                                                    tACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGtt                                              >  1:594494/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTAc                                            >  1:459613/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTACACAACCATCATTAGCTGGCGACTGGTTAc                                            >  1:407763/1‑229 (MQ=255)
                                                                                                                                                                       ggagagggtaaaaTATATTAAATTAAATAAATTTTGTAGGCCGAAAAGGGCCTCCACCACGCATCGGGAAAAAAAAACGAACACTTGAGCAACAAACTGAGTAGCTCAAGGAAACCCCAAGGAACACCAACAGGAGTCTGGTCGTCATAGTGTATCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACt                                           <  2:499886/218‑1 (MQ=255)
                                                                                                                                                                              gTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGACAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTACCTGGCGACTGGTTACTGAcccc                                     >  1:629757/1‑228 (MQ=255)
                                                                                                                                                                               tacAAGATATTGAATTTAAGCAATTTGGAGGCGCGAAAAAGGTTTGCACGCAGCAGCGGGAAAAAACAACAAACCTTTTGCAAACAAACGGAGGAGCTAAAGGAAAGCCCAATGAAACCCAACGGGAGTCGGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:214086/227‑1 (MQ=255)
                                                                                                                                                                                gaaaataatTGAATTTAGGCGATTGTAAAGCCCGAAAAAGGCGTTCCCGCCGCATCCGGCAAAAACAACGACCACTTTGTAAACAAGCTGAGTAGCTCAAGGAAATCCCAATGAAGCCCACCAGGGGTCTGGTCGTCATGTTTTCCCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:535208/221‑1 (MQ=255)
                                                                                                                                                                                   atatatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAgc                               >  1:371009/1‑228 (MQ=255)
                                                                                                                                                                                     atatTTAATTTAAGCGATTTTGTAGGCCGAATAAGGCGTTAACTCCGCAACCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTGGGCATCGTGTCACTAATTACCGCCACATTCCAGCTGGCTCCAGGCGCCCAACCTGAACTAGCATCAGTAGCTAGCGCCCGGTAACTGACCCCAAAAGaaa                             >  2:690943/1‑227 (MQ=255)
                                                                                                                                                                                       atataaatataacggttttGTAGGCGCGATAAGGCATTTCCCCCCACCACCAGAAAAAACAACGAGCCTTTTGACAACAGCCTGAGTATGTCAAGGAAATCCAAATGAAGACAAGCAGGGGTCTGGCCGTCATAGCTCCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATaa                            <  2:40657/214‑1 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTTAATGCGCCCAACCTCAACAATCAGCATTAGCTGGCTACTGGTTAATGACCCCTAAAGTTAAAAc                         >  1:515147/1‑229 (MQ=39)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           <  1:8448/83‑1 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           >  2:8448/1‑83 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCTACTGTTTACTAACCCCTAATGATAAAAc                         >  1:672384/1‑229 (MQ=255)
                                                                                                                                                                                              ttAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCagag                    >  1:482828/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAGCa                                                                      >  1:430393/1‑177 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGAGTCACTAATTACAGCCACATTACTGCTGACTGCATGCGCCCAACCTGACCAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCACAGGa                  >  1:1628/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTATGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCATGCGCCCAACCTGAACAATCCTCATTAGCTGGGGACTGGTTATTGACCCCTAAAGATaaaaacctgtc                   >  2:535208/1‑221 (MQ=255)
                                                                                                                                                                                                   gtgttttGGGGGCCGGAAAGGGCGTCCCGCCCGCAGCCGGCAAAAACAACAAACCGTTGTTCAACAAACTGTGTAGCCCAAGAAACCCCCAAGGAAGACCAGCGGGTGGCTGGTCGTGATAGTGTCACTAATTTCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:515147/226‑1 (MQ=255)
                                                                                                                                                                                                   ctaTTTTATAGGCGGAATAAGGGTTTCCGCCCCCATCGGGAAAAAAAAACGAACCCTTTGCCAAAAAACTGGTTAGCCAAAGAAACCCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTTTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:594493/227‑1 (MQ=255)
                                                                                                                                                                                                               ccGGAAAAGGGTTTCCCGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGATTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGTTCGTCATAGTGTCATTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCAtt   <  1:280876/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCag                      >  1:551446/1‑209 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCCTGCGCCCAACCTGAACAATCATCATTAGCTGGCGCCCGGTTCCTGACCCCTAAACATAAAACAAGAGGATTAAAgggaaaagttg  >  2:648586/1‑221 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  NC_000913/715128‑715562

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