Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,201:1 +C 8.0% intergenic (‑143/+33) PP_4606 ← / ← PP_4607 ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,229,2011.C 8.0% 344.8 / 24.5 88intergenic (‑143/+33)PP_4606/PP_4607ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base . (42/39);  new base C (3/4);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTG  >  NC_002947/5229056‑5229348
                                                                                                                                                  |                                                                                                                                                   
cATGGACGCAGACCTGTGATGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACTGGGTTGCCACGTGGTGTTACCTCTCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCGTCTTCGCGGGTATACCCGCCCCTACggggg                                                                                                                                                   <  2:64695/148‑1 (MQ=255)
cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                                   <  7:115420/148‑1 (MQ=255)
cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                                   <  3:273900/148‑1 (MQ=255)
 aTGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgggggg                                                                                                                                                  <  8:397975/148‑1 (MQ=255)
 aTGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGt                                                                                                                                                 <  4:214644/149‑1 (MQ=255)
  tGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGtt                                                                                                                                                >  5:154477/1‑149 (MQ=255)
     acacAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTAc                                                                                                                                             <  8:147354/146‑1 (MQ=255)
      cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTAcc                                                                                                                                            <  2:4289/149‑1 (MQ=38)
        cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCAGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCg                                                                                                                                          <  6:387908/149‑1 (MQ=37)
         aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         <  2:177222/149‑1 (MQ=255)
         aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         >  5:387372/1‑149 (MQ=255)
         aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCgc                                                                                                                                         <  6:35441/149‑2 (MQ=38)
         aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACCCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         >  4:278433/1‑149 (MQ=255)
         aGACCTGGGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGTTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         <  3:209696/149‑1 (MQ=255)
          gACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTg                                                                                                                                        <  7:271922/149‑1 (MQ=255)
           aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTg                                                                                                                                        >  5:221896/1‑148 (MQ=255)
             cTGTGAGGGGGTGATTATTGCGAATCAATTGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGCCAGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAag                                                                                                                                     <  2:29747/149‑1 (MQ=255)
               gtgCGGGGGTGATTCTTGCGACTCACTCGCACTTGCACGCCTTCCACGGATGCCCCGTGCTGGTCCCTCCCCGGCGCCGGGAAATATTGTTGCTGTGGCGGCCTCTGCGCGGGGAAACCCGCCCCTACGGGGGGTTTACCCGGGAagag                                                                                                                                   <  2:28490/149‑1 (MQ=255)
                tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAgg                                                                                                                                  >  6:84426/1‑149 (MQ=255)
                tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCCGTTACCTCACCTGTGCCCGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAagag                                                                                                                                   <  1:431053/148‑1 (MQ=255)
                tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACCCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAaga                                                                                                                                    <  3:278433/147‑1 (MQ=255)
                 ctggggggCATTATTCCGAACCAATCGCACTTCCACTCATTCCACGCATGCGACGTCCTGTTCCCTCACCTGTGCCTGCAAATATTTTTGCTGTGGGGGCCTCTTCGCGGTTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGc                                                                                                                                 <  4:23611/147‑1 (MQ=255)
                  aGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGc                                                                                                                                 <  5:114380/148‑1 (MQ=255)
                    ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                                              >  4:78639/1‑149 (MQ=255)
                    ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                                              <  1:429704/149‑1 (MQ=255)
                     gggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGc                                                                                                                             <  5:410509/149‑1 (MQ=255)
                     gggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGc                                                                                                                             <  6:221896/149‑1 (MQ=255)
                          attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                        >  8:394812/1‑149 (MQ=255)
                           ttattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGc                                                                                                                       <  4:324241/149‑1 (MQ=255)
                              ttccGAAACAATCGCAATTTAACTCATAACACCGATGCCACGTGCTGTTAACACACCTGTCCCTGGAAATATTTTTCCTCTGGCGGCATCTTAGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGAGAAGAGGCTCGCACAGGCaa                                                                                                                     >  4:297111/4‑148 (MQ=255)
                              ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGa                                                                                                                                       >  3:342314/1‑130 (MQ=255)
                              ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGa                                                                                                                                       <  4:342314/130‑1 (MQ=255)
                              ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCACCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                        <  1:255791/145‑1 (MQ=255)
                              ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCACCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                        >  2:255791/1‑145 (MQ=255)
                               tGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGCCAGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGa                                                                                                                                       >  7:126918/1‑129 (MQ=21)
                               tGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGCCAGCCTCTTCACGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGa                                                                                                                                       <  8:126918/129‑1 (MQ=21)
                                  gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTGACCCGGGAAGAGGCTGGCACAGGCAACTAcc                                                                                                                >  7:324821/1‑149 (MQ=255)
                                    atcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGa                                                                                                              <  7:36737/149‑1 (MQ=255)
                                       aatcGCACTTTCACTCATTACACGGATGCCACGTGATGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAAc                                                                                                            >  1:283611/1‑148 (MQ=255)
                                          cGCACTTTCACTCATTACACGGATGCCACGTGCTGTTAGCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAg                                                                                                         <  6:3040/148‑1 (MQ=255)
                                            cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAAcc                                                                                                           >  1:164107/1‑144 (MQ=255)
                                            cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAAcc                                                                                                           <  2:164107/144‑1 (MQ=255)
                                                  cACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAAc                                                                                                 <  8:7321/148‑1 (MQ=255)
                                                     tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         >  3:152093/1‑105 (MQ=255)
                                                     tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                         <  4:152093/105‑1 (MQ=255)
                                                       aTTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAATGcccc                                                                                            <  6:154477/148‑1 (MQ=255)
                                                       aTTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                            >  2:379124/1‑148 (MQ=255)
                                                           cacGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAaca                                                                                        <  8:344498/148‑1 (MQ=255)
                                                            acGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACg                                                                                      >  4:283270/1‑149 (MQ=255)
                                                              ggATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGt                                                                                    <  4:390167/149‑1 (MQ=255)
                                                                   ctaCGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCACAACACGGTCTGcc                                                                               >  8:247186/3‑149 (MQ=255)
                                                                    cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                              >  6:415999/1‑149 (MQ=255)
                                                                         gCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                            <  8:32872/130‑1 (MQ=255)
                                                                         gCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                            >  7:32872/1‑130 (MQ=255)
                                                                            gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAaca                                                                                        <  6:392050/131‑1 (MQ=255)
                                                                            gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAaca                                                                                        >  5:392050/1‑131 (MQ=255)
                                                                            gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                                                      >  8:95042/1‑149 (MQ=255)
                                                                               aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCa                                                                   >  4:147767/1‑149 (MQ=255)
                                                                               aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCa                                                                   >  4:319416/1‑149 (MQ=255)
                                                                                 cTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGCACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                  >  7:346583/1‑148 (MQ=255)
                                                                                         tgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGcc                                                         <  6:387372/149‑1 (MQ=255)
                                                                                            cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTg                                                           >  4:276609/1‑144 (MQ=255)
                                                                                            cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTg                                                           <  3:276609/144‑1 (MQ=255)
                                                                                             tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                  >  3:24157/1‑136 (MQ=255)
                                                                                             tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                  <  4:24157/136‑1 (MQ=255)
                                                                                               gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                     <  3:78639/147‑1 (MQ=255)
                                                                                                 tatatTTTTGCCGGGGCGGCCGCTTCGGGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAt                                                 <  2:336236/148‑1 (MQ=255)
                                                                                                   tatTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCAcc                                                              <  4:8039/134‑1 (MQ=255)
                                                                                                   tatTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCAcc                                                              >  3:8039/1‑134 (MQ=255)
                                                                                                   tatTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATcc                                               <  5:84426/149‑1 (MQ=255)
                                                                                                     tttttGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCCCTACGGGGGGTTTAACCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCAGCCGTGAAAAACCCATAACCATTCCCAGCCAAACCAg                                             >  6:369681/1‑149 (MQ=255)
                                                                                                           cTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca                                        <  1:379124/148‑1 (MQ=255)
                                                                                                            tgtgGCGGCCGCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCcagc                                      <  8:324821/149‑1 (MQ=255)
                                                                                                             gtgGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaaaa                                                                       <  6:196329/115‑1 (MQ=255)
                                                                                                             gtgGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGaaaaa                                                                       >  5:196329/1‑115 (MQ=255)
                                                                                                                           gcgGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAcc                       >  5:303171/1‑149 (MQ=255)
                                                                                                                             gggTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                     >  3:393085/1‑149 (MQ=255)
                                                                                                                              ggTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAACCCCATCACCTTGCCCAGCGAAGCCAGCGCCAGCAGCGGCTCGGCCGCCgg                     >  3:285525/1‑148 (MQ=255)
                                                                                                                              ggTAAACCCCCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCACCTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAACCCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                     >  7:191931/1‑148 (MQ=38)
                                                                                                                              ggTAAACCCCCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACTACAGAACCCGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGACCCCCgg                     >  3:284866/1‑148 (MQ=37)
                                                                                                                               gTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                    >  3:140245/1‑148 (MQ=255)
                                                                                                                                  aaCCCGCCCCTACGCGTGGTTTACCCGTGAAGAGGCGGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAg                 <  4:402679/148‑1 (MQ=255)
                                                                                                                                   aCCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGACCCAGGGTCAAACGCCCCAACACGGGCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCACCGGCTCGGCCACCGGGAAGCt               >  3:277896/1‑149 (MQ=255)
                                                                                                                                    cccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTg              >  4:284164/1‑149 (MQ=255)
                                                                                                                                     ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc             >  1:124017/1‑149 (MQ=255)
                                                                                                                                     ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc             >  6:9498/1‑149 (MQ=255)
                                                                                                                                     ccGCCCCTACGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGACAAACGCACCAACACGGTCTGCCGGGAAAAACCCATCACCTTGACCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGTGAAGCTGc             >  6:130105/1‑149 (MQ=255)
                                                                                                                                         cccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCt         >  8:370157/1‑149 (MQ=255)
                                                                                                                                          ccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                            >  6:273965/1‑129 (MQ=255)
                                                                                                                                          ccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                            <  5:273965/129‑1 (MQ=255)
                                                                                                                                           cTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCGGCGGCTCGGCCACCGGGAAGCTGCCGCTGa       >  2:195456/1‑149 (MQ=255)
                                                                                                                                              cGGGGGGTTTACCCGTGAAGCGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCACCACGGGCGGCCGTGCAAAACCCATCACCTGGCAGAGCCAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGATGCCGCTGACCt    >  8:270047/1‑149 (MQ=255)
                                                                                                                                               ggggggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCtt   >  3:298853/1‑149 (MQ=255)
                                                                                                                                               ggggggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCtt   >  4:103639/1‑149 (MQ=255)
                                                                                                                                                gggggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTg  >  6:186380/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTG  >  NC_002947/5229056‑5229348

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