Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 6.0% intergenic (+73/‑32) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG 6.0% 172.9 / 2.4 54intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base G (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACGACCACAAGCGCCCCC  >  NC_000913/3156294‑3156878
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               
aTCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCt                                                                                                                                                                                                                                                                                                 >  1:50982/1‑298 (MQ=255)
     tGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCg                                                                                                                                                                                                                                                                                           >  1:187729/1‑299 (MQ=255)
           ttCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                                                                                                                                                                                                                       <  1:270810/297‑1 (MQ=255)
                    tgaGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaa                                                                                                                                                                                                                                                                             <  1:381545/298‑1 (MQ=255)
                      aGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                                         <  2:33108/284‑1 (MQ=255)
                         gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                                                                                                                                                                         <  1:548035/297‑1 (MQ=255)
                                 gccgcGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGACGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGCCTTAAGTTCCCAACACCTCCACGGAGCCTGCTCCGGTGAGTTCATATaa                                                                                                                                                                                                                                                                             >  2:225826/1‑285 (MQ=255)
                                               aaCCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAACGGAGGACCGTATg                                                                                                                                                                                                                                                               >  2:38327/1‑285 (MQ=255)
                                                  ccGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                                                                                                                                                                                                            <  2:43636/285‑1 (MQ=255)
                                                    gCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACCGttatta                                                                                                                                                                                                                                            >  1:454226/1‑299 (MQ=255)
                                                                 gCAATTAGGCGTGCCGACCCACCTCCCCCACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAGAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatt                                                                                                                                                                                                                                             <  2:307855/285‑1 (MQ=255)
                                                                         gCGTGACGACCCACCTCTCCGACGACGGTCTGGAAGGCAGATCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACGGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                                                                                                                                                                                                                     <  2:227769/285‑1 (MQ=255)
                                                                               cGACCCACCTCTCCGACTGCGGTCTGGGCGGCAGCTCCATCCCTGCATTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGGTATTAAGCTACAGGATg                                                                                                                                                                                                                                <  2:435271/284‑1 (MQ=255)
                                                                                           ccGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                                                     <  1:55964/299‑1 (MQ=255)
                                                                                               cTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAgtcccagctggggca                                                                                                                                                                                                 <  1:135518/299‑15 (MQ=255)
                                                                                                     tCTGGACGGCAGGTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCAGGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                                                                                                                                                                         <  2:165263/285‑1 (MQ=255)
                                                                                                         gACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAAGTATGGCTAATCCAACCGTTATTAAGCTACCGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                                                     >  2:7822/1‑285 (MQ=255)
                                                                                                                  tCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                                                <  1:75682/297‑1 (MQ=255)
                                                                                                                     aTCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                                                           >  1:301788/1‑299 (MQ=255)
                                                                                                                                 cTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAgg                                                                                                                                                                             <  2:79871/285‑1 (MQ=255)
                                                                                                                                  tGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAgatcggaagagc                                                                                                                                                               >  1:91532/1‑287 (MQ=255)
                                                                                                                                   gAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAt                                                                                                                                                             >  1:335554/1‑299 (MQ=255)
                                                                                                                                   gAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAt                                                                                                                                                             >  1:532118/1‑299 (MQ=255)
                                                                                                                                           tGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                                                                                                                                                                    >  2:68292/1‑284 (MQ=255)
                                                                                                                                                   cACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTACAAAAGCGTTAGaa                                                                                                                                             >  1:317225/1‑299 (MQ=255)
                                                                                                                                                        cATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTgg                                                                                                                                         <  1:256350/298‑1 (MQ=255)
                                                                                                                                                        cATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTgg                                                                                                                                         <  1:352133/298‑1 (MQ=255)
                                                                                                                                                           gACCCAACTGTGCTAAAATCATGACAGTACGTTGGATGGCAGCCGCCGTAAATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                                                                                                                                                    <  2:480308/284‑1 (MQ=255)
                                                                                                                                                            aCCCCACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTGTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCg                                                                                                                                                   <  2:454886/284‑1 (MQ=255)
                                                                                                                                                                     ggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGGTCAGACGTAAGTTAAAAAAAACAAACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGAtt                                                                                                                           >  1:313748/1‑299 (MQ=255)
                                                                                                                                                                         aaaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTa                                                                                                                                     <  2:475464/285‑1 (MQ=255)
                                                                                                                                                                           aaTCATGACATTACGTTAGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTAt                                                                                                                                    <  2:294639/284‑1 (MQ=255)
                                                                                                                                                                                   tcttccgatctaTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGc                                                                                                             <  1:440736/288‑1 (MQ=255)
                                                                                                                                                                                    aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTGCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGAt                                                                                                                            <  2:464171/283‑1 (MQ=255)
                                                                                                                                                                                         gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTcagatcggaag                                                                                                                     >  2:317809/1‑275 (MQ=255)
                                                                                                                                                                                            ggATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGa                                                                                                     >  1:377408/1‑298 (MQ=255)
                                                                                                                                                                                               tGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCcgcg                                                                                                               >  2:385383/1‑285 (MQ=255)
                                                                                                                                                                                                gTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCgg                                                                                                              >  2:238371/1‑285 (MQ=255)
                                                                                                                                                                                                  cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCCGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGcc                                                                                                            >  2:75111/1‑285 (MQ=255)
                                                                                                                                                                                                      cgccgTATATACGGAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGaagaag                                                                                            <  1:78310/297‑1 (MQ=255)
                                                                                                                                                                                                      cgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGaagaag                                                                                            <  1:522811/297‑1 (MQ=255)
                                                                                                                                                                                                                        gcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGaaa                                                                        >  1:411982/1‑299 (MQ=255)
                                                                                                                                                                                                                          ccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGaaaaa                                                                      <  1:9184/299‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCTCAAACTTTCGTTTTCGGGCCTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGCAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCACCaga                                                      >  1:129247/1‑299 (MQ=255)
                                                                                                                                                                                                                                          gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATg                                                                    >  2:318920/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                       tCCAGACTTAAGTTCACAACACCTCACCGGAGCCGGCTTCGGTGAGTTCATAGAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTtcat                                        <  2:266507/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGAGTTAGAAGTGGGTTATCGCGCGATGGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTAAACAGAGAAGCAATGTTCAGCACCACTa                                >  2:512902/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                 aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACgacgac               <  1:270412/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 aGTTCACAACACCGAACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACgacgac               <  1:86317/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATAACCACTAAGCTGTGGa                       >  2:298212/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                           acCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACGACCACAAGCGcc     <  1:363558/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAacga                   >  2:255441/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                              tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACGACCACAAGCGcccc   >  1:471272/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                              tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGCACGACGACCACAAGCGcccc   >  1:158700/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                               cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTACTGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTTGATTGATACCGCCGCGGCCTACAAGCACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCAACACTCAGCTGTGGAACGACGACCACAAGCGccccc  >  1:311552/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               
ATCACTGAAGGTTCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACGACCACAAGCGCCCCC  >  NC_000913/3156294‑3156878

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