Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,817,177 C→T 5.1% 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,817,1770CT 5.1% 973.7 / 25.0 312intergenic (+103/‑290)PP_4238/PP_5665hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base C (151/143);  new base T (9/7);  total (161/151)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 316 total aligned reads displayed.

GAGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAG  >  NC_002947/4817028‑4817326
                                                                                                                                                     |                                                                                                                                                     
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATGCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                    <  2:2762102/150‑1 (MQ=255)
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                    <  1:6050190/150‑1 (MQ=255)
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                    <  2:1599061/150‑1 (MQ=255)
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                    >  1:2850080/1‑150 (MQ=255)
   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCGGATGTCaaa                                                                                                                                                    <  2:2549541/150‑1 (MQ=255)
      caTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACAt                                                                                                                                                 <  2:6836462/150‑1 (MQ=255)
        ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                               <  1:4004724/150‑1 (MQ=255)
          ggTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGc                                                                                                                                             >  1:2537726/1‑150 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCCCGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCGAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCa                                                                                                                                            <  1:2408421/150‑1 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCa                                                                                                                                            <  2:1061428/150‑1 (MQ=255)
           gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCa                                                                                                                                            <  1:4810293/150‑1 (MQ=255)
            tGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCAt                                                                                                                                           <  2:1781419/150‑1 (MQ=255)
               ttGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCt                                                                                                                                        >  1:3566165/1‑150 (MQ=255)
                    cGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCgcgc                                                                                                                                   <  2:729995/150‑1 (MQ=255)
                     gCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGcc                                                                                                                                  >  1:375973/1‑150 (MQ=255)
                        ggCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCt                                                                                                                               >  1:4720957/1‑150 (MQ=255)
                          ccGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCCAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtta                                                                                                                             >  1:5737820/1‑150 (MQ=255)
                           cGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttat                                                                                                                            >  1:5845872/1‑150 (MQ=255)
                            gTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatt                                                                                                                           >  1:6602717/1‑150 (MQ=255)
                             ttCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatta                                                                                                                          >  2:3033076/1‑150 (MQ=255)
                                tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                       <  2:4987854/150‑1 (MQ=255)
                                tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                       <  1:5280243/150‑1 (MQ=255)
                                tttCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaaaa                                                                                                                       >  2:5861524/1‑150 (MQ=255)
                                         cAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACAcc                                                                                                              <  1:89328/150‑1 (MQ=255)
                                              aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                         <  1:3312514/150‑1 (MQ=255)
                                              aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                         >  2:2474557/1‑150 (MQ=255)
                                              aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGAATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                         <  1:6334341/150‑1 (MQ=255)
                                              aaTCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTAAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa                                                                                                         <  2:6057859/150‑1 (MQ=255)
                                                 ctcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATaa                                                                                                      >  1:5025697/1‑150 (MQ=255)
                                                  tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                     <  2:635656/150‑1 (MQ=255)
                                                  tcCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAAc                                                                                                     >  2:4460709/1‑150 (MQ=255)
                                                    ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                                   >  2:1170102/1‑150 (MQ=255)
                                                    ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                                   >  1:3039114/1‑150 (MQ=255)
                                                    ccGAAATCAATCACACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCACAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCAGAACaa                                                                                                   >  2:5691020/1‑150 (MQ=255)
                                                     cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                  <  2:3614091/150‑1 (MQ=255)
                                                     cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                  <  1:5519318/150‑1 (MQ=255)
                                                     cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                  <  1:5933157/150‑1 (MQ=255)
                                                     cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                  <  1:6197867/150‑1 (MQ=255)
                                                     cGAAAACAATCTCACGACCCGGCAAATCTATCCCAACGCAGGTCTAACCATGAAACTCAAAATATCTATAAAAATTTCTTCGGAGCATGCTGATGTCAAAAATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                  <  1:2903880/150‑1 (MQ=255)
                                                        aatcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                               <  2:2287364/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                                                                                       <  2:3328243/150‑1 (MQ=255)
                                                                tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                       >  2:142314/1‑150 (MQ=255)
                                                                tcACGACCCGGCAAATCTATCCCCACGCAGGCCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcaca                                                                                       >  2:5286323/1‑150 (MQ=255)
                                                                 cACGACCCGGCAAATCTATCCCCACGCATGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCacac                                                                                      <  2:2067305/150‑1 (MQ=255)
                                                                  aCGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATAAAAACTTCACACCGCGCATAACAAAGAACACTCACACg                                                                                     >  2:2538486/1‑150 (MQ=255)
                                                                    gACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCt                                                                                   >  2:5645360/1‑150 (MQ=255)
                                                                      cccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTca                                                                                 <  1:510864/150‑1 (MQ=255)
                                                                       ccGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcac                                                                                >  2:384709/1‑150 (MQ=255)
                                                                        cGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacc                                                                               >  2:4159114/1‑150 (MQ=255)
                                                                         ggCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacca                                                                              <  1:1812272/150‑1 (MQ=255)
                                                                         ggCAAATCTATACCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGACACAATTCTTCGGTGCATGCTGATGTCATACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacca                                                                              <  1:2123907/150‑1 (MQ=255)
                                                                          gCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcaccac                                                                             >  2:4540410/1‑150 (MQ=255)
                                                                           cAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACcaca                                                                            <  1:4860566/150‑1 (MQ=255)
                                                                               tCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGc                                                                        >  2:3679543/1‑150 (MQ=255)
                                                                                 tATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCtt                                                                      >  1:2455404/1‑150 (MQ=255)
                                                                                   tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGTAGCATGATGATGTCAAACTTCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                    <  1:5227243/150‑1 (MQ=255)
                                                                                   tCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                    <  2:4720957/150‑1 (MQ=255)
                                                                                    ccccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCg                                                                   >  1:4230722/1‑150 (MQ=255)
                                                                                      ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                               >  2:2242405/1‑72 (MQ=255)
                                                                                      ccACGAAGTTCTAACCCTGCAAATCAAATTATCTAGATACAATTCTTCGTAGCATGCTGTTGTCAAACATCAGCATGCTCTCGACCCTTATTAATACTTAACACCGCGCTTATCAATGAACACTCACACGCTCACCACAGAGCTTTCGCa                                                                 <  1:4432326/150‑1 (MQ=255)
                                                                                           cAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCa                                                            <  1:2720651/150‑1 (MQ=255)
                                                                                              gTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGaa                                                         >  1:2693865/1‑150 (MQ=255)
                                                                                               tCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAAt                                                        <  2:6514093/150‑1 (MQ=255)
                                                                                                cTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATc                                                       >  2:2879980/1‑150 (MQ=255)
                                                                                                 tAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCa                                                      >  2:5118559/1‑150 (MQ=255)
                                                                                                    cccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGaa                                                   <  2:350894/150‑1 (MQ=255)
                                                                                                         cAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc                                              <  1:4947032/150‑1 (MQ=255)
                                                                                                         cAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc                                              <  1:1170102/150‑1 (MQ=255)
                                                                                                             tCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTAc                                          >  2:6143989/1‑150 (MQ=255)
                                                                                                               aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACAt                                        <  2:5554836/150‑1 (MQ=255)
                                                                                                               aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACAt                                        <  2:1211142/150‑1 (MQ=255)
                                                                                                                aaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATa                                       >  2:6737061/1‑150 (MQ=255)
                                                                                                                  atatCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATaaa                                     <  1:192607/150‑1 (MQ=255)
                                                                                                                   tatCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAat                                    <  1:5329259/150‑1 (MQ=255)
                                                                                                                       tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAATACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                                >  2:6607800/1‑150 (MQ=255)
                                                                                                                       tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaa                                >  2:5398259/1‑150 (MQ=255)
                                                                                                                        aGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaag                               >  2:4043699/1‑150 (MQ=255)
                                                                                                                        aGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaag                               <  2:5755507/150‑1 (MQ=255)
                                                                                                                         gAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaaga                              <  1:2549202/150‑1 (MQ=255)
                                                                                                                          aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaa                             >  2:285768/1‑150 (MQ=255)
                                                                                                                          aaaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaa                             >  1:1569005/1‑150 (MQ=255)
                                                                                                                               aTTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACt                        >  1:4876203/1‑150 (MQ=255)
                                                                                                                                ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTg                       >  1:1334845/1‑150 (MQ=255)
                                                                                                                                ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTg                       >  1:996262/1‑150 (MQ=255)
                                                                                                                                  cttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaa                     <  2:2673659/150‑1 (MQ=255)
                                                                                                                                   ttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaa                    >  1:4498870/1‑150 (MQ=255)
                                                                                                                                    tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                   <  1:3890018/150‑1 (MQ=255)
                                                                                                                                    tcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGaaaa                   <  2:6279735/150‑1 (MQ=255)
                                                                                                                                      ggAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACt                 >  2:3211107/1‑150 (MQ=255)
                                                                                                                                           aTGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaaa            >  1:3374560/1‑150 (MQ=255)
                                                                                                                                             gCTGATGTTTGACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAcc          >  1:102367/1‑150 (MQ=255)
                                                                                                                                             gCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAg                                                     >  1:4116469/1‑107 (MQ=255)
                                                                                                                                               tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt        >  2:836177/1‑150 (MQ=255)
                                                                                                                                               tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCAt        >  1:6822081/1‑150 (MQ=255)
                                                                                                                                                gATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTa                                           <  2:4263886/114‑1 (MQ=255)
                                                                                                                                                    tCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCa   >  2:3310972/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATCAG  >  NC_002947/4817028‑4817326

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