Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,388,176 A→G 14.1% S5P (TCC→CCC)  aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,1760AG14.1% 190.9 / 3.8 77S5P (TCC→CCC) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (13/51);  new base G (8/4);  total (22/55)
Fisher's exact test for biased strand distribution p-value = 2.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGT  >  NC_000913/3387947‑3388412
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
ccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATa                                                                                                                                                                                                                                            <  1:474482/237‑1 (MQ=255)
 cGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTCCCGCCAAACGCAAATTGTGGGTTGCAATGGGGAAAAAAcc                                                                                                                                                                                                                                          >  1:70298/1‑238 (MQ=255)
 cGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGGTTTACCGCCAAGCGCAAATGTTGGGTTAGAATCGGGACAATAccc                                                                                                                                                                                                                                         >  1:1068908/1‑239 (MQ=255)
   gaaaTACAGCCAATAAATGTTCCGGTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCCCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAt                                                                                                                                                                                                                                       <  2:285007/237‑1 (MQ=255)
    aaaTACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAt                                                                                                                                                                                                                                       <  1:525686/237‑1 (MQ=255)
        accagcAATAAATGGCCCAATGATGCGCAAAACGCCACGATAGCTAATAGCGCCAGAATCCCGTTCACCCCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGa                                                                                                                                                                                                                                <  2:728986/234‑1 (MQ=255)
                tAAAGGTGCCGATGATGCGCAAAAAGCCACGTGAGCGAATAGCGCCAAAATTCGGTTCACCGCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCaa                                                                                                                                                                                                                        <  2:905185/239‑1 (MQ=255)
                   aTGGGCCTAGGATGCGCAAAAAGCCACGATAGCGAATGGCCCCAGAATCCGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACt                                                                                                                                                                                                                     <  2:1025499/239‑1 (MQ=255)
                       gCCGATGATGCCCAAAAAGCCAGGAGAGCGAATAGCCCCAGAAAACGGTTCACCTCCCGCGGCAAAGGCCGTACCGGCGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACg                                                                                                                                                                                                                 <  2:323069/239‑1 (MQ=255)
                        ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGCACAATGGTGATCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCCATGGAGAAAAGACCCATCGTGACTAACCAAAATCaacc                                                                                                                                                                                                                >  1:125249/1‑236 (MQ=255)
                           atgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCCGGGAGTGAAAGCCAACAAACAGCGCCAGGACAATGGCGGGCCCCCGTTTTATCCCCAAAAGGAGATGTTGGTGTGCAATGGGGGAAAAACACACACCGATTAACCaaaaacccccaag                                                                                                                                                                                                             >  2:138442/1‑229 (MQ=255)
                                   cAAAAAGCCACCAGAGCGAATAGCGCCAGAATCCGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGc                                                                                                                                                                                                     <  2:520678/239‑1 (MQ=255)
                                     aaagGCCACGGAAGCGAAGGGCCCCAAAATACGGTTCGCCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCTGTACCGCCCAGCGTGGCGTTTCCGGCTGGAAGTGAAAGCCAATAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCa                                                                                                                                                                                                   <  2:309397/239‑1 (MQ=255)
                                     aaaaGCCAAGGAAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCa                                                                                                                                                                                                   <  1:545308/239‑1 (MQ=255)
                                              gatgggtaaTAGCCCCAGAATGCGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGACGCTGTCAGTACCGGCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGga                                                                                                                                                                                          <  1:944194/232‑1 (MQ=255)
                                              gATTGCGAACAGCCCCAGAATCCCGTTCCCCTCCCGCAGCAAAGGCCGTCCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGga                                                                                                                                                                                          <  2:355386/239‑1 (MQ=255)
                                                   cGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGATTTTCCCCCAACCCGAAAATTTGGTTTGCGATGGGGAAAAAACCCACCGGGGTTAACAAAAATCACCCAGGCGAAGGGGCATTTTACgggaaaaa                                                                                                                                                                                     >  2:811268/1‑233 (MQ=255)
                                                   cGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTTTCCCCCAAAGAATATGTTGGTTAGCACATGGGAAAATACCCCACTTGATTTACCAAACTCCCGCAGGGGATGGGCCACTTTACGGAAGAAc                                                                                                                                                                                     >  1:143474/1‑239 (MQ=255)
                                                          cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCCTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTGACCAGTTTTATCGCCAACGGGAAATGATTGTTAGAAATAGTGAGAATAAACCTCTCGTGTAAACCAAACCCACACAGAGACGGGCCCATTTATGgaaaaaaacaaaac                                                                                                                                                                              >  2:99040/1‑231 (MQ=255)
                                                            ccAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGGGATTAACCCAACTCAACCAGGCGAGTGGGCATTTTTCGGAAGAAAGAAACCTg                                                                                                                                                                             >  1:638336/1‑238 (MQ=255)
                                                                  tccGGTTCCCCTCCCGCAGGAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGACAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcg                                                                                                                                                                      <  1:752009/237‑1 (MQ=255)
                                                                         cACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCCCAAAACGAATAAGTTGGGTAGCAATGGAGAAAATACCCATCGTGATTAACCCAAACCCCCCAAGCGATGGGCCATTTTTCGGAAGAAAGAAACCTGGCTTTCTTcgcggt                                                                                                                                                               >  1:744543/1‑237 (MQ=255)
                                                                                    ccaAAGGCCCGACCGGGGCAAACAAGCGCCGCGTTCAGTACCGCCCAGGGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCa                                                                                                                                                   <  2:293681/238‑1 (MQ=255)
                                                                                           ccGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGac                                                                                                                                              <  1:709126/237‑1 (MQ=255)
                                                                                           ccGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGac                                                                                                                                              <  1:916284/237‑1 (MQ=255)
                                                                                           ccGTACCGGCGGCAAAAATCGCCCCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGac                                                                                                                                              <  1:486635/237‑1 (MQ=255)
                                                                                            cGTCCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccac                                                                                                                                           <  1:872682/239‑1 (MQ=255)
                                                                                            cGTCCCGGCGGCAACAATCGCCGCTGGCAGGACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccac                                                                                                                                           <  1:493265/239‑1 (MQ=255)
                                                                                            cGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGTAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccac                                                                                                                                           <  1:731711/239‑1 (MQ=255)
                                                                                               cccGGGGGAAACAATCGCCGCTGGCAGTCCCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  2:129507/238‑1 (MQ=255)
                                                                                               cccGGCGGCAACAATCGGCGCTGTCAGTACCGCCCAGCGTGGCGTTGCCAGCTGGAAGTTAAAGCCAACAAACAGCGGCAGGACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  2:166017/238‑1 (MQ=255)
                                                                                               aCCGGGGGCAACCAACGCCGGTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  1:994842/239‑1 (MQ=255)
                                                                                               aCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  1:22940/239‑1 (MQ=255)
                                                                                                ccGGCGGCAACAATCGCCGCTGTCAGTTCCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  1:907201/238‑1 (MQ=255)
                                                                                                                ccgctgctcgTACCCCCCAGCGTGTGGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGGGCCAGTACAATGGCGGGCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCg                                                                                                                       <  2:492368/230‑1 (MQ=255)
                                                                                                                acagtGTCAGTCCCGCCCCGGGTGGCGTTTCCCGCTGGAAGTGAACGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGccc                                                                                                                        <  2:364/234‑1 (MQ=255)
                                                                                                                     ggcAGTACCGCCCGGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGGCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGc                                                                                                                  <  1:203800/237‑1 (MQ=255)
                                                                                                                          tACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac                                                                                                                                        <  1:1063685/212‑1 (MQ=255)
                                                                                                                             cGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGacc                                                                                                                                             <  1:820796/204‑1 (MQ=255)
                                                                                                                               cccAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGGTAGCAATGGGGAAAATACCCATCGGGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCCATGACCACCACTGTAGCAATGGTGCCCCCAGGcccagttcgt                                                                                                         >  1:364/1‑230 (MQ=255)
                                                                                                                               cccAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGGTGGGTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtct                                                                                                        <  2:516745/239‑1 (MQ=255)
                                                                                                                                ccAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAACGAATATGTTGGGTAGCAATGGAGAAAATACCCATCGTGATTAAACAAACTCCAGCAAGCGGGGGGCCCTTTTTAGGAAGAAACAACCCTGGCCTTTGTCCCGAACTTGTTTTCCCTGGGCCCCCACTTCAACGAGGGTGCCCCCAGCCCAGATGGtctc                                                                                                       >  2:324411/1‑239 (MQ=255)
                                                                                                                                  aGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAAAACCCAt                                                                                                                                                                                                                                       >  2:965512/1‑112 (MQ=255)
                                                                                                                                  aGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAt                                                                                                                                                                                                                                       <  1:965512/112‑1 (MQ=255)
                                                                                                                                              ccAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGCTTTACCCCAAAGCGAATATGTTGGTTAGGCATTGAGAAAATAACCATCCGTGTTAACCACACCCACGCCGGCGATGGGCCCTTTTTCGGCAAAGCGAATCCCTGCGTTCTTCTCAACCTTGTTTGCCTGTGCCCACCACCCAAACAGGGGTTCCCGCCGCGCAGAAGTTCCCTTGCCGGtttgctg                                                                                         >  2:641788/1‑234 (MQ=255)
                                                                                                                                                 gCTGGAAGTGAAAGCCAACAAACAGCGCCAGTCCAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCCCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                     <  1:251390/239‑1 (MQ=255)
                                                                                                                                                          gAAAACCAACAAACCGCGCCAGTACAATGGCGGTCGCCAGTTTTACTGCAAAGCGAATATGCTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGTGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACg                                                                              <  1:1074152/237‑1 (MQ=255)
                                                                                                                                                                    aaacaGCGCCCTTACAATGGGGGGCGCCAGTTTTACCTCAAAGCGAATATGTGGGTTAGCAATGGAAAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCCCCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACg                                                                 <  2:629506/239‑1 (MQ=255)
                                                                                                                                                                    aaacaGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACg                                                                 <  1:138442/239‑1 (MQ=255)
                                                                                                                                                                           cccAGTACAATGGCGGTCGCCAGTTTTAACGCAAAGGGAATATGTTGGTTAGCAATGGAGAAAATACCCCTCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTgcg                                                            <  1:1081066/236‑1 (MQ=255)
                                                                                                                                                                                tACAATGGCGGTCGCCAGTTTTTCCCCAAAACGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCCAACTCACGCCGGCGATGGGCCATTTTACGGGAGAACGAATACTGGCTTTCGTCGCACGCTTTTTTGCCCGGGAACCCCCCCTTTGCAGGGGTGGCCCCCGGCCCAATGTTCTCTTTCGGGTTTTGCGAGCCGAAAGCGACCCGGAAAAGCGTTgccgacac                                                       >  2:910544/1‑233 (MQ=255)
                                                                                                                                                                                 aCAATGGGGGTCGCCAGGTTTACCGAAAAGCGAAAAGTTTGGTTAGCAATGGGGAAAATCCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGc                                                    <  2:614655/239‑1 (MQ=255)
                                                                                                                                                                                  cAATGGCGGTGGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAAGGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCa                                                   <  1:288642/239‑1 (MQ=255)
                                                                                                                                                                                  cAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGGTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCa                                                   <  1:641788/239‑1 (MQ=255)
                                                                                                                                                                                  cAAAGGCGGTCCCCCAGTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATCCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCa                                                   <  1:749217/239‑1 (MQ=255)
                                                                                                                                                                                    aTGGCGGGAGCAAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATCCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCAcc                                                 <  1:715316/239‑1 (MQ=255)
                                                                                                                                                                                          gTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCCTTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCCGATGTTCTCTTTCTGGTTGTCGAGACGAATACGAACCGGGACACGTTTGGCAAGACGCAccccctca                                           >  2:222346/1‑234 (MQ=255)
                                                                                                                                                                                          gTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGATCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                           >  1:148484/1‑239 (MQ=255)
                                                                                                                                                                                          gTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGGGACGAATACGCACCGGAAAACCGTGCGCAAGACGCACCCATTCa                                           >  1:964215/1‑239 (MQ=255)
                                                                                                                                                                                          gTCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGATTAACCAAAATCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCCAGGGACACCCCAGTAGCAGTGGTTCCCGCGGGCCACAGTTTCTCTTTGTGGTTTTGCGGGGGAATACGAAACCGGACACGTTGTGGAAGACGAGCCCaagca                                           >  2:974534/1‑235 (MQ=255)
                                                                                                                                                                                           gcGCCAGTTTTACCCGAAAGCGAATCTGTTGGGTAGAAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTCACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCa                                           <  1:525353/237‑1 (MQ=255)
                                                                                                                                                                                               ctgTTTTGCCGCAAAGCGAATATTGTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGtt                                      <  2:1003833/237‑1 (MQ=255)
                                                                                                                                                                                               ccgTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGGTATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTTGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGtt                                      <  2:226783/237‑1 (MQ=255)
                                                                                                                                                                                               cAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCTCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGt                                       <  1:1003603/238‑1 (MQ=255)
                                                                                                                                                                                               cAGTTTTACCCCAAAGCAAAATTGTTGGTTAGGAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGtt                                      <  2:785235/239‑1 (MQ=255)
                                                                                                                                                                                                 gTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGAccacca                                                                                                                                         <  2:281568/140‑1 (MQ=255)
                                                                                                                                                                                                 gTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGAccacca                                                                                                                                         >  1:281568/1‑140 (MQ=255)
                                                                                                                                                                                                    tttCCCGCAAAGCGAATATGTTGGGTAGGAATGGAGAAAATTCCCAGCGTGATTAACCAAACCCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                     <  1:112126/236‑1 (MQ=255)
                                                                                                                                                                                                    tttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                     <  1:699562/236‑1 (MQ=255)
                                                                                                                                                                                                    tttACCGCAAAGCGAATATGGTGGGTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGGGACCGCCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTa                                     <  1:50227/236‑1 (MQ=255)
                                                                                                                                                                                                          gCAAAGGGAATATGTTTGTTGGAAATGGGAAAAATCCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAGCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtata                            <  2:143474/239‑1 (MQ=255)
                                                                                                                                                                                                                           ggTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGCCGAATACGAACCGGAACACGTTGCGCAGGACGCACCCATTACAGGTTAGAGTCTAGAGTCCCCAACCCTTTg             >  1:1055980/1‑238 (MQ=255)
                                                                                                                                                                                                                           ggTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTg             >  1:17277/1‑238 (MQ=255)
                                                                                                                                                                                                                               tgCAATGGGGAAAATACCCCGCGGTGATAATCAAGCTCCAGCGGGGTGTGGGCCATTTTACGGAAGAACGAATCCTGGCTTGCGTCGCGATCTTGTTTGCCAGTGACCACCCCTGTTGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtcgt         <  2:929784/237‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGCCCAGAGTTTCTCTTGCTGGATGTCGAGACGAGTAGCATCCGGAAACCGTTGCGCAAGCCGCCCCCTTTCAAGGTTATAGACCATAGACCGAACCCCTTTgtcggacccc    >  2:763497/1‑233 (MQ=255)
                                                                                                                                                                                                                                     tGGGGAAAATGCCCAGTGTGATTAACCAAACTCACGCGGGGGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcgcg   <  2:223231/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcgcg   <  1:324411/239‑1 (MQ=255)
                                                                                                                                                                                                                                         gagaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGt  >  2:200369/1‑238 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
CCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGT  >  NC_000913/3387947‑3388412

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