Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,338 T→A 15.5% 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,3380TA15.5% 154.3 / 4.1 63intergenic (+140/‑74)pgm/ybfPphosphoglucomutase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (40/12);  new base A (2/9);  total (42/21)
Fisher's exact test for biased strand distribution p-value = 4.36e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  NC_000913/715128‑715562
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                
gttaagaaCATCGGAAGCAGATTGAAAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACCCA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCAt                                                                                                                                                                                                                <  1:96546/226‑1 (MQ=255)
 tGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAAGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGACTGATGCGTTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGGATTTAAGCGATTTTGTAAGC‑CGGATAAGGCGTTAACGCCGCATc                                                                                                                                                                                                               >  1:448155/1‑229 (MQ=255)
      aaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAACGATTTTGTAGCC‑CGGATAAGGCGTCCACGCAGCATCCGGcc                                                                                                                                                                                                          >  1:432867/1‑228 (MQ=255)
         aTCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTGAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCaaaa                                                                                                                                                                                                       <  1:571026/229‑1 (MQ=255)
           cGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAaac                                                                                                                                                                                                     >  1:651338/1‑229 (MQ=255)
             cAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCGGATGGTGCAATATATTGACTTTAAGCGTTTTTGTAGGC‑CGGATAAGGCGTTCACGCGCCATCCGGCCAAaacaa                                                                                                                                                                                                   >  1:473433/1‑229 (MQ=255)
                    aTTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGGTAAGGCGTTCACGCCGAATCCGGCAAAAAAAACGAaaac                                                                                                                                                                                            >  2:64285/1‑229 (MQ=255)
                              aagCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGCTTTTGTAGGC‑CGGTAAAGGCGTTCAGGCCGCATCCGGCAAAAAAAACGAAAAATTTTAcaac                                                                                                                                                                                   >  2:286820/1‑228 (MQ=38)
                               agCTGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑ATTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑AGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  1:375314/229‑1 (MQ=255)
                                agggTTTAGATTTGTTAGCGAAGTTTGTAAAACCGGGTAAACTCA‑TTTAAAAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGAAAAACACAGAATTGCCTGATGGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGGTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                 <  2:52135/227‑1 (MQ=38)
                                      gaga‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGAAAAGGCGTTCACGCCGCATCCGGCAAAAGCAACGAACACTTTGTCAACAAAccaaga                                                                                                                                                                          >  2:445466/1‑224 (MQ=255)
                                              ttAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAACTCaa                                                                                                                                                                   >  1:256011/1‑229 (MQ=255)
                                                 ccGAATTTCTCAAACACGCTTAAA‑ACATTTTAAAAAAAAAAGGGCCGGCGTAAGATCCACATTTTTTACTTAGGACAAACACAGAAATGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGTAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGa                                                                                                                                                                <  2:568844/228‑1 (MQ=21)
                                                           gAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACTCTTAGCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATg                                                                                                                                                      <  1:624072/229‑1 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCTGCAAAAACACCGAACACTTTGTACAAAAACTGAGTAGATCAACGAAAACCAAATaaaaaac                                                                                                                                                >  2:214086/1‑225 (MQ=37)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGa                                                                                                                                                  >  1:346693/1‑227 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAATAGCCCAAGGAAAACCCAATGAAAAca                                                                                                                                                >  2:280876/1‑228 (MQ=255)
                                                                 cgcgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTAAGTAGCCCAAGGAAATCCCAAAGAAAAcc                                                                                                                                                >  2:422247/1‑229 (MQ=255)
                                                                   cgTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCaa                                                                                                                                              >  1:499886/1‑229 (MQ=255)
                                                                     ttaaCACA‑TTGAAAAAAAAAAGGGCTGTTCCAAGATCGCTCTTTTTTCCGGAGTCCAAACATATAATTGCTTGATTGGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGTCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACa                                                                                                                                            <  1:744314/227‑1 (MQ=37)
                                                                               ttAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTCCGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGCTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGCGGAGCTCAAGGAAAACCAAAGAAAGACCAACAagagaagtg                                                                                                                                   >  2:756026/1‑224 (MQ=37)
                                                                                     aaaaaaGGGCGGTCGCAAGGTCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCa                                                                                                                             <  1:411911/229‑1 (MQ=255)
                                                                                       aaaaGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCACCAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGACTGGCCGTCCTa                                                                                                                           >  1:334582/1‑229 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       >  1:625598/1‑211 (MQ=255)
                                                                                             gCGGTCGCAAGATCGCCCTTTTTTACGGATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGt                                                                                                                                       <  2:625597/211‑1 (MQ=255)
                                                                                              cgttCGCAAGAGCGCCTTTTTTTACGTATGACAAACACAGATTTGGCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGAAAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCAc                                                                                                                    <  2:334582/226‑1 (MQ=255)
                                                                                               ggTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGCCAAAAGGTCACt                                                                                                                   >  2:340060/1‑229 (MQ=255)
                                                                                                gTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAg                                                                                                                                                   >  1:8900/1‑196 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATAAAGACCAACAGGAGTCTGGGCGTCATAGTGTCACTaa                                                                                                                 >  1:311813/1‑229 (MQ=255)
                                                                                                 tCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCAATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCACCAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATATGGTCACTac                                                                                                                 >  1:338225/1‑228 (MQ=255)
                                                                                                         tccgctttttttACATTTGACAAACAAAGAATTGCTGGATGCGCTGCGCTTATGAGGCATACGAGGAGGGGGAAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGAAAAAGCAACGAACACTTTGCCAACAAACTGAGTAGCCCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGc                                                                                                          <  1:445466/222‑1 (MQ=37)
                                                                                                                 ttttACGTATGACAAACACAGAATTGCATGATGCGCTACGCTTATCAGGCCTACGAGGTTGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGCTAAGGCGTTCAGGCCGCATCCGGCAAAAACAACGAACACTTTGTCACAAAACTGAGTCGCTCAAGGAAATCCCAATGACAACCACCAGGAGTATGGTCGTCATAGTGCAAGTACTTACCGcaccatgat                                                                                                  >  1:315891/1‑220 (MQ=25)
                                                                                                                           gACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCACCAGGAGTCTGGTCGTCATAGTGTAACAAATTACCGCCACATTACTGCTGACttc                                                                                        >  1:205877/1‑226 (MQ=255)
                                                                                                                                            cTGATGGGCTACGTTTACCGGGCCTAAGGGGATGGTGCAATTAATGGAATTAAAGCTATTTTGAAGGC‑CGGAAAAGGCGTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      <  1:64285/229‑1 (MQ=255)
                                                                                                                                               aTGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGACGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACa                                                                      >  1:543375/1‑226 (MQ=255)
                                                                                                                                                tGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGTCATACTGTCACTAATTACCGCCCCATTACTGCTGACTGCATGCGCCCAACCTGAACAatc                                                                   >  1:116748/1‑228 (MQ=255)
                                                                                                                                                   ccTACGTTTAACAGGCCTACGGGGATGGTGCAATATATGAAATTTAAGCGATTTTGTAGGC‑CGTAAAAGGCTTTCACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAAtcatca                                                               <  1:422247/228‑1 (MQ=255)
                                                                                                                                                            aTCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGCCTGCATGCGCCCAACCTGAAAAATCACCATTAGCTg                                                        >  2:481984/1‑227 (MQ=255)
                                                                                                                                                                ggCCTACGGGGGGGGCGAAATATATTAATTTAAAGCGATTTGGTAGCC‑CGGAAAAGGTGTACACGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  <  1:340060/229‑1 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACt                                                  >  1:133682/1‑228 (MQ=255)
                                                                                                                                                                 gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGACTGGTCGTCATAGTGTCACTAATTACCGCCACATAACTTCTGACTGCATCCGCCCAACCTGAACAATCATCATAAACTGGCGAcg                                                  >  2:655613/1‑227 (MQ=255)
                                                                                                                                                                    tACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGtt                                              >  1:594494/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTAc                                            >  1:459613/1‑229 (MQ=255)
                                                                                                                                                                      cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTACACAACCATCATTAGCTGGCGACTGGTTAc                                            >  1:407763/1‑229 (MQ=255)
                                                                                                                                                                       ggagagggtaaaaTATATTAAATTAAATAAATTTTGTAGGC‑CGAAAAGGGCCTCCACCACGCATCGGGAAAAAAAAACGAACACTTGAGCAACAAACTGAGTAGCTCAAGGAAACCCCAAGGAACACCAACAGGAGTCTGGTCGTCATAGTGTATCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACt                                           <  2:499886/218‑1 (MQ=255)
                                                                                                                                                                              gTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGGTCACGCCGCATCCGACAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCAACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTACCTGGCGACTGGTTACTGAcccc                                     >  1:629757/1‑228 (MQ=255)
                                                                                                                                                                               tacAAGATATTGAATTTAAGC‑AATTTGGAGGCGCGAAAAAGGTTTGCACGCAGCAGCGGGAAAAAACAACAAACCTTTTGCAAACAAACGGAGGAGCTAAAGGAAAGCCCAATGAAACCCAACGGGAGTCGGGTTGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:214086/227‑1 (MQ=255)
                                                                                                                                                                                gaaaataatTGAATTTAGGCGATTGTAAAGCC‑CGAAAAAGGCGTTCCCGCCGCATCCGGCAAAAACAACGACCACTTTGTAAACAAGCTGAGTAGCTCAAGGAAATCCCAATGAAGCCCACCAGGGGTCTGGTCGTCATGTTTTCCCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTa                                   <  1:535208/221‑1 (MQ=255)
                                                                                                                                                                                   atatatTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAgc                               >  1:371009/1‑228 (MQ=255)
                                                                                                                                                                                     atatTTAATTTAAGCGATTTTGTAGGC‑CGAATAAGGCGTTAACTCCGCAACCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGAAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTGGGCATCGTGTCACTAATTACCGCCACATTCCAGCTGGCTCCAGGCGCCCAACCTGAACTAGCATCAGTAGCTAGCGCCCGGTAACTGACCCCAAAAGaaa                             >  2:690943/1‑227 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCACAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTTAATGCGCCCAACCTCAACAATCAGCATTAGCTGGCTACTGGTTAATGACCCCTAAAGTTAAAAc                         >  1:515147/1‑229 (MQ=39)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           <  1:8448/83‑1 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                           >  2:8448/1‑83 (MQ=255)
                                                                                                                                                                                         tGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCTACTGTTTACTAACCCCTAATGATAAAAc                         >  1:672384/1‑229 (MQ=255)
                                                                                                                                                                                              ttAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCagag                    >  1:482828/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAGCa                                                                      >  1:430393/1‑177 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGAGTCACTAATTACAGCCACATTACTGCTGACTGCATGCGCCCAACCTGACCAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAACCCACAGGa                  >  1:1628/1‑229 (MQ=255)
                                                                                                                                                                                                aaGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTATGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCATGCGCCCAACCTGAACAATCCTCATTAGCTGGGGACTGGTTATTGACCCCTAAAGATaaaaacctgtc                   >  2:535208/1‑221 (MQ=255)
                                                                                                                                                                                                   gtgttttGGGGGC‑CGGAAAGGGCGTCCCGCCCGCAGCCGGCAAAAACAACAAACCGTTGTTCAACAAACTGTGTAGCCCAAGAAACCCCCAAGGAAGACCAGCGGGTGGCTGGTCGTGATAGTGTCACTAATTTCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:515147/226‑1 (MQ=255)
                                                                                                                                                                                                   ctaTTTTATAGGC‑GGAATAAGGGTTTCCGCCCCCATCGGGAAAAAAAAACGAACCCTTTGCCAAAAAACTGGTTAGCCAAAGAAACCCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTTTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa               <  2:594493/227‑1 (MQ=255)
                                                                                                                                                                                                               c‑cGGAAAAGGGTTTCCCGCCGCATCCGGAAAAAACAACGAACACTTTGTCAACAAACTGATTAGCTCAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGTTCGTCATAGTGTCATTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCAtt   <  1:280876/229‑1 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCag                      >  1:551446/1‑209 (MQ=255)
                                                                                                                                                                                                                 cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATACCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTCCTGACTGCCTGCGCCCAACCTGAACAATCATCATTAGCTGGCGCCCGGTTCCTGACCCCTAAACATAAAACAAGAGGATTAAAgggaaaagttg  >  2:648586/1‑221 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                
GTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGA‑TTGTTAGCGAAGTTCTGAAAAACGCGTAAACACA‑TTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGC‑CGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTG  >  NC_000913/715128‑715562

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