Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,812,230 A→G 6.5% intergenic (+105/‑288) PP_4238 → / → PP_5665 hypothetical protein/McbA family microcin B17‑like processing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,812,2300AG 6.5% 227.1 / 4.4 77intergenic (+105/‑288)PP_4238/PP_5665hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base G (3/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATT  >  NC_002947/4812082‑4812372
                                                                                                                                                    |                                                                                                                                              
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                               >  2:322388/1‑150 (MQ=255)
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                               >  1:2003749/1‑150 (MQ=255)
 cacaTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAAc                                                                                                                                              <  1:322388/150‑1 (MQ=255)
  acaTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACa                                                                                                                                             >  1:959832/1‑150 (MQ=255)
    aTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATc                                                                                                                                           >  1:1125623/1‑150 (MQ=255)
        gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAGTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAACTATCTAGAAACAATTCTTCGGAGCATGGTGATGTTTGACATCAGCa                                                                                                                                       <  1:2186887/150‑1 (MQ=255)
          gATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACACAGATATAACACTGCCGGCGAAAATATCTAGAAACAAGACTTCGGAGCAGGCTGAGGTCGAACATCAGCATg                                                                                                                                     >  2:1092521/1‑150 (MQ=255)
           aTTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAAATCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGc                                                                                                                                    >  1:2124090/1‑150 (MQ=255)
           aTTGGACGCAGGCCGTTCTTTCGCAAACCGTAAAAGTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGc                                                                                                                                    <  2:843342/150‑1 (MQ=255)
               gACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATg                                                                                                                                     <  1:161901/145‑1 (MQ=255)
               gACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATg                                                                                                                                     >  2:161901/1‑145 (MQ=255)
                 cGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                              <  2:1637130/150‑1 (MQ=255)
                    aGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGcccc                                                                                                                           <  1:1377537/150‑1 (MQ=255)
                       ccGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGATTGACATCAGCAt                                                                                                                                      >  1:1925775/1‑136 (MQ=255)
                        cGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACCTCAGCATGCTCGCGCCCCttat                                                                                                                       <  2:574769/150‑1 (MQ=255)
                         gTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatt                                                                                                                      >  2:401015/1‑150 (MQ=255)
                          gtCTGTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatta                                                                                                                     >  1:81701/2‑150 (MQ=255)
                             tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                  <  2:1125623/150‑1 (MQ=255)
                             tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                  <  1:488592/150‑1 (MQ=255)
                             tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCCTATTaaaa                                                                                                                  <  1:2045480/150‑1 (MQ=255)
                              ttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCGGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAAc                                                                                                                 >  2:1978306/1‑150 (MQ=255)
                                cGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAAATCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACtt                                                                                                               <  2:2124090/150‑1 (MQ=255)
                                  cAAGCAGTAAATCTCCCGAAATCAATCTCACGTCCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                             <  2:609664/150‑1 (MQ=255)
                                  cAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCGGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                             <  1:1978306/150‑1 (MQ=255)
                                     gCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCacac                                                                                                          >  2:2097747/1‑150 (MQ=255)
                                        gTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACcgc                                                                                                       <  1:926395/150‑1 (MQ=255)
                                         tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACGCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                             >  2:1180569/1‑143 (MQ=255)
                                         tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACGCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                             <  1:1180569/143‑1 (MQ=255)
                                         tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACcgcg                                                                                                      >  2:1679982/1‑150 (MQ=255)
                                          aaaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCgcgc                                                                                                     <  2:959832/150‑1 (MQ=255)
                                                cccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                               <  2:1646307/150‑1 (MQ=255)
                                                 ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                              >  2:741312/1‑150 (MQ=255)
                                                 ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                              >  1:1315576/1‑150 (MQ=255)
                                                 ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATACTCGCACCCCGTATTAAAACTTCACACCGCGCATAACaa                                                                                              >  2:1938626/1‑150 (MQ=255)
                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                             <  1:1991323/150‑1 (MQ=255)
                                                       tcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAaca                                                                                        >  1:89860/1‑150 (MQ=255)
                                                       tcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCACCTATTAAAACTTCACACCGCGCATAACAAAGAaca                                                                                        <  1:1938626/150‑1 (MQ=255)
                                                        caatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAacac                                                                                       <  1:2097747/150‑1 (MQ=255)
                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                  >  2:1751896/1‑150 (MQ=255)
                                                                  aCCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                       >  2:1456471/1‑92 (MQ=255)
                                                                  aCCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                       <  1:1456471/92‑1 (MQ=255)
                                                                    ccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                           <  1:1092521/150‑1 (MQ=255)
                                                                     cGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                  <  1:749062/142‑1 (MQ=255)
                                                                     cGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                  >  2:749062/1‑142 (MQ=255)
                                                                          aaTCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACGCTCACACGCTCACCACaga                                                                     <  2:1315576/150‑1 (MQ=255)
                                                                              tgtCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCtt                                                                 >  1:2089549/3‑150 (MQ=255)
                                                                                   ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCTGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCATTATTAAAAATCCACACCGCGCACAACAAAAAACACGCACACGCTCACCACAGAGCTTCCGCa                                                            >  2:2259810/1‑150 (MQ=255)
                                                                                           gTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGaa                                                    <  1:741312/150‑1 (MQ=255)
                                                                                              tAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCa                                                 <  1:401015/150‑1 (MQ=255)
                                                                                                aCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGa                                               >  1:1130671/1‑150 (MQ=255)
                                                                                                  ccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAg                                             >  1:87729/1‑150 (MQ=255)
                                                                                                   cTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTCTTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGc                                            >  1:2247307/1‑150 (MQ=255)
                                                                                                       aaCTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCACGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                          <  1:1662587/100‑1 (MQ=255)
                                                                                                       aaCTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCACGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                          >  2:1662587/1‑100 (MQ=255)
                                                                                                          tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                       >  2:741604/1‑52 (MQ=255)
                                                                                                          tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                       <  1:741604/52‑1 (MQ=255)
                                                                                                          tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCCCAACTCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTAc                                     >  1:1210526/1‑150 (MQ=255)
                                                                                                            aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACCCAGCGCGCATACAAGCACAATCACACGCTCACTACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACAt                                   >  1:862142/1‑150 (MQ=255)
                                                                                                            aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACAt                                   <  2:377764/150‑1 (MQ=255)
                                                                                                                  tCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAatat                             <  2:2003749/150‑1 (MQ=255)
                                                                                                                    tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTCCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                           <  1:2259810/150‑1 (MQ=255)
                                                                                                                    caGAAACAATACTTCATACCATACTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                           <  2:2247307/149‑1 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                  <  2:318776/60‑1 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                  >  1:318776/1‑60 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaa                        <  1:1751896/150‑1 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCAAGAATCAGAAGCTTCTTACATAAATATaagaa                        <  1:383062/150‑1 (MQ=255)
                                                                                                                        aaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaag                       <  1:1679982/150‑1 (MQ=255)
                                                                                                                        aaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaag                       >  2:1936596/1‑150 (MQ=255)
                                                                                                                            aTTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACt                   >  2:1202520/1‑150 (MQ=255)
                                                                                                                                 tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACACAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa              <  1:691924/150‑1 (MQ=255)
                                                                                                                                  cGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAAc             <  1:1202520/150‑1 (MQ=255)
                                                                                                                                   ggAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCAAAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACt            <  2:81701/150‑1 (MQ=255)
                                                                                                                                      gCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaa         >  1:1190850/1‑150 (MQ=255)
                                                                                                                                       cATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAAtata                            >  2:433921/1‑130 (MQ=255)
                                                                                                                                       cATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAAtata                            <  1:433921/130‑1 (MQ=255)
                                                                                                                                        aTGCTGATGTTTGACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaaa       >  1:340312/1‑150 (MQ=255)
                                                                                                                                         tGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAc      <  2:1130671/150‑1 (MQ=255)
                                                                                                                                         tGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAc      <  2:1190850/150‑1 (MQ=255)
                                                                                                                                             gATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAtt  >  2:2270990/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATT  >  NC_002947/4812082‑4812372

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