Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,301,992 A→T N271Y (AAT→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,9920AT98.5% 221.0 / NA 68N271Y (AAT→TAT) oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (32/35);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAA  >  NC_000913/1301753‑1302235
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
cACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAATAATGGTCCCAGCCTGTTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTTAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCTAAAACCTTTATTAAGCAGGTAACCTATTTGCCGATTGCTTCTGAAGGTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCagg                                                                                                                                                                                                                                          <  2:633335/251‑3 (MQ=255)
    ggTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAGCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCg                                                                                                                                                                                                                                       >  1:81896/1‑250 (MQ=255)
     gTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAATGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGAt                                                                                                                                                                                                                                     >  1:145938/1‑251 (MQ=255)
       gcaaaaaGCCGCTATCGAGAAATGCGGCGAAAAATGGACCCAGCCTGGTAATAGCGTCACCAACGGTGCCTAGACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAAAAACGCGAAAACCGTTATTAACAAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGGGAAATCGACATGACGTATAACAGCATGCCGATCg                                                                                                                                                                                                                                   <  2:273296/248‑1 (MQ=255)
                cGCTATCGAGAAATTCGTCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTcc                                                                                                                                                                                                                          <  1:686091/251‑1 (MQ=255)
                      cGAGAAATTCGGCGAAAAATGGTCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGc                                                                                                                                                                                                                    <  1:606370/251‑1 (MQ=255)
                      cGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAg                                                                                                                                                                                                                     <  1:534279/250‑1 (MQ=255)
                          aaaTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGa                                                                                                                                                                                                                 >  1:578865/1‑250 (MQ=255)
                              tCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGAATGCCGCTCGAATTGTTCCAGAAGCTGaacaa                                                                                                                                                                                                             >  2:556867/1‑247 (MQ=255)
                              tCGGCGAAAAATGGACCCAGCCTGGTAATATAGTCCCCAGCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCCACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGTATTGTTCCAGAAGCTGaaaaa                                                                                                                                                                                                             <  2:81896/250‑1 (MQ=255)
                                                aGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGtt                                                                                                                                                                                          >  1:352877/1‑251 (MQ=255)
                                                     ggTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGTTCTTAAACGAACGAATCGTTCTTGAACGCAGCCCGAGCTACTGGAACAACGCGAAAACCGTTATTAACCATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCAc                                                                                                                                                                                       <  1:663062/249‑1 (MQ=255)
                                                        aaTATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACATATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGa                                                                                                                                                                                  <  1:351104/251‑1 (MQ=255)
                                                         atatCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAt                                                                                                                                                                                 <  1:538363/251‑1 (MQ=255)
                                                           atCGTCAACAACGGTGACTATACCTTAAAAGCTTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCAGACCTACTGGAACAACGCGAAAACAGTTCTTAACCAGGTAACATATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGAAATGACGTATAACAGCATGCCGATCGAGTTGTTCCAGCAGCTGAAAAAAGAGATCCCGGACGAAGttatagttgatac                                                                                                                                                                               >  2:44179/1‑240 (MQ=255)
                                                                 aCCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                  <  2:276589/226‑1 (MQ=255)
                                                                 aCCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                  >  1:276589/1‑226 (MQ=255)
                                                                  ccAACGGTGCCTATAACTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCGTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATGAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGCCGAAGTTCACGTTGATCCATACCtg                                                                                                                                                                        <  2:125865/251‑1 (MQ=255)
                                                                      cGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCa                                                                                                                                                                    <  1:676398/251‑1 (MQ=255)
                                                                        gTGCCTATACCTTAAAAGATTGGGTAGGAAACTAACGAATCGTTCTTGAACGCAGCCCGTCCTACTGGAACAACGAGAAAACCTTTATTAACCAGGTAACCTATTTGCCTAGTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAAGCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACt                                                                                                                                                                  <  2:184826/251‑1 (MQ=255)
                                                                          gCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACATACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTTATCCATACCTGTGCACtt                                                                                                                                                                 >  2:291546/1‑250 (MQ=255)
                                                                            cTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACt                                                                                                                                                              >  1:492352/1‑251 (MQ=255)
                                                                             tataCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATaa                                                                                                                                                                                                                                                 >  1:348640/1‑167 (MQ=255)
                                                                             tataCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATaa                                                                                                                                                                                                                                                 <  2:348640/167‑1 (MQ=255)
                                                                                aCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGACCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAtt                                                                                                                                                           <  1:565433/250‑1 (MQ=255)
                                                                                  cctAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACg                                                                                                                                                        <  2:662844/249‑1 (MQ=255)
                                                                                            tGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGa                                                                                                                                                                                  <  1:480594/215‑1 (MQ=255)
                                                                                            tGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACCCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGACTTGTTCAAGAAGCTGAAAAAAGAGATCCCGGACGAAGATCACGTTGa                                                                                                                                                                                  >  2:480594/1‑215 (MQ=255)
                                                                                                    aaaCGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAAc                                                                                                                                      <  1:441716/251‑1 (MQ=255)
                                                                                                    aaaCGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGaaa                                                                                                                                       <  1:324210/250‑1 (MQ=255)
                                                                                                     aacgaacgCATCGTTCTTGAACGCAGCCCGCCCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAAcc                                                                                                                                     >  2:569771/1‑251 (MQ=255)
                                                                                                                 gTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACg                                                                                                                          >  1:437730/1‑250 (MQ=255)
                                                                                                                  ttCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGAAGAAGTTCACGTTGATCCATACCTGTGCCCTTACTATTACGAAATTAACACCCAGAAACCGCCATGCAACGAt                                                                                                                        >  2:370160/1‑251 (MQ=255)
                                                                                                                                 cGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAGAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTAACCGATgtgcgtgtgcgtccc                                                                                                         >  2:286053/1‑248 (MQ=255)
                                                                                                                                       aCTGGCACAACGCGAAAACCGCTATTAACCAGGCAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATAGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCTCGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACAACCTGAACCAGCAATTCAACGACTTGCGTGTGAGTACAGCACTg                                                                                                   >  2:331585/1‑251 (MQ=255)
                                                                                                                                            aacaccGCGAAAACCGGTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTAACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAAAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                              <  1:627461/251‑1 (MQ=255)
                                                                                                                                            aacaacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACTTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTCCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACACCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                              <  1:417480/251‑1 (MQ=255)
                                                                                                                                            aacaacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATccc                                                                                                                                                                                                    <  2:14159/149‑1 (MQ=255)
                                                                                                                                            aacaacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATccc                                                                                                                                                                                                    >  1:14159/1‑149 (MQ=255)
                                                                                                                                              caacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                              >  1:396921/1‑249 (MQ=255)
                                                                                                                                               cactcgAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTAAAGCTACCGATGTCAACCGCTACCGTATTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAATTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAATTTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAgg                                                                                           <  2:100865/247‑1 (MQ=255)
                                                                                                                                                   cgAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTAAACAGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACAGCCCTGACACTAGGTATg                                                                                       >  2:632305/1‑251 (MQ=255)
                                                                                                                                                      aaaCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAaga                                                                                                                                                                                                          >  1:414953/1‑133 (MQ=255)
                                                                                                                                                      aaaCCGTTATTAACCAGGTAACCTAGTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAaga                                                                                                                                                                                                          <  2:414953/133‑1 (MQ=255)
                                                                                                                                                        aCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACAGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcg                                                                                  >  2:303860/1‑251 (MQ=255)
                                                                                                                                                           gTTATTACCCAGGTAACCTAGTTGCCTATTTCTGCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACCGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTCCGAAATTAACAACCAGAAACCGCCATTCAACGAtgt                                                                                                                      <  2:390768/212‑1 (MQ=255)
                                                                                                                                                           gTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTCGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTACACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgcg                                                                                >  2:557447/1‑250 (MQ=255)
                                                                                                                                                           gTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgt                                                                                                                      >  1:390768/1‑212 (MQ=255)
                                                                                                                                                                  aCCAGGTAACCTATTTGCTTATGGCTTATAAAGTTACCGATGTCACCCGCTACCGTAGTGGTGAAACCGACAGGACGTAAAACAGCATGCCGATCGATTTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTTCCTATTACGAAATTACCAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCAtt                                                                        <  2:578865/251‑1 (MQ=255)
                                                                                                                                                                    cAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAATTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCTTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGt                                                                      <  1:643562/251‑1 (MQ=255)
                                                                                                                                                                        aaaCCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGACACCTGAGTGGTGAAATCGACATGAAGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAAATCCCGGACGAAGTCCACGTTGATCCATACCTGTGCCCTTTCTATTACGAAATTAACAACCATAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTtaat                                                                  <  2:437730/250‑1 (MQ=255)
                                                                                                                                                                           ccTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCa                                                                                                                             >  1:72177/1‑189 (MQ=255)
                                                                                                                                                                             gagttGCCTTTTGCTGCTGAAGTTACTGATGTCAACCGCTACCGTAGTGGTAAAATCGACAGGACGTATAACATCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGAACCCGGACGAATTTCACGTTGATCCATACCTGGGCACTTCCTATTAGGAAATTAACAACCAAAATCCGCCATTCAACGATGTGCGTGTGCTTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTtaatagag                                                              <  2:235457/247‑4 (MQ=255)
                                                                                                                                                                                 ggcttATTTCTTCTGAAGTTACGGATGGCAACTGCTACCTTAGTTGTAAAATCGACAGGAGATATAACTGCACTCCGTTCGATTTTTTCCAGAATCGGAAAAAAGAGAGTCCGGACGAATTTCACGTGGATCCATACCTGTGCTCTTACTATTACGAATTTAACACCCAGAAACCGCCATGCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGaaa                                                         <  2:558029/247‑1 (MQ=255)
                                                                                                                                                                                   ccTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGGGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAgc                                                       <  1:115686/251‑1 (MQ=255)
                                                                                                                                                                                     gtgttcttctaaaggtgcagATGTCATCCGCTACCGTAGGGGTAAAATCTACATGCCGTATACAAGCAGGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGCTCCCGGACGAAGGTCACGTTGATCTATCCCTGTGCCCTTACTACTACGAAATTAAAAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGTACCGCGATATCATTGTTAATAAAGTGAAAgcg                                                      <  2:542562/231‑1 (MQ=255)
                                                                                                                                                                                      aTTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAgcg                                                      <  1:556867/249‑1 (MQ=255)
                                                                                                                                                                                       ttGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAgcgc                                                     >  1:169836/1‑249 (MQ=255)
                                                                                                                                                                                         gCTTCTGAAGTTACCGATGTCAACCGCTCCCGAAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAggg                                                 >  1:512114/1‑251 (MQ=255)
                                                                                                                                                                                                aaGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATg                                          >  2:148183/1‑251 (MQ=255)
                                                                                                                                                                                                    tACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                                                                                                    <  2:558267/109‑1 (MQ=255)
                                                                                                                                                                                                    tACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                                                                                                    >  1:558267/1‑109 (MQ=255)
                                                                                                                                                                                                     aCCGATGTCAACCGTTACCGTAGTGGTGAAATCGACGTGACGTATAACAGCATGCCTATCGAATTGTCCCAAAAGCTGAAAAAAGAGATCCCGGACGAATTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACGACCAAAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCg                                      <  2:246575/250‑1 (MQ=255)
                                                                                                                                                                                                                   cTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCAc                                                                                                     <  2:552686/173‑1 (MQ=255)
                                                                                                                                                                                                                   cTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCAc                                                                                                     >  1:552686/1‑173 (MQ=255)
                                                                                                                                                                                                                   cTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACCCTc                       >  2:633160/1‑251 (MQ=255)
                                                                                                                                                                                                                              gtgAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACt            >  1:475129/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aTGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCaa  >  1:211785/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
CACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAA  >  NC_000913/1301753‑1302235

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