Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,341,024 G→A 9.6% R184H (CGC→CAC)  yhbJ → adaptor protein RapZ for GlmZ/GlmY sRNA decay; glucosamine‑6‑phosphate regulated; NTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,341,0240GA 9.6% 158.2 / 5.6 63R184H (CGC→CAC) yhbJadaptor protein RapZ for GlmZ/GlmY sRNA decay; glucosamine‑6‑phosphate regulated; NTPase
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base A (2/4);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATACCTCAGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGCTATGGTTACC  >  CP009273/3340877‑3341165
                                                                                                                                                   |                                                                                                                                             
aTACCTCAGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGGCTTTGACCGgcgc                                                                                                                                              >  2:255126/1‑149 (MQ=255)
    cTCAGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGGACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCActtct                                                                                                                                          >  1:645388/1‑149 (MQ=255)
     tCAGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGcttctt                                                                                                                                         >  2:7907/1‑149 (MQ=255)
     tCAGAAATGTCCGTTCACGAGCGGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTAACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTTGCTCACTCAATGCAGACTAAGGCTTTGACGTGCGcttctt                                                                                                                                         >  1:827637/1‑149 (MQ=255)
     gcgggaaTGTCCGTTGACGGGGTGGCTGAAATGCGGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGcttctt                                                                                                                                         <  1:206133/144‑1 (MQ=255)
       agcaaTGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGc                                                                                                                                       <  2:354135/146‑1 (MQ=255)
       aGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCCAAAACGATAACCCTATCGATGCAGAATACGTCTTTGACGTGCGCTTCTTGc                                                                                                                                       >  1:398786/1‑149 (MQ=255)
       aGAAAGGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGGAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGc                                                                                                                                       <  1:854314/149‑1 (MQ=255)
        gAAATGGCCGTTCACGAGCTGGCTGAAATGCGGCGTCCCCGTCCGTTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCCTTCGGCCTCAAACACGGTACCCCTATGGATGCAGATTACGTCTTTGACGTGCGCTTCTTGcc                                                                                                                                      <  1:543068/149‑1 (MQ=255)
              tCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAAccc                                                                                                                                >  2:86326/1‑149 (MQ=255)
                 gTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCa                                                                                                                             >  1:820709/1‑149 (MQ=255)
                        aGCTGGCTGAAATGCTGCGTACCCGTCGGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGAt                                                                                                                      <  2:605911/149‑1 (MQ=255)
                             gctgAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGGATCAGTATCGATGCAGATGACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGaa                                                                                                                 >  1:650693/1‑149 (MQ=255)
                                 aaaTGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTg                                                                                                             <  1:474982/149‑1 (MQ=255)
                                    tgctgcGTACCCGTCTGCTGGGTAAACGGGAACGCGAACTCACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGt                                                                                                          <  1:687914/149‑1 (MQ=255)
                                    tgctgcGTACCAGACTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTATTTCGGCCGCAAACACGGTATCCTTATCGATGCAGATTACGGCTTTGACGTGCCCTTCTTGCCGAACCCGCACTGGGATGCGACACTGAGt                                                                                                          >  1:185666/1‑149 (MQ=255)
                                     gctgcGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTc                                                                                                         <  1:663608/149‑1 (MQ=255)
                                        gcGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCaa                                                                                                      <  1:634347/149‑1 (MQ=255)
                                              ccGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAg                                                                                                >  2:370619/1‑149 (MQ=255)
                                                  ctgctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACAAGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCt                                                                                            >  2:383189/1‑149 (MQ=255)
                                                   tgctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCtt                                                                                           <  2:65664/149‑1 (MQ=255)
                                                    gctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGTCAGGTCTTg                                                                                          <  2:78447/149‑1 (MQ=255)
                                                     ctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACGCGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGa                                                                                         >  1:555592/1‑149 (MQ=255)
                                                     ctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGa                                                                                         <  2:364189/149‑1 (MQ=255)
                                                     ctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGa                                                                                         >  2:368401/1‑149 (MQ=255)
                                                     ctgGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGa                                                                                         >  2:799752/1‑149 (MQ=255)
                                                          tAAACGTGAACGCGAACTGACCATGGTTTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAAc                                                                                    <  1:691690/149‑1 (MQ=255)
                                                          tAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAAc                                                                                    <  2:29844/149‑1 (MQ=255)
                                                             aCGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTg                                                                                 >  2:14771/1‑149 (MQ=255)
                                                                    cgcgAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGATACTGCGTCCAATGACAGGTCTTGATAAACCTGTcgccgc                                                                          <  2:172078/149‑1 (MQ=255)
                                                                    cgcgAACTGACCATGGTCGTTGCGTCCTTCTGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACTCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTcgccgc                                                                          <  1:306957/149‑1 (MQ=255)
                                                                         aCTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTcc                                                                     <  2:318605/149‑1 (MQ=255)
                                                                           tGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCCCACTGGGATCCGAAACTCCGTCCAATCACAGGGCTTGATAAACCTTTCTCCGCGTTCCTc                                                                   >  1:440501/1‑149 (MQ=255)
                                                                               ccttgTCTTTCAGTCTTTCGGCTTCAAACACGGTATCCCCATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGAcc                                                               <  2:135513/145‑1 (MQ=255)
                                                                               cATTGTCTTTTAGTCTTTCGGCTTCAAACCCGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGAcc                                                               <  1:683139/149‑1 (MQ=255)
                                                                               cATGGTCTTTGGGTCTGCCGGCCTCCAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGAcc                                                               <  2:318883/149‑1 (MQ=255)
                                                                                aTGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCg                                                              <  1:707910/149‑1 (MQ=255)
                                                                                 tGGTCTTTGAGTCTTTCGTCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTAGATAAACCTGTCGCCGCCGTCCACGACCGc                                                             >  1:557586/1‑149 (MQ=255)
                                                                                           gTCTTTCGGTTTCAAACACGGTACCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAg                                                   <  1:425572/149‑1 (MQ=255)
                                                                                              tttCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCACTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTac                                                <  2:103269/149‑1 (MQ=255)
                                                                                              tttCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCACTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTac                                                <  2:621619/149‑1 (MQ=255)
                                                                                              tgtcgccctCAAACACGGTCTCCCTATCGATGCAGTTTACGGCTTCGACGTGCACTTCTGGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATCACAGGCCTTGATACACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTac                                                <  1:277285/141‑1 (MQ=255)
                                                                                                     ttCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAAttt                                         >  1:671998/1‑149 (MQ=255)
                                                                                                     ttCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAAttt                                         >  1:547296/1‑149 (MQ=255)
                                                                                                     ttCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAAttt                                         >  2:570425/1‑149 (MQ=255)
                                                                                                     ttCAAACACGGTATCCCTATCGATGCAGATTACCTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAAttt                                         >  2:556334/1‑149 (MQ=255)
                                                                                                      tCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTCAACCAtttt                                        >  1:155240/1‑149 (MQ=255)
                                                                                                       cAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCACCACCCAGAAGTAAACCAttcac                                       >  1:162717/1‑146 (MQ=255)
                                                                                                        aaaCACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTTATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTAt                                      >  1:482262/1‑149 (MQ=255)
                                                                                                          acacGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCt                                    >  2:2212/1‑149 (MQ=255)
                                                                                                          acacGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCACTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAAACTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCt                                    >  2:283522/1‑149 (MQ=255)
                                                                                                            acGGTATCCCTTTCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTAc                                  <  1:504295/149‑1 (MQ=255)
                                                                                                            acGGTACCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTAc                                  <  2:283442/149‑1 (MQ=255)
                                                                                                                 acccctacgaaTGCTGTTTATGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGTATCCGAAACTGCGTCCAATGACAGGTCTTTATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGAc                             <  1:629301/139‑1 (MQ=255)
                                                                                                                 aTCCGTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAATCTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGAc                             <  1:647586/149‑1 (MQ=255)
                                                                                                                  ccccTTTCGATGCCGATTACGTCCTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATCACCGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGAcg                            <  1:696328/148‑1 (MQ=255)
                                                                                                                   cccTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGAcgc                           <  2:737097/149‑1 (MQ=255)
                                                                                                                     cTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGa                         <  2:180250/149‑1 (MQ=255)
                                                                                                                                  ttACGTCTTGGACGTGCACTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGc            <  2:683698/149‑1 (MQ=255)
                                                                                                                                      gTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGCCTTGATAAACCTGTCGCCGCGTTCATAGACTGACACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTACCTTGAGCTagg        >  2:812539/1‑147 (MQ=255)
                                                                                                                                       ccTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGCTATgg       >  1:518707/2‑149 (MQ=255)
                                                                                                                                          ttGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGCTATGGTTa    >  2:35644/1‑149 (MQ=255)
                                                                                                                                           tGACGTGCGCCTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGCTATGGTTAc   <  1:126058/149‑1 (MQ=255)
                                                                                                                                            gACGTGCGCTCCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAAGTTTATCTACCAGACGCGAAGCTATCTTGAGCTATGGTTAcc  <  1:745531/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ATACCTCAGAAATGTCCGTTCACGAGCTGGCTGAAATGCTGCGTACCCGTCTGCTGGGTAAACGTGAACGCGAACTGACCATGGTCTTTGAGTCTTTCGGCTTCAAACACGGTATCCCTATCGATGCAGATTACGTCTTTGACGTGCGCTTCTTGCCGAACCCGCACTGGGATCCGAAACTGCGTCCAATGACAGGTCTTGATAAACCTGTCGCCGCGTTCCTCGACCGCCACACAGAAGTACACAATTTTATCTACCAGACGCGAAGCTATCTTGAGCTATGGTTACC  >  CP009273/3340877‑3341165

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