Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,914,857 C→A 7.5% intergenic (‑71/‑175) PP_3454 ← / → PP_3455 transcriptional regulator/multidrug RND transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,914,8570CA 7.5% 210.3 / 6.5 67intergenic (‑71/‑175)PP_3454/PP_3455transcriptional regulator/multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

GAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATA  >  NC_002947/3914709‑3915005
                                                                                                                                                    |                                                                                                                                                    
gAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTc                                                                                                                                                      >  5:25602/1‑149 (MQ=255)
gAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTc                                                                                                                                                      >  8:309845/1‑149 (MQ=255)
 aaGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTcc                                                                                                                                                     <  2:217353/149‑1 (MQ=255)
      gttgcgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGt                                                                                                                                                 <  8:303689/148‑1 (MQ=255)
      gttgcgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGt                                                                                                                                                 <  2:368234/148‑1 (MQ=255)
          cgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAg                                                                                                                                            <  2:43611/149‑1 (MQ=255)
           gttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           >  3:247670/1‑149 (MQ=255)
           gttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           <  7:181031/149‑1 (MQ=255)
            ttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCg                                                                                                                                          >  3:170273/1‑149 (MQ=255)
               cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCggg                                                                                                                                        <  7:309845/148‑1 (MQ=255)
                aGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGtt                                                                                                                                      >  6:34776/1‑149 (MQ=255)
                  gTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCa                                                                                                                                    >  6:299337/1‑149 (MQ=255)
                    aGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                                  >  6:327624/1‑149 (MQ=255)
                      cTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCcg                                                                                                                                <  5:3258/149‑1 (MQ=25)
                       tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                               >  1:378441/1‑149 (MQ=255)
                       tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                               >  8:329547/1‑149 (MQ=255)
                       tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                               >  5:120691/1‑149 (MQ=255)
                       nggCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCACCGTAGGAGCGGGTTCACCcgc                                                                                                                               >  2:109996/2‑149 (MQ=255)
                         gCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGGTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTGAACCcgctc                                                                                                                             <  3:295854/149‑3 (MQ=255)
                           gATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           >  5:1780/1‑133 (MQ=255)
                           gATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           <  6:1780/133‑1 (MQ=255)
                            aTCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                          <  6:54412/149‑1 (MQ=255)
                             tCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                         <  3:71068/149‑1 (MQ=255)
                                  tCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                    >  8:4/1‑149 (MQ=255)
                                       cAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAg                                                                                                               <  5:299337/149‑1 (MQ=255)
                                         gggCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGc                                                                                                             <  5:327624/149‑1 (MQ=255)
                                               tGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGGTCACCCGCGAAGCGGCCAGCACAGGCCCAGGc                                                                                                       >  4:322143/1‑149 (MQ=255)
                                                   gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGaa                                                                                                   >  4:351573/1‑149 (MQ=255)
                                                   gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGa                                                                                                    <  6:92653/148‑1 (MQ=255)
                                                    tcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCGGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAat                                                                                                  <  1:109996/149‑1 (MQ=255)
                                                      ttCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                         <  4:237545/108‑1 (MQ=255)
                                                      ttCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                         >  3:237545/1‑108 (MQ=255)
                                                          aCCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGaa                                                                                            <  6:120691/149‑1 (MQ=255)
                                                            cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                                          <  8:158411/149‑1 (MQ=255)
                                                             aGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                       >  6:139837/1‑119 (MQ=255)
                                                               gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                       >  6:279372/1‑149 (MQ=255)
                                                               gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                       >  3:14085/1‑149 (MQ=255)
                                                                aaaaTGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAc                                                                                      <  7:257642/149‑1 (MQ=255)
                                                                    tGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACt                                                                                  <  7:329547/149‑1 (MQ=255)
                                                                                gTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAt                                                                      >  4:82415/1‑149 (MQ=255)
                                                                                 tGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           <  7:28520/79‑1 (MQ=255)
                                                                                 tGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                                           >  8:28520/1‑79 (MQ=255)
                                                                                    ttGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAttctt                                                                  <  5:34776/149‑1 (MQ=255)
                                                                                      gggCGTCTGGGTTTAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAttctt                                                                  <  7:165714/147‑1 (MQ=255)
                                                                                               ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                        <  7:246178/148‑1 (MQ=255)
                                                                                               ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                        >  6:38989/1‑148 (MQ=255)
                                                                                                gTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTca                                                      >  2:174805/1‑149 (MQ=255)
                                                                                                      gTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCACCCACTCGACAGGATCATTCTTCACACTTTTTCACACTCg                                                >  1:57757/1‑149 (MQ=255)
                                                                                                       tAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGc                                               >  2:296188/1‑149 (MQ=255)
                                                                                                            cTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGTAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAGCACTCGCCAGGATCATTGTTCACACTTTTTCACACTCTCCGTa                                           >  5:186927/1‑148 (MQ=38)
                                                                                                            cTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTAc                                          <  4:247670/149‑1 (MQ=255)
                                                                                                                  cGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCCCACTTTTTCACACTCGCCCTACAGGTAt                                    <  1:214481/149‑1 (MQ=255)
                                                                                                                   gAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                   <  6:25602/149‑1 (MQ=255)
                                                                                                                     gcgcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTaa                                 <  2:201104/149‑1 (MQ=255)
                                                                                                                        cTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGCAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACCCTTTTTCACACTCGCCCTACAGGTATTAACAg                              >  3:2950/1‑149 (MQ=255)
                                                                                                                           cgcgGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgc                             <  5:38989/147‑1 (MQ=255)
                                                                                                                              gggTAAACCCGTTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATACAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAg                                        <  3:105468/133‑1 (MQ=255)
                                                                                                                               ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                   <  4:226516/137‑1 (MQ=25)
                                                                                                                               ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAtt                                   >  3:226516/1‑137 (MQ=25)
                                                                                                                               ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGCTCGCCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa                               <  5:221871/141‑1 (MQ=18)
                                                                                                                               ggTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGCGCTCGCCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa                               >  6:221871/1‑141 (MQ=18)
                                                                                                                                 tAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATACTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAgg                     >  2:77688/1‑149 (MQ=255)
                                                                                                                                    aCCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACCCTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg                  >  5:151266/1‑149 (MQ=255)
                                                                                                                                     cccGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGc                 <  3:351573/149‑1 (MQ=255)
                                                                                                                                        gctgctGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTc              <  1:296188/149‑1 (MQ=255)
                                                                                                                                         ctgctgCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCt             <  2:378441/149‑1 (MQ=255)
                                                                                                                                           gctgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                              <  6:42516/114‑1 (MQ=255)
                                                                                                                                           gctgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                              >  5:42516/1‑114 (MQ=255)
                                                                                                                                              gcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt         >  6:136237/1‑148 (MQ=255)
                                                                                                                                                  gTCCCCGTAGCAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATATAACACAACCACTCGACAGCATCATTCTTCACACTTTTTAACACACGCCCTACAGGTATTAACAGCCGCCAAGGCAGCGTCTGCCATAGCa    >  4:311111/1‑149 (MQ=38)
                                                                                                                                                    ccccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCata  <  4:170273/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATA  >  NC_002947/3914709‑3915005

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