Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,812,228 C→T 7.2% intergenic (+103/‑290) 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,2280CT 7.2% 240.9 / 7.1 83intergenic (+103/‑290)PP_4238/PP_5665hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base C (37/40);  new base T (3/3);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GAGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCA  >  NC_002947/4812079‑4812376
                                                                                                                                                      |                                                                                                                                                    
gagaCACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTc                                                                                                                                                      >  2:1755015/1‑150 (MQ=255)
gagaCACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTc                                                                                                                                                      >  1:1394252/1‑150 (MQ=255)
  gaCACATTAGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaa                                                                                                                                                    >  2:1420729/1‑150 (MQ=255)
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTACCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGTAGCATGCTGATGTCaaa                                                                                                                                                   <  2:1233674/150‑1 (MQ=255)
        ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTACATCTTACGCCATCAAACCCACGTAGCGGGCCACTTAACCCCACGCATGTGTAACCCTGCAACGCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                              <  2:1477884/150‑1 (MQ=255)
        ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                              <  2:1410879/150‑1 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                              >  2:1442138/1‑147 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                              <  1:1442138/147‑1 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                           <  2:2008973/150‑1 (MQ=255)
             gATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAACTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                           >  2:1617523/1‑148 (MQ=255)
             gATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCa                                                                                                                                           <  1:1617523/148‑1 (MQ=255)
                tGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATGCTc                                                                                                                                      <  1:1512252/150‑1 (MQ=255)
                   aCGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTcgcg                                                                                                                                   <  2:2182616/150‑1 (MQ=255)
                     gCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAACCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTGACCGACCGACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCGGATGTCAAACATCAGCATGCTCGCGcc                                                                                                                                 >  1:1022441/1‑150 (MQ=255)
                       aGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGcccc                                                                                                                               <  2:1839931/150‑1 (MQ=255)
                        ggCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCt                                                                                                                              <  1:2241156/150‑1 (MQ=255)
                         gCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtt                                                                                                                             <  1:683513/150‑1 (MQ=255)
                           cGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTc                                                                                                                                                      >  2:1205002/1‑123 (MQ=255)
                           cGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTc                                                                                                                                                      <  1:1205002/123‑1 (MQ=255)
                           cGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttat                                                                                                                           >  2:370357/1‑150 (MQ=255)
                               cTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGc                                                                                                                                            <  2:484715/129‑1 (MQ=255)
                               cTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGc                                                                                                                                            >  1:484715/1‑129 (MQ=255)
                               cTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaa                                                                                                                       <  2:1394252/150‑1 (MQ=255)
                                 ttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAAc                                                                                                                     <  1:1755015/150‑1 (MQ=255)
                                    gCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTc                                                                                                                  >  2:99409/1‑150 (MQ=255)
                                            tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCCGCATGCTCGCGCCCGTAAATAAAAGATCACACcgcg                                                                                                          >  2:1719953/1‑150 (MQ=255)
                                                   cccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAcg                                                                                                   <  2:286686/150‑2 (MQ=255)
                                                   cccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                                   >  2:409347/1‑150 (MQ=255)
                                                     cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                 <  2:1853967/150‑1 (MQ=255)
                                                            aatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGGGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACt                                                                                          >  2:1336346/1‑150 (MQ=255)
                                                             atcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGACCACTc                                                                                         <  2:953327/150‑1 (MQ=255)
                                                              tgtcACTACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTca                                                                                        >  2:648084/3‑150 (MQ=255)
                                                              tctcACGACCCGGCAAATCTATCCCCACGCATGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTAATGTCAAACATCAGCATGATCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAacaatga                                                                                        >  1:443681/1‑146 (MQ=255)
                                                              tctcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGc                                                                                                                                        <  2:688449/102‑1 (MQ=255)
                                                              tctcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGc                                                                                                                                        >  1:688449/1‑102 (MQ=255)
                                                                    gACCCGGCAAATCTATCCCCACGCAGGTCTAACACTGCAACACAAAATCTCTAGAAACAATTCTTCGGAGCATGCTGATATCAAACCTCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCt                                                                                  >  2:275565/1‑150 (MQ=255)
                                                                       cccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                      >  2:1045389/1‑112 (MQ=255)
                                                                       cccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                      <  1:1045389/112‑1 (MQ=255)
                                                                           gCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCACAACAAAGAACACTCACACGCTcaccac                                                                            >  1:568401/1‑150 (MQ=255)
                                                                                tCTATCACCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGc                                                                       >  2:150301/1‑150 (MQ=255)
                                                                                 cTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCt                                                                      >  1:2086512/1‑150 (MQ=255)
                                                                                  tATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCtt                                                                     <  1:1719953/150‑1 (MQ=255)
                                                                                    tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGGGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                   <  1:1336346/150‑1 (MQ=255)
                                                                                    tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                               <  1:1782568/138‑1 (MQ=255)
                                                                                    tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                               >  2:1782568/1‑138 (MQ=255)
                                                                                    tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                   <  1:370357/150‑1 (MQ=255)
                                                                                      cccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                                  <  1:2126557/85‑1 (MQ=255)
                                                                                      cccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                                  >  2:2126557/1‑85 (MQ=255)
                                                                                      cccACGCAGCAATCACCCTGCAACGCAAAATATCTAGAAACAATTCTTCGTAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGc                                                                 >  2:101712/1‑150 (MQ=255)
                                                                                       ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATg                                                                                                                                         <  1:1032295/77‑1 (MQ=37)
                                                                                       ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATg                                                                                                                                         >  2:1032295/1‑77 (MQ=37)
                                                                                       ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCa                                                                >  1:795348/1‑150 (MQ=255)
                                                                                       ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCa                                                                >  1:2182415/1‑150 (MQ=255)
                                                                                         aCGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGccc                                                                                                                                <  2:155830/84‑1 (MQ=255)
                                                                                         aCGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGccc                                                                                                                                >  1:155830/1‑84 (MQ=255)
                                                                                                tCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAacac                                                                                           <  2:727366/114‑1 (MQ=255)
                                                                                                tCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAacac                                                                                           >  1:727366/1‑114 (MQ=255)
                                                                                                    aCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATGCTc                                                                                                                                      >  2:137586/1‑67 (MQ=25)
                                                                                                    aCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCATGCTc                                                                                                                                      <  1:137586/67‑1 (MQ=25)
                                                                                                    aCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGa                                                   <  1:231886/150‑1 (MQ=255)
                                                                                                      ccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAg                                                 >  2:966199/1‑150 (MQ=255)
                                                                                                          cAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc                                             <  1:1420729/150‑1 (MQ=255)
                                                                                                             cTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTa                                          <  1:528565/150‑1 (MQ=255)
                                                                                                              tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTAc                                         <  1:409347/150‑1 (MQ=255)
                                                                                                              tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCACGCATAACAAAGaa                                                                                              <  2:1058921/97‑1 (MQ=255)
                                                                                                              tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCACGCATAACAAAGaa                                                                                              >  1:1058921/1‑97 (MQ=255)
                                                                                                                     atCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAata                                  <  2:2182415/150‑1 (MQ=255)
                                                                                                                      tCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGctt                                              >  2:1492147/1‑137 (MQ=255)
                                                                                                                      tCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGctt                                              <  1:1492147/137‑1 (MQ=255)
                                                                                                                       cTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                               <  2:1821258/103‑1 (MQ=255)
                                                                                                                       cTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                               >  1:1821258/1‑103 (MQ=255)
                                                                                                                        tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                               >  1:579147/1‑150 (MQ=255)
                                                                                                                          gAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaaga                             >  2:1860239/1‑150 (MQ=255)
                                                                                                                           aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                                 <  1:1829423/65‑1 (MQ=255)
                                                                                                                           aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                                 >  2:1829423/1‑65 (MQ=255)
                                                                                                                            aaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaag                           <  1:99409/150‑1 (MQ=255)
                                                                                                                                 ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATaaa                                    <  1:1492760/136‑1 (MQ=255)
                                                                                                                                 ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATaaa                                    >  2:1492760/1‑136 (MQ=255)
                                                                                                                                    ttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaa                   >  2:1559781/1‑150 (MQ=255)
                                                                                                                                     tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGAGCCTTATACACATCACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                  <  1:150301/150‑1 (MQ=255)
                                                                                                                                      cGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAAc                 <  2:579147/150‑1 (MQ=255)
                                                                                                                                           cATGCTGATGTTTGACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaa            >  1:33725/1‑150 (MQ=255)
                                                                                                                                             tGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACACACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCCTACAAAAATACAAGACGAACTCAAAACAAAAAcc         >  2:992067/1‑150 (MQ=255)
                                                                                                                                                tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt       >  1:2379167/1‑150 (MQ=255)
                                                                                                                                                 gATCTCAAACATCAACCTGCTCGCGCCCTATATGAAAACTTCACACCGCTCATAACAAAGCACACTCACACGCTCACCACACAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAATAAGCACTGAAACCTAAAAACCAtt      >  2:31636/1‑150 (MQ=255)
                                                                                                                                                     tCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCa  <  2:2379167/150‑1 (MQ=255)
                                                                                                                                                     tCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCACAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCa  <  2:568401/150‑1 (MQ=255)
                                                                                                                                                      cAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                  >  2:2291458/1‑133 (MQ=255)
                                                                                                                                                      cAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                  <  1:2291458/133‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCA  >  NC_002947/4812079‑4812376

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