Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,385,152 T→C 5.6% intergenic (‑122/+39) proS ← / ← oprD proline‑‑tRNA ligase/basic amino acid specific porin OprD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,385,1520TC 5.6% 206.8 / 2.2 71intergenic (‑122/+39)proS/oprDproline‑‑tRNA ligase/basic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGC  >  NC_002947/1385004‑1385300
                                                                                                                                                     |                                                                                                                                                    
ggCGAGCAAATATTGCCTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGTGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGAtt                                                                                                                                                      <  3:348150/149‑1 (MQ=255)
    aGCAAATATTGACTTGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCg                                                                                                                                                  <  1:86173/149‑1 (MQ=255)
     gCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCgg                                                                                                                                                 >  3:371783/1‑149 (MQ=255)
     gCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCgg                                                                                                                                                 <  3:81696/149‑1 (MQ=255)
          tatTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCttt                                                                                                                                            >  5:59357/1‑149 (MQ=255)
          tatTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATACCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCttt                                                                                                                                            <  5:126139/149‑1 (MQ=255)
            ttGACTGTTGCGCATGGTTAGCCGTGTCGATTGCCTAAGTCGTTGAAATCGCCCCGCATTGTACGGTGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt                                                                                                                                           <  5:363856/148‑1 (MQ=255)
            ttGACTGGTGCGCATGGTTAGCCGTGTCTGATTGCCTAAGACGTTGAAATAGCCCCGCATTTTACGGGGGGGAAACGAAGGCGTACAGCATGCCCCAGGGCGTGACACTTGTCAGCCAACAAAAAGCCCGGCGTGATTGCCGGGCttt                                                                                                                                            >  1:27693/1‑148 (MQ=255)
               aCTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTctt                                                                                                                                       <  1:36987/149‑3 (MQ=255)
                cTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAg                                                                                                                                       <  4:395284/148‑1 (MQ=255)
                 tGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGt                                                                                                                                     >  1:44093/1‑149 (MQ=255)
                    tGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCg                                                                                                                                  >  8:352408/1‑149 (MQ=255)
                        cATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCgggg                                                                                                                               >  8:222977/1‑148 (MQ=255)
                         aTGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCa                                                                                                                             <  6:104978/149‑1 (MQ=255)
                            gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACt                                                                                                                           <  5:296282/148‑1 (MQ=255)
                            gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACt                                                                                                                           <  7:111396/148‑1 (MQ=255)
                             ttAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGc                                                                                                                         <  7:352408/149‑1 (MQ=255)
                                gCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCaa                                                                                                                       <  6:197273/148‑1 (MQ=255)
                                gCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAg                                                                                                                      <  3:410748/149‑1 (MQ=255)
                                  cGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGc                                                                                                                     >  4:276189/1‑148 (MQ=255)
                                       cGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCAAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCt                                                                                                               >  3:225008/1‑149 (MQ=255)
                                           ttGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTAc                                                                                                            <  1:187205/149‑1 (MQ=255)
                                           ttGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTAc                                                                                                            <  2:89896/149‑1 (MQ=255)
                                            tGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTAc                                                                                                            <  4:371783/148‑1 (MQ=255)
                                             gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                           <  2:44093/148‑1 (MQ=255)
                                               cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa                                                                                                        >  7:98715/1‑149 (MQ=255)
                                                  aGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATg                                                                                                      <  7:222977/148‑1 (MQ=255)
                                                   gACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCt                                                                                                    <  6:229984/149‑1 (MQ=255)
                                                      gTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAg                                                                                                 >  4:329914/1‑149 (MQ=255)
                                                      gTTGAAATAGCCCCGCATTGTACGGGGGGGAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                                >  1:886/1‑149 (MQ=255)
                                                       ttGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                                <  7:243553/149‑1 (MQ=255)
                                                        tGAAATAGCCCCGCATTGTACGGGGAGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTcttg                                                                                                                                      <  2:79117/110‑4 (MQ=37)
                                                        tGAAATAGCCCCGCATTGTAAGGGGAGGAAACGAAGGAGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTcttg                                                                                                                                      >  1:79117/1‑107 (MQ=25)
                                                                ccccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTcc                                                                                       <  3:167906/149‑1 (MQ=255)
                                                                ccccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTcc                                                                                       <  8:187536/149‑1 (MQ=255)
                                                                  ccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGAAAGCCAAGAAAAAGCCCGGCGTCATTGCCGGGCTTTTTCAGGTCCGGGGCACTCCAACCCAGGCTAAAAGGATGCTCAGCGGGAACTCCAc                                                                                     >  4:93609/1‑149 (MQ=255)
                                                                     tacaatgcGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGAt                                                                                  >  4:339006/8‑149 (MQ=255)
                                                                             gggggggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa                                                                          <  3:282043/149‑1 (MQ=255)
                                                                                  ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcg                                                                     >  8:150390/1‑149 (MQ=255)
                                                                                    aaaCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgtc                                                                   <  4:141266/149‑1 (MQ=255)
                                                                                      aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcga                                                                 <  7:150390/149‑1 (MQ=255)
                                                                                        gAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatc                                                               >  4:241929/1‑149 (MQ=255)
                                                                                        gAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatc                                                               <  6:59357/149‑1 (MQ=255)
                                                                                          aGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatcga                                                             >  2:233392/1‑149 (MQ=255)
                                                                                               tACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg                                                                              >  6:209763/1‑127 (MQ=255)
                                                                                               tACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg                                                                              <  5:209763/127‑1 (MQ=255)
                                                                                                   ggctGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATc                                                    <  2:194462/146‑1 (MQ=255)
                                                                                                         cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                                >  2:385009/1‑147 (MQ=255)
                                                                                                         cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                                >  5:411616/1‑147 (MQ=255)
                                                                                                                  tGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATc                                     >  4:152905/1‑149 (MQ=255)
                                                                                                                    acacTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                                   <  3:268100/149‑1 (MQ=255)
                                                                                                                      acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                                   >  4:76474/1‑147 (MQ=255)
                                                                                                                      acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                                   >  8:71808/1‑147 (MQ=255)
                                                                                                                      acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaa                                 <  3:339006/149‑1 (MQ=255)
                                                                                                                          gTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAAc                             >  3:429791/1‑149 (MQ=255)
                                                                                                                             aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTa                                                                                                             >  4:236683/1‑66 (MQ=255)
                                                                                                                             aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTa                                                                                                             <  3:236683/66‑1 (MQ=255)
                                                                                                                             aGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACgg                           >  8:131627/1‑148 (MQ=255)
                                                                                                                              gCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACgg                           >  1:64555/1‑147 (MQ=255)
                                                                                                                               gaaaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAg                        <  3:93609/147‑1 (MQ=255)
                                                                                                                               ccAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                                            >  7:180380/1‑129 (MQ=255)
                                                                                                                               ccAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                                            <  8:180380/129‑1 (MQ=255)
                                                                                                                                  aGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTg                     <  6:411616/149‑1 (MQ=255)
                                                                                                                                  aGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTg                     <  3:152905/149‑1 (MQ=255)
                                                                                                                                   gAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg                    <  3:241929/149‑1 (MQ=255)
                                                                                                                                   gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg                    >  2:350430/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc                  >  6:204280/1‑149 (MQ=255)
                                                                                                                                        aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGa               >  3:428515/1‑149 (MQ=255)
                                                                                                                                         gCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac              <  3:76474/149‑1 (MQ=255)
                                                                                                                                         gCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac              <  3:329914/149‑1 (MQ=255)
                                                                                                                                               cGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGacgcacgc        <  1:385009/149‑1 (MQ=255)
                                                                                                                                                   aTTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGaa    <  2:886/149‑1 (MQ=255)
                                                                                                                                                     tGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGc  >  1:169300/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGAGCAAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGC  >  NC_002947/1385004‑1385300

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