Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,817,178 A→T 5.2% intergenic (+104/‑289) 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,817,1780AT 5.2% 1028.8 / 30.7 307intergenic (+104/‑289)PP_4238/PP_5665hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base A (153/138);  new base T (9/7);  total (162/145)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 323 total aligned reads displayed.

ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAGA  >  NC_002947/4817031‑4817327
                                                                                                                                                   |                                                                                                                                                     
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATGCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                     <  2:2762102/150‑1 (MQ=255)
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                     <  1:6050190/150‑1 (MQ=255)
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                     >  1:2850080/1‑150 (MQ=255)
acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                     <  2:4758452/150‑1 (MQ=255)
    aTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATc                                                                                                                                                 >  1:5037976/1‑150 (MQ=255)
     ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                                <  2:6210894/150‑1 (MQ=255)
     ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                                >  1:5641587/1‑150 (MQ=255)
       ggTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGc                                                                                                                                              >  1:2537726/1‑150 (MQ=255)
         tGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCAt                                                                                                                                            <  2:5699205/150‑1 (MQ=255)
            ttGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCt                                                                                                                                         >  1:3566165/1‑150 (MQ=255)
             tGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTc                                                                                                                                        >  1:4975965/1‑150 (MQ=255)
                 cGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                                    <  1:491807/150‑1 (MQ=255)
                 cGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                                    <  2:729995/150‑1 (MQ=255)
                  gCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGcc                                                                                                                                   >  1:375973/1‑150 (MQ=255)
                      gCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtt                                                                                                                               >  2:6613948/1‑150 (MQ=255)
                       ccGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCCAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtta                                                                                                                              >  1:5737820/1‑150 (MQ=255)
                         gTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatt                                                                                                                            >  1:6602717/1‑150 (MQ=255)
                           tCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaa                                                                                                                          >  2:5413541/1‑150 (MQ=255)
                             tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                        <  1:5279854/150‑1 (MQ=255)
                             tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                        <  1:4065660/150‑1 (MQ=255)
                                 gCAAGCAGTAAATCTCCCGAAATCGATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTc                                                                                                                    <  2:6295392/150‑1 (MQ=255)
                                  cAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTca                                                                                                                   <  1:3510667/150‑1 (MQ=255)
                                   aaGCAGTAAATCTCCCGAAATCAATCTCCCGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTcac                                                                                                                  <  1:5188531/150‑1 (MQ=255)
                                      cAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACAcc                                                                                                               <  1:89328/150‑1 (MQ=255)
                                           aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                          <  1:3312514/150‑1 (MQ=255)
                                           aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGAATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                          <  1:6334341/150‑1 (MQ=255)
                                               tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                      <  2:1198737/150‑1 (MQ=255)
                                               tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                      >  2:4460709/1‑150 (MQ=255)
                                               tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                      >  1:5497074/1‑150 (MQ=255)
                                               tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                      >  2:2549202/1‑150 (MQ=255)
                                                cccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                                     >  2:5945643/1‑150 (MQ=255)
                                                 ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                                    >  1:3039114/1‑150 (MQ=255)
                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                   <  1:6197867/150‑1 (MQ=255)
                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                   <  2:251003/150‑1 (MQ=255)
                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                   <  2:3614091/150‑1 (MQ=255)
                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                   <  1:5933157/150‑1 (MQ=255)
                                                  cGAAAACAATCTCACGACCCGGCAAATCTATCCCAACGCAGGTCTAACCATGAAACTCAAAATATCTATAAAAATTTCTTCGGAGCATGCTGATGTCAAAAATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                   <  1:2903880/150‑1 (MQ=255)
                                                     aatcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                                <  1:4790236/150‑1 (MQ=255)
                                                     aatcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                                <  2:2287036/150‑1 (MQ=255)
                                                     aatcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                                <  2:5025697/150‑1 (MQ=255)
                                                        caatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAacac                                                                                             <  2:375973/150‑1 (MQ=255)
                                                         aatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACt                                                                                            <  1:1475684/150‑1 (MQ=255)
                                                         aatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACt                                                                                            >  1:764163/1‑150 (MQ=255)
                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAAc                                                                                                                                                    <  2:882004/90‑1 (MQ=255)
                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                        >  1:4866956/1‑150 (MQ=255)
                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                        <  2:3328243/150‑1 (MQ=255)
                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGCCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                        >  2:5286323/1‑150 (MQ=255)
                                                                    ccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                                 >  2:384709/1‑150 (MQ=255)
                                                                     cGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacc                                                                                >  2:2925340/1‑150 (MQ=255)
                                                                       gCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcaccac                                                                              >  2:4540410/1‑150 (MQ=255)
                                                                         aaaTCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACag                                                                            <  1:1901369/150‑1 (MQ=255)
                                                                           aTCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACagag                                                                          <  2:5826878/150‑1 (MQ=255)
                                                                            tCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGc                                                                         >  2:1776879/1‑150 (MQ=255)
                                                                            tCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGc                                                                         >  2:3679543/1‑150 (MQ=255)
                                                                              tATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCtt                                                                       >  1:2455404/1‑150 (MQ=255)
                                                                                tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                     <  1:1956659/150‑1 (MQ=255)
                                                                                 ccccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCg                                                                    <  2:5597519/150‑1 (MQ=255)
                                                                                 ccccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCg                                                                    >  1:4230811/1‑150 (MQ=255)
                                                                                   ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                                <  1:2242405/72‑1 (MQ=255)
                                                                                   ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCa                                                                  >  2:5647449/1‑150 (MQ=255)
                                                                                   ccACGAAGTTCTAACCCTGCAAATCAAATTATCTAGATACAATTCTTCGTAGCATGCTGTTGTCAAACATCAGCATGCTCTCGACCCTTATTAATACTTAACACCGCGCTTATCAATGAACACTCACACGCTCACCACAGAGCTTTCGCa                                                                  <  1:4432326/150‑1 (MQ=255)
                                                                                           gTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGaa                                                          >  1:2693865/1‑150 (MQ=255)
                                                                                            tCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCCCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAAt                                                         <  2:6514278/150‑1 (MQ=255)
                                                                                                 cccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGaa                                                    <  2:350773/150‑1 (MQ=255)
                                                                                                 cccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGaa                                                    <  2:1052841/150‑1 (MQ=255)
                                                                                                    tGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGct                                                 >  2:6388783/1‑150 (MQ=255)
                                                                                                    tGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGct                                                 >  2:1078814/1‑150 (MQ=255)
                                                                                                     gCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGctt                                                <  2:1513263/150‑1 (MQ=255)
                                                                                                      cAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc                                               <  1:1170102/150‑1 (MQ=255)
                                                                                                           cAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACa                                          >  2:2409738/1‑150 (MQ=255)
                                                                                                            aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACAt                                         >  2:287825/1‑150 (MQ=255)
                                                                                                             aaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATa                                        >  2:6737061/1‑150 (MQ=255)
                                                                                                                    tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                                 <  2:4186887/150‑1 (MQ=255)
                                                                                                                      gAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaaga                               <  1:2549202/150‑1 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaa                              >  1:1569005/1‑150 (MQ=255)
                                                                                                                       aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaa                              >  2:2172093/1‑150 (MQ=255)
                                                                                                                        aaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaag                             <  1:2710569/150‑1 (MQ=255)
                                                                                                                          cAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                                 <  2:1906698/144‑1 (MQ=255)
                                                                                                                             ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTg                        >  1:1334845/1‑150 (MQ=255)
                                                                                                                               cttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaa                      <  2:2673659/150‑1 (MQ=255)
                                                                                                                                ttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaa                     >  1:4498870/1‑150 (MQ=255)
                                                                                                                                ttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaa                     >  1:4839071/1‑150 (MQ=255)
                                                                                                                                 tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                    <  1:211161/150‑1 (MQ=255)
                                                                                                                                 tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                    <  1:3890018/150‑1 (MQ=255)
                                                                                                                                   ggAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACt                  >  2:3211107/1‑150 (MQ=255)
                                                                                                                                    gAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTa                 >  1:4142226/1‑150 (MQ=255)
                                                                                                                                        aTGCTGATGTTTGACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaata             >  1:5299612/1‑148 (MQ=255)
                                                                                                                                        aTGCTGATGTTTGACATCAGCATGCTCGCGCCCCTAATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaaa             >  2:5204459/1‑150 (MQ=255)
                                                                                                                                          gCTGATGTTTGACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAcc           >  1:102367/1‑150 (MQ=255)
                                                                                                                                          gCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAg                                                      <  2:4116469/107‑1 (MQ=255)
                                                                                                                                            tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGACCACTAAAAACCAt         >  2:6126389/1‑150 (MQ=255)
                                                                                                                                            tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt         >  1:6822081/1‑150 (MQ=255)
                                                                                                                                            tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt         <  2:1384275/150‑1 (MQ=255)
                                                                                                                                            tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt         >  2:836177/1‑150 (MQ=255)
                                                                                                                                             gATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTa                                            <  2:4263886/114‑1 (MQ=255)
                                                                                                                                                gTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATc     <  2:2850080/150‑1 (MQ=255)
                                                                                                                                                 tCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCa    >  2:3310972/1‑150 (MQ=255)
                                                                                                                                                   aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAGa  >  1:5830087/1‑150 (MQ=255)
                                                                                                                                                   aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAGa  >  1:994738/1‑150 (MQ=255)
                                                                                                                                                   aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAGa  >  2:689198/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAGA  >  NC_002947/4817031‑4817327

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