Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,557,486 G→A L185L (CTG→TTG)  ddpF ← putative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,4860GA98.0% 164.9 / NA 51L185L (CTG→TTG) ddpFputative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/33);  minor base C (1/0);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACACCGGTAAGCGCGGATTAAGCGAAGAGAGCGGGTC  >  NC_000913/1557216‑1557777
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   
aaCAGCCTTGCTGCAACGTGCGGTTTCCCTGCAGATCCGTTTTACTTAATGCCCATTCTTCCTCCAGCTGTTTATCAATTGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCTTCGCCCAGTTCTACAACCCGCCCTAGATACCTCACCGCCACCCGATCGCTCACATTACGTACCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGgaga                                                                                                                                                                                                                                                                        <  2:100251/300‑1 (MQ=255)
    gCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATc                                                                                                                                                                                                                                                                    >  1:205480/1‑300 (MQ=255)
    gCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCACGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATc                                                                                                                                                                                                                                                                    >  2:697409/1‑300 (MQ=255)
     ccTTTCGGCAACATGCGGTTTCCCTGCAGATCCGTTTTACTTAATGCCCATTCCTCCTCTAGCGGTCTATCAATTGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCTCCACCCGATCGCTCATATCACCTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATcc                                                                                                                                                                                                                                                                   <  2:212530/300‑1 (MQ=255)
     ccTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATcc                                                                                                                                                                                                                                                                   <  1:179049/300‑1 (MQ=255)
          cGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACTTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAgcgc                                                                                                                                                                                                                                                              <  1:698513/300‑1 (MQ=255)
          cGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAgcgc                                                                                                                                                                                                                                                              <  1:698521/300‑1 (MQ=255)
                                        tttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTc                                                                                                                                                                                                                                <  1:260475/300‑1 (MQ=255)
                                        tttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTc                                                                                                                                                                                                                                <  2:450348/300‑1 (MQ=255)
                                                   cccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACATATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAag                                                                                                                                                                                                                     <  2:682903/300‑1 (MQ=255)
                                                                 aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTACCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCa                                                                                                                                                                                                       <  1:371738/300‑1 (MQ=255)
                                                                            aaTGGCGGGGAGGGAATCCAGCAATAATCGTGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACACCACCGCCACCCGATCGCTCATACGACCTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGTAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACg                                                                                                                                                                                            <  2:279272/300‑1 (MQ=255)
                                                                                gCGGGTAGGGAATCCAGCAATAATCGGGTGTATGGGTGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCTAGATACATCACCGCCACCCGATCGCCCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGc                                                                                                                                                                                        <  2:256457/300‑1 (MQ=255)
                                                                                               aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCTAGATACATCACCGCCACCCGATCGCTCATANGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGTAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGgag                                                                                                                                                                         <  2:621362/300‑1 (MQ=255)
                                                                                                           gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGAACTTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTATCGCGTTACCAACACAGTGAGGATCTGCGCATGAACGGAGATATCCAGCACAGAGGCTGACTACTCACGCACCATCACGCCAGGCTGCGAAGAGAGTGTTTTGGCCATGGCGATACGTTGCCgtcgcacgtaggagaacgaatgcggc                                                                                                                                                             >  2:274451/1‑275 (MQ=255)
                                                                                                                   aTGTACAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACACTGTGTGAAATCAGCACATGGGTCAGCCTGTGGTTTTCCTTTAGCGTTACCCATAAATTGCGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTCGGCTCATCAAGCACAGTCACGTCAAGCTGCGAAGAGAGTGCTCNGGCAATGGCGATACATTGCCGCTGCCCGCATGAGAAAGCAtacgacaggcggcc                                                                                                                                                     >  2:525911/1‑287 (MQ=255)
                                                                                                                                  tcgtagatgtgtataagagacaggTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGNTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGATGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACg                                                                                                                                      <  1:503247/277‑1 (MQ=255)
                                                                                                                                   aGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGa                                                                                                                                     >  1:147292/1‑300 (MQ=255)
                                                                                                                                     cACCTGCTGCGCGTCGCCCAGTCCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAGTTTCCTGTAGAGTTACCCCTAAAGTGAGGATCTGCGCCCGACCGACGATATCCAGCGCAGAGGTTGGCTCATCAAGCACCATCACGCCAAGCCACGCAGAGAGTGTCCCGTCACTGGCGATCCGTTGCCGCCGTCCGCCGGCGAACGCATGAGGCCGGCGGTCGAGGTCCTACGGACGagt                                                                                                                                   >  2:351013/1‑298 (MQ=255)
                                                                                                                                                 gTCGCCCAGTTCTACAATCTGCCCTAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGcac                                                                                                                       <  1:6335/300‑1 (MQ=255)
                                                                                                                                                      ccAGTTCTACAATCTGCCCTAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGTAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCa                                                                                                                  <  1:190671/300‑1 (MQ=255)
                                                                                                                                                         gTTCTACAATCTGCCCGAGATACATCAACGCCACCCGATCGCTCATATGACGTATCACCGAGACGTTGTGTAAAATCAGCACATAGGTAAGACCGTGGTTTTCCTATAGCTTTAGCAATACCTTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGCGGTTAGCTCATCAAGCACAATCACGTCCGGCTGCGACGGGAGTGTTCTAGCACTAGCGACACGTTACCGCTGCCCGCCGAAGAACGCATGAGGCAGGCAGTCGAGATACTCCGGACGAACACCCACTTGcccatcaaatt                                                                                                               >  2:16528/1‑290 (MQ=255)
                                                                                                                                                         gTTCTACAATCTGCCCGAGATACATCAACGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAAct                                                                                                               >  1:337731/1‑300 (MQ=255)
                                                                                                                                                           tCTACAATCTGCCCGAGATACATCACCGCAACCCGATCGCTCATATGACGTATCCACGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGACGGTCGAGATACTCCGGACGAATACCCACCTGCCCAGCCCActcc                                                                                                             >  1:358726/1‑300 (MQ=255)
                                                                                                                                                                 aTCTGCCCTAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGcc                                                                                                       <  2:146693/300‑1 (MQ=255)
                                                                                                                                                                            aTACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTcgccg                                                                                            <  1:521034/300‑1 (MQ=255)
                                                                                                                                                                              aCATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCNATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCt                                                                                          >  1:478518/1‑300 (MQ=255)
                                                                                                                                                                                  cACCGCCACCCGATCGCTCATATGACGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTGGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCtgttg                                                                                      <  2:330714/300‑1 (MQ=255)
                                                                                                                                                                                   cccaccccccGCTCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCCGTCACGTTACCAATCAATTGAGGACCTGCGCCTGAGCGGAGATATCCAGTGCAGAGGTTGGCTCCTCAAACACAATTACGTCAGGCTGCGAAGAAAGTGCTCTGGCAATGGTGATACGTTGCCGCTGCACTCCAGAGAACACATGAGGCAGGCGGTCGAAATACTCAGGACGAGTACCCACCTACCCAGCCAACTCCTCTGcaaaagcttgcggtcgttgt                                                                                     >  2:197674/8‑281 (MQ=255)
                                                                                                                                                                                     cGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCTAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTc                                                                                   <  1:441056/300‑1 (MQ=255)
                                                                                                                                                                                               tCGCTCATATGACTTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCNCTTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATANCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGc                                                                         <  2:644116/300‑1 (MQ=255)
                                                                                                                                                                                                 gCTCAGATGACTTATCACCGATACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCNAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCtt                                                                       <  2:150607/300‑1 (MQ=255)
                                                                                                                                                                                                  cTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTa                                                                      >  1:554932/1‑300 (MQ=255)
                                                                                                                                                                                                    cATATGACGTATCACCGAGACATTGTTTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGc                                                                    <  1:176857/300‑1 (MQ=255)
                                                                                                                                                                                                      tatGACGTATCAACGAGACATTGTGTGAAATTAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCATGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCAAGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACCCCGGACGAATACCCACCTGCACCGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTa                                                                  >  2:196218/1‑300 (MQ=255)
                                                                                                                                                                                                         taCGTATCACCTAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGTAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAANCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATcc                                                               <  2:166380/299‑1 (MQ=255)
                                                                                                                                                                                                              aTCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGTAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGc                                                          <  1:34612/300‑1 (MQ=255)
                                                                                                                                                                                                              aTCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGc                                                          <  1:93132/300‑1 (MQ=255)
                                                                                                                                                                                                                  ccGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGtt                                                      >  2:245621/1‑300 (MQ=255)
                                                                                                                                                                                                                   cGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTc                                                     <  1:117287/300‑1 (MQ=255)
                                                                                                                                                                                                                      gaCATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCtgt                                                  >  2:538291/1‑300 (MQ=255)
                                                                                                                                                                                                                        cATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCNAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGtga                                                <  1:511686/300‑1 (MQ=255)
                                                                                                                                                                                                                        cATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGtga                                                <  2:385989/300‑1 (MQ=255)
                                                                                                                                                                                                                          ttGTTTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCGACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGtgatg                                              <  2:248890/300‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtgAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTgatgat                                            >  1:248333/1‑300 (MQ=255)
                                                                                                                                                                                                                             tgtgAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGatgatg                                           <  1:249222/300‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACAcc                                  <  1:246955/300‑1 (MQ=255)
                                                                                                                                                                                                                                         acaTAGGTCAGCCCGTGATTTTCCTGTAGCGTTCCCAATAAATTGAGGATCTGCGCCCGAACGGGGGTATCCAGAGCAAAGGTTGGCTCATCAAGCCCAATCACGTCGGGCTGCGAAGGGAGTGCTCTGGGAATGGCGATACGTTGCCGCTACCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCAGGACAAATACCCACCTGCACAGCCAACTCCTCTGCCGACGCTCGCCGCTGTTGTTCACTACTAAGCTTAGCTATCCCGAGCGGTTCTATGATGACGCGCCACAACGGt                               >  2:283360/1‑300 (MQ=255)
                                                                                                                                                                                                                                              ggTCAGCCCGTGATTTTCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACACCGGTAAgcg                          >  1:616101/1‑300 (MQ=255)
                                                                                                                                                                                                                                                              tCCTGTAGCGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACACCGGTAAGCGCGGATTAAGCGAagag          <  1:434474/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cGTTACCAATAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACACCGGTAAGCGCGGATTAAGCGAAGAGAGCGGGTc  <  2:173257/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   
AACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCCTGAACGGAGATATCCAGCGCAGAGGTTGGCTCATCAAGCACAATCACGTCAGGCTGCGAAGAGAGTGCTCTGGCAATGGCGATACGTTGCCGCTGCCCGCCGGAGAACGCATGAGGCAGGCGGTCGAGATACTCCGGACGAATACCCACCTGCACAGCCAACTCCTCTGCCAACGCTCGCCGCTGTTGTTCACTACTACGCTTAGCTATCCAGAGCGGTTCTGTGATGATGCGCCACACCGGTAAGCGCGGATTAAGCGAAGAGAGCGGGTC  >  NC_000913/1557216‑1557777

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