Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,894 A→G 6.2% intergenic (‑19/+167) PP_3629 ← / ← opdT‑II hypothetical protein/tyrosine‑specific outer membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,8940AG 6.2% 239.0 / 2.8 79intergenic (‑19/+167)PP_3629/opdT‑IIhypothetical protein/tyrosine‑specific outer membrane porin D
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base G (3/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGG  >  NC_002947/4125746‑4126036
                                                                                                                                                    |                                                                                                                                              
cTGAAGGGTGAGTGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGaa                                                                                                                                               <  2:436453/150‑1 (MQ=255)
cTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGCAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGaa                                                                                                                                               >  1:1366051/1‑150 (MQ=255)
cTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGaa                                                                                                                                               <  1:939743/150‑1 (MQ=255)
cTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGaa                                                                                                                                               >  1:1791834/1‑150 (MQ=255)
cTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGaa                                                                                                                                               >  2:363362/1‑150 (MQ=255)
          aTCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCtt                                                                                                                                     <  1:2287751/150‑1 (MQ=255)
          aTCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCtt                                                                                                                                     >  2:1294680/1‑150 (MQ=255)
            cggcggGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGa                                                                                                                                   >  1:2089205/1‑150 (MQ=255)
                   aGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGc                                                                                                                            <  2:2089205/150‑1 (MQ=255)
                         gTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATa                                                                                                                      <  1:924502/150‑1 (MQ=255)
                         gTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATa                                                                                                                      <  2:2337963/150‑1 (MQ=255)
                           gCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGGAGGTGAGCCTTGAGTCTAGCCCCATAGc                                                                                                                    >  2:892702/1‑150 (MQ=255)
                           gCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGGAGGTGAGCCTTGAGTCTAGCCCCATAGc                                                                                                                    >  2:892712/1‑150 (MQ=255)
                           gCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGc                                                                                                                    >  2:892741/1‑150 (MQ=255)
                             acaGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAg                                                                                                                                         >  1:1020153/1‑127 (MQ=255)
                             acaGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAg                                                                                                                                         <  2:2277449/127‑1 (MQ=255)
                             acaGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAg                                                                                                                                         <  2:1020153/127‑1 (MQ=255)
                             acaGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAg                                                                                                                                         >  1:2277449/1‑127 (MQ=255)
                              caGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGcccc                                                                                                                 <  1:1674150/150‑1 (MQ=255)
                                gCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCgg                                                                                                               >  1:2209467/1‑150 (MQ=255)
                                       cccGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTa                                                                                                                              <  1:896821/128‑1 (MQ=255)
                                       cccGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTa                                                                                                                              >  2:896821/1‑128 (MQ=255)
                                       cccGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATa                                                                                                                      <  2:2309692/136‑1 (MQ=255)
                                       cccGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATa                                                                                                                      >  1:2309692/1‑136 (MQ=255)
                                            cTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGc                                                                                                   <  1:849204/150‑1 (MQ=255)
                                              aTCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGt                                                                                                                                 >  2:1577003/1‑118 (MQ=255)
                                              aTCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGt                                                                                                                                 <  1:1577003/118‑1 (MQ=255)
                                                     aaaGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGacca                                                                                          >  1:2248762/1‑150 (MQ=255)
                                                         ttCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTg                                                                                      >  1:1151459/1‑150 (MQ=255)
                                                            gggCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTa                                                                                                     <  1:2073310/132‑1 (MQ=255)
                                                            gggCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTa                                                                                                     >  2:2073310/1‑132 (MQ=255)
                                                                   gAAGTTCGTCGCCTCTTCGGCATCTTTGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGc                                                                            <  1:2335004/150‑1 (MQ=255)
                                                                     aGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAg                                                                          <  1:676359/150‑1 (MQ=255)
                                                                      gTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAg                                                                                                                                              <  1:1657034/81‑1 (MQ=255)
                                                                      gTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAg                                                                                                                                              >  2:1657034/1‑81 (MQ=255)
                                                                              ccTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAg                                                                 >  1:1261620/1‑150 (MQ=255)
                                                                                tctTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATTCAAATGCAGGATTGATAGGc                                                               >  1:546945/1‑150 (MQ=255)
                                                                                tctTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGc                                                               <  1:2426651/150‑1 (MQ=255)
                                                                                  ttCGGCATCGTGGTCGTACCACAGGCAGATGGTAATCATGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGa                                                             >  1:1145466/1‑150 (MQ=255)
                                                                                     ggCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAAt                                                                              <  1:1671899/130‑1 (MQ=255)
                                                                                     ggCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAAt                                                                              >  2:1671899/1‑130 (MQ=255)
                                                                                      gCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGAt                                                                                    >  1:754669/1‑123 (MQ=255)
                                                                                      gCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGAt                                                                                    <  2:754669/123‑1 (MQ=255)
                                                                                                     cacaGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAg                                                                          >  1:136933/1‑118 (MQ=255)
                                                                                                     cacaGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAg                                                                          <  2:136933/118‑1 (MQ=255)
                                                                                                       caGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAg                                                        >  1:132207/1‑134 (MQ=255)
                                                                                                       caGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAg                                                        <  2:132207/134‑1 (MQ=255)
                                                                                                        aGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGGAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGcc                                       >  1:421091/1‑150 (MQ=255)
                                                                                                        aGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGcc                                       >  1:421204/1‑150 (MQ=255)
                                                                                                            aGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCaa                                   >  1:1738605/1‑150 (MQ=255)
                                                                                                             gATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGGAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAg                                  <  1:892702/150‑1 (MQ=255)
                                                                                                             gATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGGAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAg                                  <  1:892712/150‑1 (MQ=255)
                                                                                                             gATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGa                                                             >  2:1171277/1‑123 (MQ=255)
                                                                                                             gATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGa                                                             <  1:1171277/123‑1 (MQ=255)
                                                                                                             gATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAg                                  <  1:892741/150‑1 (MQ=255)
                                                                                                                ggTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCgg                               <  2:2248762/150‑1 (MQ=255)
                                                                                                                   gtTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAgc                            <  1:1294680/150‑1 (MQ=255)
                                                                                                                       ttGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCCAGCGTCGAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCa                        >  1:553714/1‑150 (MQ=255)
                                                                                                                          gCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCgg                                                                                                       <  1:1800209/68‑1 (MQ=255)
                                                                                                                          gCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCgg                                                                                                       >  2:1800209/1‑68 (MQ=255)
                                                                                                                              tCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATc                 <  1:401382/150‑1 (MQ=255)
                                                                                                                                   gTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAg                                                                 <  2:2106372/97‑1 (MQ=255)
                                                                                                                                   gTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAg                                                                 >  1:2106372/1‑97 (MQ=255)
                                                                                                                                         cTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCa      >  2:53533/1‑150 (MQ=255)
                                                                                                                                          tCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCaa     >  2:869644/1‑150 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGAt                                                                   <  1:1434968/87‑1 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGAt                                                                   >  2:1434968/1‑87 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCaaa    <  2:1261620/150‑1 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCaaa    >  1:2245657/1‑150 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCaaa    <  1:363362/150‑1 (MQ=255)
                                                                                                                                           ccAAGATCGAAGGTGAGCCTTGAGTCTACCCCCATAGCCCTGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCaaa    <  1:217554/150‑1 (MQ=255)
                                                                                                                                            cAAGATCGAAGGTGAGCCTTGATTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGa                                                                                     >  1:1856444/1‑68 (MQ=255)
                                                                                                                                            cAAGATCGAAGGTGAGCCTTGATTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGa                                                                                     <  2:1856444/68‑1 (MQ=255)
                                                                                                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCTGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGGTCGCCAGAAAATCCACCTAAACCAAAg   >  1:1567484/1‑150 (MQ=255)
                                                                                                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGAt                                                                       >  2:1164092/1‑82 (MQ=255)
                                                                                                                                            cAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGAt                                                                       <  1:1164092/82‑1 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCAc              <  1:2335141/138‑1 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAgg  >  1:1409366/1‑150 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAgg  >  1:870869/1‑150 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAgg  >  1:1745187/1‑150 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAgg  <  2:870869/150‑1 (MQ=255)
                                                                                                                                             aaGATCGAAGGTGAGACTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCAc              >  2:2335141/1‑138 (MQ=255)
                                                                                                                                                 tCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAGCGCCCCAAGc                                 <  2:4008/114‑1 (MQ=255)
                                                                                                                                                 tCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAGCGCCCCAAGc                                 >  1:4008/1‑114 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTTCTTGGCAGTCATGGTCTTGCTCCAAGATCGAAGGTGAGCCTTGAGTCTAGCCCCATAGCCCCGGCGTACCGGTAGCACCTGACCACCTGATGCAAATGCAGGATTGATAGGCGATCAAGCCGATAAAAAAAGCGCCCCAAGCGGAGCGCCAGAAAATCCACCTAAACCAAAGG  >  NC_002947/4125746‑4126036

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