Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,152,038 A→C 5.3% intergenic (+71/‑326) PP_3654 → / → PP_3655 leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,152,0380AC 5.3% 397.5 / 2.6 143intergenic (+71/‑326)PP_3654/PP_3655leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base A (75/53);  new base C (5/5);  total (82/61)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGC  >  NC_002947/4151889‑4152187
                                                                                                                                                          |                                                                                                                                                       
cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                                                         <  2:4541227/150‑1 (MQ=255)
cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                                                         <  1:869141/150‑1 (MQ=255)
cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                                                         <  1:7297220/150‑1 (MQ=255)
 aGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAg                                                                                                                                                        <  2:6815676/150‑1 (MQ=255)
  gTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGCCCCTCGCAGGAGc                                                                                                                                                      >  2:3862951/1‑150 (MQ=255)
   tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCGGCAAAGGTGACGTACCCCGTCGCAGGAGCg                                                                                                                                                     >  1:9051518/1‑150 (MQ=255)
    aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCgg                                                                                                                                                    <  2:6229589/150‑1 (MQ=255)
     ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCggg                                                                                                                                                   <  1:5998802/150‑1 (MQ=255)
      tCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCCCCAGCAAAAGTGGCGTAGCCCGGCGCACGCGCGGGt                                                                                                                                                  >  1:8819213/1‑150 (MQ=255)
      tCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAACTGATCCAGGATCTGTTGCTCCCCATGAAGCAGTATCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGAAGAAGCCCGTCGCAGGatcctgt                                                                                                                                                   >  2:7620560/1‑144 (MQ=255)
        ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCgg                                                                                                                                                    >  1:1982903/1‑146 (MQ=255)
        ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCgg                                                                                                                                                    <  2:1982903/146‑1 (MQ=255)
        ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGttt                                                                                                                                                >  2:3669258/1‑150 (MQ=255)
        ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCCAAGGGGACGGAGCCCGGCGCACGAGCGGGgtt                                                                                                                                                >  1:5234190/1‑150 (MQ=255)
         cAATGTGAAACGCTTGAGCTCGACCTAGGTGATGAAGGAAGTGAACCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCCCCAGCAAAGGTGACGTAGCACGTCCCCGGAGCGGGTTTa                                                                                                                                               >  2:6382727/1‑150 (MQ=255)
          aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTAc                                                                                                                                              <  1:4989152/150‑1 (MQ=255)
                aaaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGTAAACCCCCCCGGCCAAAGTGGCGTAGCCCCTCGCAGGAGCGGGGTTAGCCCCGa                                                                                                                                        >  1:1002742/1‑150 (MQ=255)
                aaaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGa                                                                                                                                        >  2:6624366/1‑150 (MQ=255)
                 aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGGAAACCCCCCCGGCAAAAGTGTCGGAGCCCCGCGCAGGCGGGGGGTTTAccccgac                                                                                                                                        >  1:2445018/1‑146 (MQ=255)
                 aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTATCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTTACGTAGCCCCTCGCAGGAGCGGGTTTACCCGCGac                                                                                                                                       >  2:6109632/1‑149 (MQ=255)
                          gCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCAGCAAACGTGACGTAGCCCCTCGCAGGAGCCGGTTTACCCCCCACGAggccggc                                                                                                                              >  2:2713396/1‑150 (MQ=255)
                           cTCGCCCTTGGGGATGACGGACGTGATCCTGGATCGGTTGCTCCCCATGTAGCTGTTGCCTCTTGTGCCGGCCCCTTCGCGGTTAAACCCGCGCCTGCAACGGTGACTTAGCCCTTCGAAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                                                                             <  2:6973744/150‑1 (MQ=37)
                            tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTATCCCCTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGccc                                                                                                                            >  1:8215553/1‑150 (MQ=255)
                              gACCTTGTTGATGACGGAAGTGAGCCCGTATGGGTTGCCCCCCATGTCGGTGTTTCCTCTTGGGCCGGCCCCTGCGCGGGTAAACCGGCGCCTGCAAAGGTGACGTAGCCCGCCGCTGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAg                                                                                                                          <  2:293719/150‑1 (MQ=255)
                                cgttGGTGATGAAGTAAGTGATCCAGTATGGGTTGCTCCCCATGTAGCTGTTTCCTCTTGGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCaaga                                                                                                                        <  1:7520687/148‑4 (MQ=38)
                                       gatgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAAGGTGAAGTAGCCCGTCGCCGGAGCGGGTTTACCCGCGAAGCGGCCGGCCCAGGCAACCatt                                                                                                                 >  2:4244147/1‑150 (MQ=39)
                                         tgaCGGATGTGATCCCGGATCGGTTGCTCCCCATTTTGCTGTTTCCTCTTGGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAACGGTGACGTAGCCGGTCGAAGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCa                                                                                                               <  2:1150216/150‑1 (MQ=37)
                                          gaaggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCtttaaa                                                                                                              >  2:7500781/1‑147 (MQ=255)
                                              gcagTGACCCAGGATCGGTTGCGCCCCATTTAGCTGTTTCCTGTTTTGCCGGCCCCTTCGCGGGTAACCCCGCGCCTGCAAAGTTGACGTAGCCGGTCGCTGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGc                                                                                                          <  2:124105/148‑1 (MQ=255)
                                                agTGATCCAGGATGGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGccc                                                                                                        <  2:4331817/150‑1 (MQ=255)
                                                agTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGccc                                                                                                        <  2:8933203/150‑1 (MQ=255)
                                                       cAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGGAGCCCCGCGCAGGAGCGGGGTTACCCCCGCAGAGGCCGGCCCCGGCCACCAATCATCTGCCCTGaaaaa                                                                                                 >  1:5291368/1‑149 (MQ=255)
                                                         ggATCGGTTGCTCCCCCTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAACCCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgc                                                                                               <  2:4254175/150‑1 (MQ=255)
                                                         ggAGCGGGTGCGCCCCCGGTAGCGGTGTCCTCTGGGGCCGGCCCCGTCGCGGGTACCCCCGCGCCTGCAACGGTGCCGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgc                                                                                               <  2:9086547/150‑1 (MQ=37)
                                                                ttGCTCCCCTTGTGTTGTTCCCTTGTGGTGCGGCCCCCTGCGGGGGTAAACCCGCGCCTGCAAGGGTGTGGTGGCCGGTCGATGGCGGGTGTTTACCCGCGATGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATc                                                                                        <  1:7205098/150‑1 (MQ=255)
                                                                   cTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCa                                                                                     <  2:7015947/150‑1 (MQ=255)
                                                                    tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCCAAGGTGACGTAGCCCGTCGCAGGAGCGGGGTTTACCGCGAAGAGGCCGGCCCACGCAACCATTAATATGCCCTGTAAACGCGCCCATCGCAc                                                                                   >  2:6928437/1‑150 (MQ=255)
                                                                     ccccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGACGTAGCCCGTAGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCACGCAACCCTTCATCAGCCCAGAAAACGCGCACATAGCAcc                                                                                  >  1:6758812/1‑150 (MQ=38)
                                                                        cATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCCTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCACGCAACCATTCAAATGCCCAGAAAAACCGCACAACGCCCCAAc                                                                               >  2:3464141/1‑150 (MQ=255)
                                                                         aTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCGGCCAAGGTGAACGACCCCGTCGCCGGCGCGGGGTTAACCGCCAAGAGGCCGGCCCAGGCAAACGATCAACTGCCCAGAAAACGCGCCCAACGcacccacc                                                                              >  2:3225333/1‑150 (MQ=255)
                                                                             aGCTGTTTCCTCTTGGGCGGGCCCCTTCGCGGTTAAACCCGCGCCTGCAAAGGTGACGTAGCCGGTCGCTGGAGGGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGc                                                                          <  1:5142114/150‑1 (MQ=255)
                                                                                tGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCCTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCAGAAAACGCGCACATCGCACCAACCCAGCCCa                                                                       >  2:5329487/1‑150 (MQ=255)
                                                                                tGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCAGCCAAGGTGACGTAGCCCCTCGCAGGAGCCGGTTTAACCGCGCCCGGGCCGGCCCCGGCCACCATTCATCAGCCCGGAAAAAGCGCACAACGCacccacccccgccc                                                                       >  2:8419475/1‑144 (MQ=255)
                                                                                tGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGTAACCCCCCCCGGCCAAAGTTACGTTGCCCCGCGCAGGAGCGGGGTTAACCGCCCAGAGGCCGGCCCAGGCCACCAATAAACTGCCCAGAAAAAGCCCCCAACGCACCCAAccccccca                                                                        >  2:2276001/1‑145 (MQ=255)
                                                                                  tttCCTCTTGTGCCGGCCCGTTCGGGTTAAAACCCGCGCCTCAAACGGTGACGTCGCCGGTCGCCGGTGGGGTTTAACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCAcc                                                                     <  1:6956339/150‑1 (MQ=255)
                                                                                        tctTGTGCCGGCCCCTTCGCGGGTTAAACCCCCCCGGCAAAGGTGTCCTAGCCCCGCGCAGGAGCCGGGTTTACCGCCAACAGGCCGGCCCACGCAACCATTAAAATGCCCGGAAAAACCGCACAACCCACCAACCCAGcccccccgggt                                                                >  2:2577219/1‑145 (MQ=255)
                                                                                           tggcccGCCCCCTGCGGGGGTAACCCCGCGCCTGAAATGGTGTCGTCGCCCGTCGCCGGTGGGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCg                                                             <  1:4697403/146‑1 (MQ=255)
                                                                                              gccgccCCCTGCGGGGTAAACCCCGCGCCTCAACTGGTGCGGCCGCCGGCCGCTGGCGGGGTTTAACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                                                          <  1:800585/150‑1 (MQ=255)
                                                                                              gCCGGCCCCTGCGCGGGTAACCCCGCGCCGGCACCGGTGAGGTCGCCCGTCGCAGGCGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                                                          <  1:7046824/150‑1 (MQ=255)
                                                                                               ccGGCCTCTTCGCGGGTAAAACCCCAACTGCGCCGGGCGACGTCCCCGTTGCAGGCGCGGGGTTACCCGCGGAGGGGCCCGCCCAGGCAACCTTTCATCTGCCCTGCAAACGCGCCCCACGCACCAACCCAGCCCACCGGGATCCGGTTg                                                         >  2:7716469/1‑150 (MQ=12)
                                                                                               ccGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGACAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTg                                                         <  2:4588867/150‑1 (MQ=255)
                                                                                                 ggCCCCTTCGCGGGTAAACCCGCCCCAGCAAAAGTGACGTAGCCCCTCGCAGGAGCGGGTTTAACCGCGAAGAGGCCGGCCCACGCAACCAATAAAATGCCCAGAAAAAGCGCACAAAGCAACAAACCACCCCACCGTGATACCGTTTcc                                                       >  2:3256854/1‑150 (MQ=255)
                                                                                                    cccTTCGCGGGTAAACCCGCGCCTGCAAAGGTTACCTATCCCCTCGCAGGAGCCGGTTTACCCGCGAAGAGGCCCGCCCACGCAACCAATAATCTGCCCAGAAAAAGCGCACAACCCACCAACCCACCCCACACAGAACCCGGTGCCAcc                                                    >  2:8797588/1‑150 (MQ=19)
                                                                                                       ttCGCGGGTAAACCCGCTCCTGCGACGGGCTACGGCACCTTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCCACCCACCCTGTTCCGGGTGcccccccc                                                 >  1:2962201/1‑150 (MQ=9)
                                                                                                             gggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCcgccgc                                            <  2:5121867/150‑1 (MQ=255)
                                                                                                              ggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc                                           >  1:2221671/1‑150 (MQ=255)
                                                                                                               gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCGCCGccc                                          >  1:8582093/1‑150 (MQ=255)
                                                                                                                tAACCCCGCGCCTTAAATGGTGAGGTCGCCGGGCGCTGGAGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGcccc                                         <  1:120241/150‑1 (MQ=255)
                                                                                                                tAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGcccc                                         <  2:7345321/150‑1 (MQ=255)
                                                                                                                     cccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTa                                     >  1:3161330/1‑150 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAcccc                                                  >  1:4617973/1‑136 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAcccc                                                  <  2:4617973/136‑1 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                    >  1:7885986/1‑150 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                    >  2:1308513/1‑150 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                    >  1:5985763/1‑150 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCGAc                                    >  1:1170071/1‑150 (MQ=255)
                                                                                                                      ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCCCCGccccgcgac                                    >  1:9139341/1‑145 (MQ=255)
                                                                                                                       cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAcc                                   >  2:1401055/1‑150 (MQ=255)
                                                                                                                       cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCCCCGccccgcgaac                                   >  2:520378/1‑144 (MQ=255)
                                                                                                                        gcgcCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAccc                                  >  1:5377163/1‑150 (MQ=255)
                                                                                                                        gcgcCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAccc                                  >  1:5377185/1‑150 (MQ=255)
                                                                                                                         cgcCTGAAACGGTGCCGTCCCCCGGCGCCGGCGGGGGTTTACCCGCGAAGAGGCCGGCCCATGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAcccc                                 <  1:2605584/150‑1 (MQ=255)
                                                                                                                            cTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCGCCGCCCCTCTACCCCAAg                              >  1:8539405/1‑150 (MQ=255)
                                                                                                                            cTGAAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAg                              <  2:1763952/150‑1 (MQ=255)
                                                                                                                              gCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGc                            >  1:8380154/1‑150 (MQ=255)
                                                                                                                                   aGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTc                        <  2:9063049/150‑1 (MQ=255)
                                                                                                                                   aGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTc                        <  2:5080936/150‑1 (MQ=255)
                                                                                                                                    ggTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTca                       <  2:1002742/150‑1 (MQ=255)
                                                                                                                                    ggTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTca                       <  2:2213826/150‑1 (MQ=255)
                                                                                                                                    ggTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATACGGTTGCCACCCCCGCCGCCCCTCTAACCCAAGGCCTTca                       >  1:7962255/1‑150 (MQ=255)
                                                                                                                                     gTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTcag                      <  2:3597245/150‑1 (MQ=255)
                                                                                                                                      tGAGGTAGCCGGTCGCTGGAGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTcagc                     <  1:6881679/150‑1 (MQ=255)
                                                                                                                                      tGACGTCGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTcagc                     <  1:5072900/150‑1 (MQ=255)
                                                                                                                                      tGACGTCGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTcagc                     <  1:3669258/150‑1 (MQ=255)
                                                                                                                                      tGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCctc                                       <  1:4405909/132‑1 (MQ=255)
                                                                                                                                      tGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCctc                                       >  2:4405909/1‑132 (MQ=255)
                                                                                                                                      tGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCCCCGCCCCCCTACCCCAAGGCCTTcagc                     >  1:3707087/1‑150 (MQ=255)
                                                                                                                                          gTAGCCCGTCGCAGGAGCGGGTTTATCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCGTcagccagc                 >  2:53126/1‑150 (MQ=255)
                                                                                                                                          gTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTcagccagc                 >  2:4799019/1‑150 (MQ=255)
                                                                                                                                            aGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAccc                                  <  2:8014872/131‑1 (MQ=255)
                                                                                                                                            aGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTg               >  1:2589348/1‑150 (MQ=255)
                                                                                                                                            aGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTg               <  1:6039593/150‑1 (MQ=255)
                                                                                                                                                    gCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACtctc       <  1:6872496/150‑1 (MQ=255)
                                                                                                                                                      aGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCGCCGCCCCGCTACCCCAAGGCCGTCAGCCAGCTGGCACTCTCCt     >  1:3698792/1‑150 (MQ=255)
                                                                                                                                                      aGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCt     >  1:6675937/1‑150 (MQ=255)
                                                                                                                                                       ggAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCtt    >  2:7727044/1‑150 (MQ=255)
                                                                                                                                                       ggAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCtt    <  2:6327150/150‑1 (MQ=255)
                                                                                                                                                       ggAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCtt    <  2:2604612/150‑1 (MQ=255)
                                                                                                                                                       ggAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCtt    <  2:3866857/150‑1 (MQ=255)
                                                                                                                                                        gAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAccc                                                   >  2:3672648/1‑102 (MQ=255)
                                                                                                                                                          |                                                                                                                                                       
CAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGC  >  NC_002947/4151889‑4152187

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