Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,061,235 A→C Y24S (TAT→TCT)  yihM → putative sugar phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,061,2350AC98.7% 282.7 / NA 79Y24S (TAT→TCT) yihMputative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (47/31);  minor base T (1/0);  total (48/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACGAAGCC  >  NC_000913/4061008‑4061467
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
ccGGTCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCt                                                                                                                                                                                                                                     <  2:178726/233‑1 (MQ=255)
ccGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTc                                                                                                                                                                                                                                  <  1:382992/236‑1 (MQ=255)
  gggCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTTGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGGTTTTTGGTGACAATAAATAAAGCCAGAAAGATTTTACAAAGGTTCGGTACTCTTCCTCTTTATTTACAAGCATCTGCCTTTCAt                                                                                                                                                                                                                                >  2:45462/1‑236 (MQ=255)
   ggCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCA‑‑GCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTa                                                                                                                                                                                                                             >  1:429618/1‑236 (MQ=255)
                  tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTTAGGCCGTAGACCTCATCATTAATCTTTAGGGGTTTTATGGTGAAAATAAAAAACCCAAGAAAAATTTTACAAAGGTTCGACACCTTTCCTCTTTATTTACATGCTTTTGCCTTTCATTTTAATTTGCGGCTgg                                                                                                                                                                                                                >  2:571719/1‑236 (MQ=255)
                    tCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTTAATATTCGGCTGGaa                                                                                                                                                                                                              >  1:182492/1‑236 (MQ=255)
                                tGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTg                                                                                                                                                                                                     >  1:498235/1‑233 (MQ=255)
                                tGCGGGTATGCCAGCGATCTTAATTGAAAAGCA‑‑GCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                                                                                                                                                                                                 <  2:429618/235‑1 (MQ=255)
                                       aTGCCAGCGATCTTAATTGAAAAGCATCGCGACGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                            <  2:677086/235‑1 (MQ=255)
                                                            aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGTAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                      <  2:405972/236‑1 (MQ=255)
                                                            aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                      >  2:578750/1‑236 (MQ=255)
                                                                                   cAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAg                                                                                                                                                                           <  1:288491/208‑1 (MQ=255)
                                                                                   cAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAg                                                                                                                                                                           >  2:288491/1‑208 (MQ=255)
                                                                                       atTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCAt                                                                                                                                            <  2:82837/235‑1 (MQ=255)
                                                                                           tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGaa                                                                                                                                                                                                              >  1:326579/1‑165 (MQ=255)
                                                                                           tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGaa                                                                                                                                                                                                              <  2:326579/165‑1 (MQ=255)
                                                                                           tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGCGATGCCTGCTGATTTACTTGATATCGCAAGTGAACATAATCTACATGGCTTCAAATTCCCTGTtct                                                                                                                                       >  2:465698/1‑236 (MQ=255)
                                                                                            ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGCAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGCAAATAATCTACGTGGCGTGCAACTCCATGTtctc                                                                                                                                      >  2:636036/1‑236 (MQ=255)
                                                                                                     tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAg                                                                                                                             >  2:34757/1‑236 (MQ=255)
                                                                                                         gAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGtt                                                                                                                         <  1:571719/236‑1 (MQ=255)
                                                                                                               tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTtctc                                                                                                                   >  1:27274/1‑236 (MQ=255)
                                                                                                                 gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCtt                                                                                                                 >  2:46330/1‑236 (MQ=255)
                                                                                                                    aTTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGt                                                                                                              >  1:211111/1‑236 (MQ=255)
                                                                                                                      ttGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTaa                                                                                                            >  1:432253/1‑236 (MQ=255)
                                                                                                                            ccGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      <  1:546138/236‑1 (MQ=255)
                                                                                                                               tAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGa                                                                                                   >  2:363833/1‑236 (MQ=255)
                                                                                                                                      catcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATaa                                                                                                >  1:552614/1‑232 (MQ=255)
                                                                                                                                      catcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATaa                                                                                                >  1:495840/1‑232 (MQ=255)
                                                                                                                                       atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCaa                                                                                                                                                                            >  1:54515/1‑155 (MQ=255)
                                                                                                                                       atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCaa                                                                                                                                                                            <  2:54515/155‑1 (MQ=255)
                                                                                                                                              aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      <  2:706024/218‑1 (MQ=255)
                                                                                                                                              aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      >  1:706024/1‑218 (MQ=255)
                                                                                                                                                 tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCct                                                                                                                                                                                                  >  1:658097/1‑123 (MQ=255)
                                                                                                                                                 tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCct                                                                                                                                                                                                  <  2:658097/123‑1 (MQ=255)
                                                                                                                                                 tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      >  2:319157/1‑215 (MQ=255)
                                                                                                                                                 tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      <  1:319157/215‑1 (MQ=255)
                                                                                                                                                  gTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCttt                                                                                <  2:495840/236‑1 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                            >  1:611326/1‑127 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                            <  2:611326/127‑1 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      <  1:139640/213‑1 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      >  2:139640/1‑213 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTg                                                                               >  2:65707/1‑236 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTg                                                                               >  1:432397/1‑236 (MQ=255)
                                                                                                                                                   tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTg                                                                               >  2:723194/1‑236 (MQ=255)
                                                                                                                                                         ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATTTACGAGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGAt                                                                          >  1:496004/1‑235 (MQ=255)
                                                                                                                                                           ttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                                                                                                                                                                      <  2:625084/109‑1 (MQ=255)
                                                                                                                                                           ttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                                                                                                                                                                      >  1:625084/1‑109 (MQ=255)
                                                                                                                                                           ttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGATTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGGAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATaa                                                                        >  1:667327/1‑235 (MQ=255)
                                                                                                                                                                gTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCccg                                                                  <  1:636036/236‑1 (MQ=255)
                                                                                                                                                                 tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      <  1:180254/199‑1 (MQ=255)
                                                                                                                                                                 tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                      >  2:180254/1‑199 (MQ=255)
                                                                                                                                                                 tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAg                                                                      >  1:221029/1‑231 (MQ=255)
                                                                                                                                                                 tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAg                                                                      >  1:410622/1‑231 (MQ=255)
                                                                                                                                                                 tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGctc                                                                   >  2:201721/1‑232 (MQ=255)
                                                                                                                                                                           taaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTg                                                           <  1:65707/232‑1 (MQ=255)
                                                                                                                                                                             aCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCtt                                                     >  1:268360/1‑236 (MQ=255)
                                                                                                                                                                              cGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATCTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAg                                                                      <  1:698028/218‑1 (MQ=255)
                                                                                                                                                                              cGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATCTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAg                                                                      >  2:698028/1‑218 (MQ=255)
                                                                                                                                                                               gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACg                                                                                                    >  2:458210/1‑187 (MQ=255)
                                                                                                                                                                               gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACg                                                                                                    <  1:458210/187‑1 (MQ=255)
                                                                                                                                                                                    aaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTc                                              <  2:211111/236‑1 (MQ=255)
                                                                                                                                                                                            tACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTc                                              <  2:410622/228‑1 (MQ=255)
                                                                                                                                                                                                  gtgtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCa                                  >  2:524030/1‑234 (MQ=255)
                                                                                                                                                                                                   tgtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTg                                                                               >  1:467064/1‑188 (MQ=255)
                                                                                                                                                                                                   tgtCGACACTCTTCCTATGTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTg                                                                               <  2:467064/188‑1 (MQ=255)
                                                                                                                                                                                                    gtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGTTAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACAAtcgc                              >  1:659702/1‑236 (MQ=255)
                                                                                                                                                                                                    gtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAgcgc                              <  2:27274/236‑1 (MQ=255)
                                                                                                                                                                                                     gcGACCCTTTTCCTCATGATTTACATGCTTCGGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGcc                             <  2:498235/235‑1 (MQ=255)
                                                                                                                                                                                                       gACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCt                            <  1:363833/235‑1 (MQ=255)
                                                                                                                                                                                                            tctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGAt                       >  1:460959/1‑235 (MQ=255)
                                                                                                                                                                                                             ctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATaa                     >  1:565759/1‑236 (MQ=255)
                                                                                                                                                                                                              ttCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAg                    >  2:65024/1‑236 (MQ=255)
                                                                                                                                                                                                                          taCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCGATTGa         >  2:360804/1‑235 (MQ=255)
                                                                                                                                                                                                                             aTGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACGaa     >  1:698282/1‑236 (MQ=255)
                                                                                                                                                                                                                              tGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGAc        >  1:438246/1‑232 (MQ=255)
                                                                                                                                                                                                                                cTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCa                          <  1:55905/212‑1 (MQ=255)
                                                                                                                                                                                                                                cTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCa                          >  2:55905/1‑212 (MQ=255)
                                                                                                                                                                                                                                cTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACGAAGcc  <  1:578750/236‑1 (MQ=255)
                                                                                                                                                                                                                                cTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACGAAGc   >  1:212802/1‑235 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
CCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACGAAGCC  >  NC_000913/4061008‑4061467

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