Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,198 C→G 7.7% intergenic (‑140/+36) 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,1980CG 7.7% 284.1 / 11.8 91intergenic (‑140/+36)PP_4606/PP_4607ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base G (3/4);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

CATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCT  >  NC_002947/5229056‑5229346
                                                                                                                                              |                                                                                                                                                     
cATGGACGCAGACCTGTGATGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACTGGGTTGCCACGTGGTGTTACCTCTCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCGTCTTCGCGGGTATACCCGCCCCTACggggg                                                                                                                                                 <  2:64695/148‑1 (MQ=255)
cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                                 <  3:273900/148‑1 (MQ=255)
cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                                 <  7:115420/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                                                                                                                           <  6:221896/149‑1 (MQ=255)
                     gggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGc                                                                                                                           <  5:410509/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                                                                                          >  7:32872/1‑130 (MQ=255)
                                                                         gCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                          <  8:32872/130‑1 (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:319416/1‑149 (MQ=255)
                                                                               aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCa                                                                 >  4:147767/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                                                                <  4:24157/136‑1 (MQ=255)
                                                                                             tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCAt                                                                >  3:24157/1‑136 (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           >  6:9498/1‑149 (MQ=255)
                                                                                                                                     ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc           >  1:124017/1‑149 (MQ=255)
                                                                                                                                     ccGCCCCTACGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGACAAACGCACCAACACGGTCTGCCGGGAAAAACCCATCACCTTGACCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGTGAAGCTGc           >  6:130105/1‑149 (MQ=255)
                                                                                                                                         cccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCt       >  8:370157/1‑149 (MQ=255)
                                                                                                                                          ccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                          <  5:273965/129‑1 (MQ=255)
                                                                                                                                          ccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCgg                          >  6:273965/1‑129 (MQ=255)
                                                                                                                                           cTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCGGCGGCTCGGCCACCGGGAAGCTGCCGCTGa     >  2:195456/1‑149 (MQ=255)
                                                                                                                                              cGGGGGGTTTACCCGTGAAGCGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCACCACGGGCGGCCGTGCAAAACCCATCACCTGGCAGAGCCAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGATGCCGCTGACCt  >  8:270047/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCT  >  NC_002947/5229056‑5229346

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: