Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,301,919 G→T Q246Y (CAG→TAT oppA → oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,301,9190GT96.0% 263.8 / NA 74Q246Y (CAG→TAToppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (36/35);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATA  >  NC_000913/1301677‑1302160
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
gatgatCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCt                                                                                                                                                                                                                                             <  1:536599/250‑1 (MQ=255)
       acacaTTAGAAGTCAACTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                                                                                                                                                                                                                     <  2:43727/251‑1 (MQ=255)
         ataTTAGAAGTTATCTTAAGTGAACCCGTTCCGTACTTCTATCAATGCCTTGTTCACCCATCAACTTCACCGGTGCCAAAAGACGCTATCGAGAACTTCGGCGAAAAATGGACCCAGCCTGGTATTATCGTCACCAACGGTGCCTATACCTTCAACGATCGGGTCGTAAACGAACGAATCGTTATTGCCCGCAGCCCGACCTGCTGGAACAACGCGAAAACCGTTATTAACTATCTAACCTATTTGCCTa                                                                                                                                                                                                                                    >  2:484968/3‑250 (MQ=255)
          caTTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTAtt                                                                                                                                                                                                                                  >  1:562202/1‑251 (MQ=255)
                    cACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGa                                                                                                                                                                                                                          >  1:205998/1‑249 (MQ=255)
                         tAAGTGAACCCGTTCCGTTCTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCCCCAACGGTGTCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCTCGGCCTACTGGAACAACGCGAAATCCGTTATTAAATATGTAACCTATTTGCCTATTGGTTCTGAAGTTAc                                                                                                                                                                                                                    <  2:92102/250‑1 (MQ=255)
                         tAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                                                                                                                                                                                                                    <  2:201667/250‑1 (MQ=255)
                                         gTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTAcc                                                                                                                                                                                                   >  1:398452/1‑251 (MQ=255)
                                           aCTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGt                                                                                                                                                                                                 >  1:184529/1‑251 (MQ=255)
                                                       tACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATAGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATc                                                                                                                                                                                     >  2:476501/1‑251 (MQ=255)
                                                        aCTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCg                                                                                                                                                                                    >  1:566619/1‑251 (MQ=255)
                                                              tCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGTTAATATCGGCTCCCACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGGTCTTGAACGCAGCCCCACCTACTAGAACAACGCGAAAACCGTTATTAACTATTTAACCTATTTGCCTATTGCTTCTGAAGTTACCGAGGTCAGCCGCTACCGTAGTGGGGAAATTGACATg                                                                                                                                                                               <  2:505686/250‑1 (MQ=255)
                                                                ccccATCAACCTCACCGGTGCCAAAAGCCGCTATCGAGAACTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGTCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTTAACGCAGCCCGACCTACTCGAACAACTCGAAAACCGTTATTAACTATGTAACCGATTTGCCTATTGCTTCTTAAGTTACCGATGTCAACCGCTACCGGAGTGGTGAAATCGACATGACg                                                                                                                                                                            <  2:525082/250‑1 (MQ=255)
                                                                aCCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACg                                                                                                                                                                            >  1:50831/1‑251 (MQ=255)
                                                                aCCCAACAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACg                                                                                                                                                                            <  1:565282/251‑1 (MQ=255)
                                                                   cATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAAt                                                                                                                                                                         <  2:132894/251‑1 (MQ=255)
                                                                      cAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATTTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAAc                                                                                                                                                                      <  1:65585/251‑1 (MQ=255)
                                                                           tCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACTGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGc                                                                                                                                                                   <  1:76785/249‑1 (MQ=255)
                                                                                      aaaaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAGCGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGTTTAAATCGACGTGTCGAATAACGGCATGCCGATCGAAt                                                                                                                                                      <  2:127810/251‑1 (MQ=255)
                                                                                      aaaaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAAt                                                                                                                                                      >  1:28550/1‑251 (MQ=255)
                                                                                           ccGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTc                                                                                                                                                 <  2:50467/251‑1 (MQ=255)
                                                                                                   gagaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATGGTTCCAGAAGCt                                                                                                                                         <  2:241496/251‑1 (MQ=255)
                                                                                                    agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                        >  1:280452/1‑251 (MQ=255)
                                                                                                    agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGAGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                        >  2:621411/1‑251 (MQ=255)
                                                                                                       aaTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaa                                                                                                                                      >  2:77801/1‑250 (MQ=255)
                                                                                                          tCGGCGAAAAATGGACCCAGCCTGGTAATAACGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaaaaa                                                                                                                                   <  2:332054/250‑1 (MQ=255)
                                                                                                                  aaaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAATCTGAAAAAAGAGATccc                                                                                                                          <  1:128640/251‑1 (MQ=255)
                                                                                                                    aTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                        <  1:438049/251‑1 (MQ=255)
                                                                                                                      ggACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGAc                                                                                                                      >  1:336597/1‑251 (MQ=255)
                                                                                                                           cAGCCTGGTAATATCGTCACCAACGGGGACTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCGGCCCGACCTACTGGAACGAGGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGCAGTTACCAATGTCAACCGCTACCGTAGTGTGGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGCCGacgt                                                                                                                 >  1:202931/1‑248 (MQ=255)
                                                                                                                                    aaTATCGGCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTg                                                                                                         <  1:353534/250‑1 (MQ=255)
                                                                                                                                          ttCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCAt                                                                                                  >  2:395750/2‑250 (MQ=255)
                                                                                                                                             tccAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCt                                                                                              >  2:202235/2‑251 (MQ=255)
                                                                                                                                               aaaCGTTGCATATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCTACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGACATCGACATGACGAATAACAGCCTGCCGATCTAATTGTTCCAGAAGCTGAAACACGAGATCACGGCCGAAGTTCACGTTGATCCATACCtgt                                                                                            >  2:590208/2‑251 (MQ=255)
                                                                                                                                                  cGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCa                                                                                         <  1:262430/251‑1 (MQ=255)
                                                                                                                                                       ccTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTAc                                                                                    <  2:170037/251‑1 (MQ=255)
                                                                                                                                                           taCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAtt                                                                                >  2:571929/1‑251 (MQ=255)
                                                                                                                                                            aCCTTAAAAGATTGGGACGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTa                                                                               <  2:223478/251‑1 (MQ=255)
                                                                                                                                                              cTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACg                                                                             >  1:197886/1‑251 (MQ=255)
                                                                                                                                                                 aaaaGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGaaa                                                                          <  1:433301/251‑1 (MQ=255)
                                                                                                                                                                    aGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAAtt                                                                        >  1:565891/1‑250 (MQ=255)
                                                                                                                                                                     gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaa                                                                      <  1:2000/251‑1 (MQ=255)
                                                                                                                                                                     gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTa                                                                       >  2:116459/1‑250 (MQ=255)
                                                                                                                                                                       ttGGGTCGTAAACGAACGAATCGTTCTTGAACGCATCCCGACCTACTGGAACAACGCCAAAACCGTTATTAACTATGTCACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTTTTCCCGTAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGCTCCATACCTGTGCACTTACTATTACGAAAttcac                                                                     >  2:167698/1‑247 (MQ=255)
                                                                                                                                                                           gTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCTCGTTGATCCATGCCTGTGCACTTACTATTACGAAATTAACAACCa                                                                <  1:476433/251‑1 (MQ=255)
                                                                                                                                                                           gTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAAcc                                                                 >  2:397200/1‑250 (MQ=255)
                                                                                                                                                                               tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGcc                                                                                                                                                              <  1:266401/154‑1 (MQ=255)
                                                                                                                                                                               tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGcc                                                                                                                                                              >  2:266401/1‑154 (MQ=255)
                                                                                                                                                                                  acgaacgaATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCg                                                         <  2:552010/251‑1 (MQ=255)
                                                                                                                                                                                   cgaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGc                                                        <  2:572816/251‑1 (MQ=255)
                                                                                                                                                                                    gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGCACAACGCGAAAACAGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCACGGACGAAGTTCACGTTGATCAATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGcc                                                       >  2:228594/1‑251 (MQ=255)
                                                                                                                                                                                    gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAt                                                                                 <  1:626483/225‑1 (MQ=255)
                                                                                                                                                                                    gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAt                                                                                 >  2:626483/1‑225 (MQ=255)
                                                                                                                                                                                      acgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAACCCGTTATTAACTATGTAACCTTTTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGTAGTGGTGAAATCAACATGACGAATAACAGCATGCCGATCGAATTGTCCCAGAAGCTGAAAAAAGAGATCCCGGCCGAAGTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCAt                                                     <  1:621411/251‑1 (MQ=255)
                                                                                                                                                                                                   gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcg                                        >  1:436287/1‑251 (MQ=255)
                                                                                                                                                                                                   gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcg                                        >  1:470330/1‑251 (MQ=255)
                                                                                                                                                                                                   gAACGCAGCCCGACCTACGGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTAAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgc                                         <  1:592935/250‑1 (MQ=255)
                                                                                                                                                                                                   gAACGCAGCCCGACCGACTGTAAACACGCGAAAACCGTTATTAACTAGGTAGCCTATTCGCCTATTGCTTCTGAAGTTCACGATGACAACCGCGACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCTGACGAAGTCCACTTTGATCCATACCAGTTCACTTCCAATAACGAAATTAACAACCATAACCCGCCATTCAACGAtgtgc                                         <  2:9330/250‑1 (MQ=255)
                                                                                                                                                                                                      cGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtgt                                     >  1:193998/1‑251 (MQ=255)
                                                                                                                                                                                                         aGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtgtgcg                                  >  1:50502/1‑251 (MQ=255)
                                                                                                                                                                                                          gCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATgtgcgtgtgcgt                                 <  2:184529/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc                               <  1:543283/251‑1 (MQ=255)
                                                                                                                                                                                                             cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAcc                              >  1:536839/1‑251 (MQ=255)
                                                                                                                                                                                                                ccTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGc                            >  1:450468/1‑250 (MQ=255)
                                                                                                                                                                                                                   aCTGGAACAACGCAAAAACAGTTATTAACTATGTACCCTTGTTGCCTATTGCTTCTGAAGTGACGGATATCAACCGGTACCGTAGTGTTGAAAACGACATGACAAATAACAGCACGCCGGTCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTGGATCCATAACTGTGCACTTTGTATTACGCAGTTAACATCCAAAAGCCGCCATTCAACGATGTGCGTGTGCGTACCGCACTg                        <  2:406765/251‑1 (MQ=255)
                                                                                                                                                                                                                        aacaacGCGAAAACCGTTATTAACTAGGAAACCTATTGGCCTATTGCTTCTGAAGTTACCGAGTTCAACCGCTACCGTAGGGGTGAAAGCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCATTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCTCTGAAACt                   <  2:193998/251‑1 (MQ=255)
                                                                                                                                                                                                                         acaacGCGAAAACCGTTATTAACTAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCACTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCGGAAAAAAGAGACCCCGGACGAATTCCACGTTGACCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTa                  <  1:590208/251‑1 (MQ=255)
                                                                                                                                                                                                                           cacGCGAAAACAGTGATTAACTATGTAACCTATTTGCCTATTGCTTCAGTAAATACCGATGTCAACCGCTACCGTAGTTGCGAAATCGACATGCCGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCTCGTTGATACAATACTGTGCACTTCCTGTCAAGAAATTCACAACCAGAAGCGGCCATTCAACTATGTGCTTGTGCGTACCGCACTAAAACTAgg                >  2:432718/2‑251 (MQ=255)
                                                                                                                                                                                                                           aacGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAgg                >  1:294479/1‑251 (MQ=255)
                                                                                                                                                                                                                            acGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGt               <  2:623171/251‑1 (MQ=255)
                                                                                                                                                                                                                              gcgAAAACCGTTATTAACTTTGTAACCTTTTTGCCTTTAGCTTCTTAAGTTCCCGATGTCAACCTCTCCCGGAGGGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCGGAAAAAAGAGATCCCAGTCGAATTTCGCGTTGATCCATGCCTGTGCCCTCTCTATTACGAAATTAACAACCAGAAACCTCCATTCAACGATGTGCGTTTTCGTACCGCACTGAAAcgagtt               <  2:489684/249‑6 (MQ=255)
                                                                                                                                                                                                                                   aaCCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTCCTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGAcc        >  1:223046/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGAGGAGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGATGATCCATACGTGTGGACTTACTATTACGAAATTAACAACTAGAAAGCGCCATTCAACGATACTCGTGTGGGAACCGCAGTAAAACTAGGGATGGACcgcg     >  2:95755/1‑250 (MQ=255)
                                                                                                                                                                                                                                          attaACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTGCCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGata  <  1:263289/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATA  >  NC_000913/1301677‑1302160

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