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,9920AT95.0% 110.9 / NA 40N271Y (AAT→TAT) oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (20/18);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCA  >  NC_000913/1301751‑1302234
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ttCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGGGAAATCGACATGACGAATAACAGc                                                                                                                                                                                                                                             >  2:484777/1‑250 (MQ=255)
   aCCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTTCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGc                                                                                                                                                                                                                                         <  1:54505/251‑1 (MQ=255)
         gCCAAAAGCCGCTAACGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCACCGGTGCCTATACCTTAAAAGAGTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATc                                                                                                                                                                                                                                    <  2:309022/250‑1 (MQ=255)
                gCCGCTAACGAGAAATTCGGCGAAAAATTGACCCAGCCGGTTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCATCCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAAGTGtt                                                                                                                                                                                                                            <  2:165488/251‑1 (MQ=255)
                     aaTCGAGAAATTCGGCTAAACATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGACGTACACGAACTAATCGTTCGTTAACGCATCCCGACCTACTGTAATAATGCGAAAACCGTTATTAACCAGGTAACCTAGTTTCCTATTGCTTCTGAAGTTACCGATGTCACCCTCTAACGTAGTGGTTAAATCTACATTTCTTATAACAGCATTCCGATAGAATTGTTCCAg                                                                                                                                                                                                                        <  2:460123/249‑1 (MQ=255)
                                      aacaaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTACCCTATTTGACTATTGATTCTGAAGTTACCGATGTCACCCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGACTTGTTCCAGAAGCTGAACAAag                                                                                                                                                                                                           >  2:123569/4‑246 (MQ=255)
                                          agggACCCAGCCGGTTAATATCGTCACTAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCGTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACTTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                  <  2:393656/249‑1 (MQ=255)
                                          aTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTACCCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACATCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                  >  1:510740/1‑251 (MQ=255)
                                            ggACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGa                                                                                                                                                                                                 <  1:365469/250‑1 (MQ=255)
                                                ccAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTTGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGATCCCGGACGAAg                                                                                                                                                                                            <  1:228071/251‑1 (MQ=255)
                                                  aGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGTGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGtt                                                                                                                                                                                          <  1:486283/251‑1 (MQ=255)
                                                   gCCTGGTAATATCGTCACCAACGGTGACTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGAAACACGACCTACTGGAACAACGAGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTCCCGATGTCAACCGCTACCGTAGAGGTGACATCAACATGAAGAATAACAGCATGACGGTCGAATTGTTCAAGAAGCTGAAAAAGGGGATCCCGGAAGAAGtt                                                                                                                                                                                          >  2:307067/1‑250 (MQ=255)
                                                         tAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGACCGAATCGTTCTTGAAAGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAAAGATCCCGGCCGAAGTTCACGTTg                                                                                                                                                                                   >  2:552586/1‑251 (MQ=255)
                                                             atCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                                                                                                    <  1:447727/246‑1 (MQ=255)
                                                             atCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                                                                                                    >  2:447727/1‑246 (MQ=255)
                                                               cGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCac                                                                                                                                                                             <  1:484777/251‑2 (MQ=255)
                                                                             ccTATACCTTAACAGATTGGGTCGTAAACGAACGACTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGAAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCACGGACGAAGTTAACGTTGATCCATACCTGTGCACTTat                                                                                                                                                              >  2:560968/1‑250 (MQ=255)
                                                                                           gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTa                                                                                                                                                 >  2:111084/1‑250 (MQ=255)
                                                                                                ggTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAAcc                                                                                                                                           >  1:197896/1‑251 (MQ=255)
                                                                                                    gTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAAGTAACAACCAGaa                                                                                                                                       <  1:123569/251‑1 (MQ=255)
                                                                                                          gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGCAGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAAAAACCAGAAACCGcc                                                                                                                                 >  2:395004/1‑251 (MQ=255)
                                                                                                          gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGcc                                                                                                                                 >  1:112250/1‑251 (MQ=255)
                                                                                                                     tCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtg                                                                                                                      >  1:216235/1‑251 (MQ=255)
                                                                                                                           aCGCAGCCCGACCTACTGGACCAACGCGAAAACCGTTATTAACCCGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCCCGTTGa                                                                                                                                                                                  >  2:215334/1‑186 (MQ=255)
                                                                                                                           aCGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTGCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGa                                                                                                                                                                                  <  1:215334/186‑1 (MQ=255)
                                                                                                                            cGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGGGCACTTACTATTACGAAATTAACAACAAGAAACCGCCATTCAACGAtgtgcgtgt                                                                                                               >  1:86596/1‑251 (MQ=255)
                                                                                                                                    gACCTCCTGGAAAAACGAGAAAACCGTTATTAACCAGGTAACCTATTTGACTATGTCTTCTGAAGTTACCGATGTCACCCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGAATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAACCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg                                                                                                       <  2:136436/251‑1 (MQ=255)
                                                                                                                                            ggAACAACGCTAAAACCGTTATTATCCAGGTATCTTAGTTGCCTATTACTTCTGAAGTTACAGACGTAAACCGCTACCGTAGTGTTGAAATCTACTTGACGTATAACATCATGCCGATCAAATTGTTCCAGAATCTGAAAAAAGAGATCCCGGACGAATTTCACGTTTATCCATACCTTTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTTCGTACCGCACTGaaa                                                                                               <  1:233333/251‑1 (MQ=255)
                                                                                                                                                 aacGCGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTAACGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACAAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGCAACTAgg                                                                                          >  2:67312/1‑251 (MQ=255)
                                                                                                                                                         aaCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGGGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTGGATCCATACCTGTGCTCTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGAcc                                                                                  <  1:67312/251‑1 (MQ=255)
                                                                                                                                                           ccGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgc                                                                                >  1:538944/1‑251 (MQ=255)
                                                                                                                                                                                             ttCTGAAGTTAACGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCAGATCGAATTGTTCTAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATACATACATGTGCACTTACTATTACGACATTAACAACCAGAAACCGCACTTCAACGATGTGCGTGTGCGTACCGCAATGAAACTAGATATGGACCACGATATCATAGATAATAAAGTGAAAGCGCAgagc                                               >  2:487907/1‑247 (MQ=255)
                                                                                                                                                                                                     ttACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGccc                                      >  1:181069/1‑251 (MQ=255)
                                                                                                                                                                                                           aagTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGCAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCAATTCAACGATGTGCGTGGGCGTACCGCACTGAAACTAGGTATGAACCGCGATATCATTTTTAATAAAGTGAAAGCGCAGGGCAACATGCCAGCCTAt                                >  2:395534/3‑251 (MQ=255)
                                                                                                                                                                                                             tgcAACCTCTACCGTTGTGGTGAAATCGATATTGCATAAAACAGCATGCGGATCGATTTGTTCCAAAAGGAAAAAAAAGAGACCCCGGAAAAATTTCAGGTTGACCATACCCTGTCCACTTCCTATTACTAAATAAAAAACCGAAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCAGTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATg                               <  1:483380/248‑1 (MQ=255)
                                                                                                                                                                                                                  ccGCTACCGTAGTGGTGAAATCAACATGACTTAAAAAAGCATGCCGATCAAATTGTTCCAAAAGCTAAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACACCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTaca                         <  1:100059/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCa                                                                                                                            >  1:75224/1‑136 (MQ=255)
                                                                                                                                                                                                                                  gAAACCGAAATGACGTATAACAGCATGCCGATCGAATTGTTCCAAAAGCTGAAAAAAGAGATCCCGGACGAAGTCCACGTTGACCCATACCTGTGCACTTACTATTACGAAATTAACACCCAGAAACCGCCATTCa                                                                                                                            <  2:75224/136‑1 (MQ=255)
                                                                                                                                                                                                                                        gACATGACGTATAACAGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGGTACACTCCACCGTATACTGTTGgcgc   >  2:483718/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCa  >  1:261640/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
TTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCA  >  NC_000913/1301751‑1302234

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