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

GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  NC_000913/715128‑715562
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
gttaagaaCATCGGAAGCAGATTGAAAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACCCA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCAt                                                                                                                                                                                                                <  1:96546/226‑1 (MQ=255)
 tGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAAGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGACTGATGCGTTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGGATTTAAGCGATTTTGTAAGC‑CGGATAAGGC‑GTTAACGCCGCATc                                                                                                                                                                                                               >  1:448155/1‑229 (MQ=255)
      aaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAACGATTTTGTAGCC‑CGGATAAGGC‑GTCCACGCAGCATCCGGcc                                                                                                                                                                                                          >  1:432867/1‑228 (MQ=255)
         aTCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTGAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCaaaa                                                                                                                                                                                                       <  1:571026/229‑1 (MQ=255)
           cGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAaac                                                                                                                                                                                                     >  1:651338/1‑229 (MQ=255)
             cAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCGGATGGTGCAATATATTGACTTTAAGCGTTTTTGTAGGC‑CGGATAAGGC‑GTTCACGCGCCATCCGGCCAAaacaa                                                                                                                                                                                                   >  1:473433/1‑229 (MQ=255)
                    aTTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGGTAAGGC‑GTTCACGCCGAATCCGGCAAAAAAAACGAaaac                                                                                                                                                                                            >  2:64285/1‑229 (MQ=255)
                              aagCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGCTTTTGTAGGC‑CGGTAAAGGC‑GTTCAGGCCGCATCCGGCAAAAAAAACGAAAAATTTTAcaac                                                                                                                                                                                   >  2:286820/1‑228 (MQ=38)
                               agCTGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑ATTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑AGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  1:375314/229‑1 (MQ=255)
                                agggTTTAGATTTGTTAGCGAAGTTTGTAAAACCGGGTAAACTCA‑TTTAAAAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGAAAAACACAGAATTGCCTGATGGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGGTAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  2:52135/227‑1 (MQ=38)
                                      gaga‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGAAAAGGC‑GTTCACGCCGCATCCGGCAAAAGCAACGAACACTTTGTCAACAAAccaaga                                                                                                                                                                          >  2:445466/1‑224 (MQ=255)
                                              ttAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAACTCaa                                                                                                                                                                   >  1:256011/1‑229 (MQ=255)
                                                 ccGAATTTCTCAAACACGCTTAAA‑ACATTTTAAAAAAAAAAGGGCCGGCGTAAGATCCACATTTTTTACTTAGGACAAACACAGAAATGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGTAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGa                                                                                                                                                                <  2:568844/228‑1 (MQ=21)
                                                        tCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCCGGCCTACGAGGATGGTGCAATATATTGAATTTAACCGATTTTGTAGTC‑CGGCTAAGGC‑GTTCACGCCCCATCCTGCAAAAACAACAAACACTTTGTCACCAACCTGAgaggccacgggaaatccaa                                                                                                                                                         >  2:611225/1‑211 (MQ=255)
                                                           gAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACTCTTAGCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATg                                                                                                                                                      <  1:624072/229‑1 (MQ=255)
                                                               caaggaggtgtttaagagaacaaaaaaaGGGGGGGCGTAAGATCGCCCTATAGTACGTAGGACAAATGTAGAATTGCCTGATGCGCTACGCTTAGCAGGCCGACGAGGATGGTGCAATATATACAATTTAAGCGATTATGTAGGC‑CGGATAAGGC‑GTTCAAGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAAGCCCAATGAAg                                                                                                                                                   <  2:8900/208‑1 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCTGCAAAAACACCGAACACTTTGTACAAAAACTGAGTAGATCAACGAAAACCAAATaaaaaac                                                                                                                                                >  2:214086/1‑225 (MQ=37)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGa                                                                                                                                                  >  1:346693/1‑227 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAATAGCCCAAGGAAAACCCAATGAAAAca                                                                                                                                                >  2:280876/1‑228 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTAAGTAGCCCAAGGAAATCCCAAAGAAAAcc                                                                                                                                                >  2:422247/1‑229 (MQ=255)
                                                                   cgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCaa                                                                                                                                              >  1:499886/1‑229 (MQ=255)
                                                                     ttaaCACA‑TTGAAAAAAAAAAGGGCTGTTCCAAGATCGCTCTTTTTTCCGGAGTCCAAACATATAATTGCTTGATTGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGTC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACa                                                                                                                                            <  1:744314/227‑1 (MQ=37)
                                                                               ttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTCCGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGCTAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGCGGAGCTCAAGGAAAACCAAAGAAAGACCAACAagagaagtg                                                                                                                                   >  2:756026/1‑224 (MQ=37)
                                                                                     aaaaaaGGGCGGTCGCAAGGTCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCa                                                                                                                             <  1:411911/229‑1 (MQ=255)
                                                                                       aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCACCAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGACTGGCCGTCCTa                                                                                                                           >  1:334582/1‑229 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       >  1:625598/1‑211 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGGATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       <  2:625597/211‑1 (MQ=255)
                                                                                              cgttCGCAAGAGCGCCTTTTTTTACGTATGACAAACACAGATTTGGCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGAAAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCAc                                                                                                                    <  2:334582/226‑1 (MQ=255)
                                                                                               ggTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGCCAAAAGGTCACt                                                                                                                   >  2:340060/1‑229 (MQ=255)
                                                                                                gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAg                                                                                                                                                   >  1:8900/1‑196 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCAACAGGAGTCTGGGCGTCATAGTGTCACTaa                                                                                                                 >  1:311813/1‑229 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCAATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACCAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATATGGTCACTac                                                                                                                 >  1:338225/1‑228 (MQ=255)
                                                                                                         tccgctttttttACATTTGACAAACAAAGAATTGCTGGATGCGCTGCGCTTATGAGGCATACGAGGAGGGGGAAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGAAAAAGCAACGAACACTTTGCCAACAAACTGAGTAGCCCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGc                                                                                                          <  1:445466/222‑1 (MQ=37)
                                                                                                             ccTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTTGAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACCGAGTAGCCCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGGCATAgtggccccaataaaagccaaca                                                                                                     >  2:563831/1‑210 (MQ=25)
                                                                                                                 ttttACGTATGACAAACACAGAATTGCATGATGCGCTACGCTTATCAGGCCTACGAGGTTGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGCTAAGGC‑GTTCAGGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGAGTCGCTCAAGGAAATCCCAATGACAACCACCAGGAGTATGGTCGTCATAGTGCAAGTACTTACCGcaccatgat                                                                                                  >  1:315891/1‑220 (MQ=25)
                                                                                                                           gACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCACCAGGAGTCTGGTCGTCATAGTGTAACAAATTACCGCCACATTACTGCTGACttc                                                                                        >  1:205877/1‑226 (MQ=255)
                                                                                                                                            cTGATGGGCTACGTTTACCGGGCCTAAGGGGATGGTGCAATTAATGGAATTAAAGCTATTTTGAAGGC‑CGGAAAAGGC‑GTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      <  1:64285/229‑1 (MQ=255)
                                                                                                                                               aTGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGACGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      >  1:543375/1‑226 (MQ=255)
                                                                                                                                                tGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGTCATACTGTCACTAATTACCGCCCCATTACTGCTGACTGCATGCGCCCAACCTGAACAatc                                                                   >  1:116748/1‑228 (MQ=255)
                                                                                                                                                   ccTACGTTTAACAGGCCTACGGGGATGGTGCAATATATGAAATTTAAGCGATTTTGTAGGC‑CGTAAAAGGC‑TTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAAtcatca                                                               <  1:422247/228‑1 (MQ=255)
                                                                                                                                                    attacctgacgaggcactcGAGGGTGGTGCAATATAGTAAATTTAAGCAATTTTTGAGGC‑CGAATAAGGC‑GTTCCAGCCGGATCCGGCAAAAACAACGAACACTTTGCCAACAAACTGAGTAGCTAAATGAAATCCCACTGAAGCCCAACAGGATGCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATcatcat                                                              <  1:286820/211‑1 (MQ=255)
                                                                                                                                                            aTCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGCCTGCATGCGCCCAACCTGAAAAATCACCATTAGCTg                                                        >  2:481984/1‑227 (MQ=255)
                                                                                                                                                                ggCCTACGGGGGGGGCGAAATATATTAATTTAAAGCGATTTGGTAGCC‑CGGAAAAGGT‑GTACACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  <  1:340060/229‑1 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  >  1:133682/1‑228 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGACTGGTCGTCATAGTGTCACTAATTACCGCCACATAACTTCTGACTGCATCCGCCCAACCTGAACAATCATCATAAACTGGCGAcg                                                  >  2:655613/1‑227 (MQ=255)
                                                                                                                                                                    tACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGtt                                              >  1:594494/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTAc                                            >  1:459613/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTACACAACCATCATTAGCTGGCGACTGGTTAc                                            >  1:407763/1‑229 (MQ=255)
                                                                                                                                                                       ggagagggtaaaaTATATTAAATTAAATAAATTTTGTAGGC‑CGAAAAGGGC‑CTCCACCACGCATCGGGAAAAAAAAACGAACACTTGAGCAACAAACTGAGTAGCTCAAGGAAACCCCAAGGAACACCAACAGGAGTCTGGTCGTCATAGTGTATCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACt                                           <  2:499886/218‑1 (MQ=255)
                                                                                                                                                                              gTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GGTCACGCCGCATCCGACAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTACCTGGCGACTGGTTACTGAcccc                                     >  1:629757/1‑228 (MQ=255)
                                                                                                                                                                               tacAAGATATTGAATTTAAGC‑AATTTGGAGGCGCGAAAAAGGT‑TTGCACGCAGCAGCGGGAAAAAACAACAAACCTTTTGCAAACAAACGGAGGAGCTAAAGGAAAGCCCAATGAAACCCAACGGGAGTCGGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:214086/227‑1 (MQ=255)
                                                                                                                                                                                gaaaataatTGAATTTAGGCGATTGTAAAGCC‑CGAAAAAGGC‑GTTCCCGCCGCATCCGGCAAAAACAACGACCACTTTGTAAACAAGCTGAGTAGCTCAAGGAAATCCCAATGAAGCCCACCAGGGGTCTGGTCGTCATGTTTTCCCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:535208/221‑1 (MQ=255)
                                                                                                                                                                                   atatatTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAgc                               >  1:371009/1‑228 (MQ=255)
                                                                                                                                                                                     atatTTAATTTAAGCGATTTTGTAGGC‑CGAATAAGGC‑GTTAACTCCGCAACCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTGGGCATCGTGTCACTAATTACCGCCACATTCCAGCTGGCTCCAGGCGCCCAACCTGAACTAGCATCAGTAGCTAGCGCCCGGTAACTGACCCCAAAAGaaa                             >  2:690943/1‑227 (MQ=255)
                                                                                                                                                                                       atataaatataacggttttGTAGGC‑GCGATAAGGCATTTCCCCCCACCACCAGAAAAAACAACGAGCCTTTTGACAACAGCCTGAGTATGTCAAGGAAATCCAAATGAAGACAAGCAGGGGTCTGGCCGTCATAGCTCCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATaa                            <  2:40657/214‑1 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCACAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTTAATGCGCCCAACCTCAACAATCAGCATTAGCTGGCTACTGGTTAATGACCCCTAAAGTTAAAAc                         >  1:515147/1‑229 (MQ=39)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           <  1:8448/83‑1 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           >  2:8448/1‑83 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCTACTGTTTACTAACCCCTAATGATAAAAc                         >  1:672384/1‑229 (MQ=255)
                                                                                                                                                                                              ttAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCagag                    >  1:482828/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAGCa                                                                      >  1:430393/1‑177 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGAGTCACTAATTACAGCCACATTACTGCTGACTGCATGCGCCCAACCTGACCAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCACAGGa                  >  1:1628/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTATGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCATGCGCCCAACCTGAACAATCCTCATTAGCTGGGGACTGGTTATTGACCCCTAAAGATaaaaacctgtc                   >  2:535208/1‑221 (MQ=255)
                                                                                                                                                                                                   gtgttttGGGGGC‑CGGAAAGGGC‑GTCCCGCCCGCAGCCGGCAAAAACAACAAACCGTTGTTCAACAAACTGTGTAGCCCAAGAAACCCCCAAGGAAGACCAGCGGGTGGCTGGTCGTGATAGTGTCACTAATTTCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:515147/226‑1 (MQ=255)
                                                                                                                                                                                                   ctaTTTTATAGGC‑GGAATAAGGG‑TTTCCGCCCCCATCGGGAAAAAAAAACGAACCCTTTGCCAAAAAACTGGTTAGCCAAAGAAACCCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTTTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:594493/227‑1 (MQ=255)
                                                                                                                                                                                                               c‑cGGAAAAGGG‑TTTCCCGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGATTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGTTCGTCATAGTGTCATTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCAtt   <  1:280876/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCag                      >  1:551446/1‑209 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCCTGCGCCCAACCTGAACAATCATCATTAGCTGGCGCCCGGTTCCTGACCCCTAAACATAAAACAAGAGGATTAAAgggaaaagttg  >  2:648586/1‑221 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGC‑GTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  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: ‑