Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,914,862 T→C 7.8% intergenic (‑76/‑170) 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,8620TC 7.8% 179.9 / 5.0 64intergenic (‑76/‑170)PP_3454/PP_3455transcriptional regulator/multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/3);  total (31/33)
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.39e-01

CAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACAC  >  NC_002947/3914724‑3915010
                                                                                                                                           |                                                                                                                                                    
cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                                             <  7:163846/148‑1 (MQ=37)
cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCggg                                                                                                                                             <  1:33343/148‑1 (MQ=37)
cAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCggg                                                                                                                                             <  2:71530/148‑1 (MQ=255)
 aGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGtt                                                                                                                                           >  8:298378/1‑149 (MQ=255)
 aGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGt                                                                                                                                            >  3:42065/1‑148 (MQ=255)
  ggTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                              <  1:181096/145‑1 (MQ=255)
  ggTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                              >  2:181096/1‑145 (MQ=255)
     aGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                                       >  7:132979/1‑149 (MQ=255)
       cTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAccc                                                                                                                                      <  1:73262/149‑1 (MQ=255)
                  cTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCa                                                                                                                          <  6:247713/149‑1 (MQ=255)
                            gcgcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                                                <  7:298378/149‑1 (MQ=255)
                              gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               <  8:146378/132‑1 (MQ=255)
                              gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               >  1:342848/1‑132 (MQ=255)
                              gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               >  7:146378/1‑132 (MQ=255)
                              gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               <  2:342848/132‑1 (MQ=255)
                              gcTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                                              >  8:83359/1‑149 (MQ=255)
                                    gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                                       >  1:201304/1‑118 (MQ=255)
                                    gtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCAcc                                                                                                                                       <  2:201304/118‑1 (MQ=255)
                                     tcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAat                                                                                                       <  4:42065/149‑1 (MQ=255)
                                      cTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               >  8:27909/1‑124 (MQ=255)
                                      cTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                               <  7:27909/124‑1 (MQ=255)
                                       ttCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAata                                                                                                      >  7:150386/1‑148 (MQ=255)
                                         cGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAatat                                                                                                     <  3:91773/147‑1 (MQ=255)
                                         cGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAata                                                                                                      <  8:150386/146‑1 (MQ=255)
                                           aCCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGaa                                                                                                 >  5:60432/1‑149 (MQ=255)
                                             cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                                               >  4:143026/1‑149 (MQ=255)
                                             cAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                                               >  1:217690/1‑149 (MQ=255)
                                              aGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAcc                                                                                               >  5:96247/1‑148 (MQ=255)
                                                gAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                            >  5:328371/1‑149 (MQ=255)
                                                 aaaaTGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAc                                                                                           <  4:48871/149‑1 (MQ=255)
                                                          ggCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACa                                                                                  <  3:328672/149‑1 (MQ=255)
                                                               aGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAAtata                                                                                                    >  5:118578/1‑126 (MQ=255)
                                                               aGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAAtata                                                                                                    <  6:118578/126‑1 (MQ=255)
                                                                 gTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAt                                                                           >  2:309340/1‑149 (MQ=255)
                                                                       gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAccc                                                                                              <  6:312541/124‑1 (MQ=255)
                                                                       gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAAccc                                                                                              >  5:312541/1‑124 (MQ=255)
                                                                       gggCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTca                                                                     <  8:132979/149‑1 (MQ=255)
                                                                         gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTc                                                                                      >  3:250908/1‑130 (MQ=255)
                                                                         gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTc                                                                                      <  4:250908/130‑1 (MQ=255)
                                                                              tGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtttt                                                              >  7:164249/1‑149 (MQ=255)
                                                                                ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                             >  5:53579/1‑148 (MQ=255)
                                                                                ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtttt                                                              <  6:328371/147‑1 (MQ=255)
                                                                                ggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTc                                                            <  8:164249/149‑1 (MQ=255)
                                                                                    gAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcaca                                                         <  1:309340/148‑1 (MQ=255)
                                                                                    gAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcaca                                                         <  6:53579/148‑1 (MQ=255)
                                                                                     aaGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtt                                                                <  8:47670/140‑1 (MQ=255)
                                                                                     aaGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtt                                                                >  7:47670/1‑140 (MQ=255)
                                                                                         aaGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                   >  7:16967/1‑149 (MQ=255)
                                                                                          aGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                   <  2:217690/148‑1 (MQ=255)
                                                                                                  gCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGt                                           <  3:143026/148‑1 (MQ=255)
                                                                                                    gAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTAt                                         >  6:161059/1‑148 (MQ=255)
                                                                                                       cgcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAAc                                     <  6:60432/149‑1 (MQ=255)
                                                                                                        gcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACa                                    <  6:96247/149‑1 (MQ=255)
                                                                                                            cgcgGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgccg                                >  6:17968/1‑149 (MQ=255)
                                                                                                                  gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                             <  6:68830/112‑1 (MQ=21)
                                                                                                                  gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                             >  5:68830/4‑115 (MQ=21)
                                                                                                                  gtgaaCCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAGTCTTCACACTTTTACACACTCGCCCTACAGGTATTAAAAGCCGCCCAg                           >  8:243953/4‑149 (MQ=25)
                                                                                                                      aCCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg                       >  1:172522/1‑149 (MQ=255)
                                                                                                                                 cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt              >  4:195366/1‑147 (MQ=255)
                                                                                                                                 cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCt              >  2:39097/1‑147 (MQ=255)
                                                                                                                                 cAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTa            <  5:161059/149‑1 (MQ=255)
                                                                                                                                    gTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCa         >  1:260594/1‑149 (MQ=255)
                                                                                                                                       cccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCata       <  2:260594/148‑1 (MQ=255)
                                                                                                                                       cccGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGAGCATTCTTCACACTTTTTCACACTCGCACTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCata       <  7:323121/148‑1 (MQ=255)
                                                                                                                                           tAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATacac  <  1:39097/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACAC  >  NC_002947/3914724‑3915010

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: