Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,152,000 A→T 6.9% intergenic (+33/‑364) 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,0000AT 6.9% 430.1 / 5.7 150intergenic (+33/‑364)PP_3654/PP_3655leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base A (62/74);  new base T (6/4);  total (72/78)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

TGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCC  >  NC_002947/4151852‑4152148
                                                                                                                                                      |                                                                                                                                                         
tGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTaa                                                                                                                                                          >  2:5844740/1‑150 (MQ=255)
tGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTaa                                                                                                                                                          <  2:2328541/150‑1 (MQ=255)
tGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTaa                                                                                                                                                          >  1:3597245/1‑150 (MQ=255)
  ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAc                                                                                                                                                       <  1:1298693/150‑1 (MQ=255)
  ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAc                                                                                                                                                       <  1:1666801/150‑1 (MQ=255)
  ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAc                                                                                                                                                       <  1:2027854/150‑1 (MQ=255)
  ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAc                                                                                                                                                       <  1:2027873/150‑1 (MQ=255)
    ggCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAccc                                                                                                                                                     >  2:4697144/1‑150 (MQ=255)
         ttCAGAAGCTGTATGACGAGCACCTCACCAGTATTCCCAATCTGAAACGCTTGAGCTCGACCTTGGAGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCAGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                                <  2:5345601/150‑3 (MQ=255)
           cAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTg                                                                                                                                              <  1:72644/150‑1 (MQ=255)
            aGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGc                                                                                                                                             >  1:1078563/1‑150 (MQ=255)
            aGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGAAGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaaggg                                                                                                                                             <  2:6973516/150‑6 (MQ=255)
            aGAAGCTGTATGACGAGCACCTGACCAGTAGTCCCAATGTGAAACGCTTTAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGc                                                                                                                                             <  1:7403918/150‑1 (MQ=255)
               aGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCaaa                                                                                                                                         <  2:4074563/150‑1 (MQ=255)
                gCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                        <  2:3141685/150‑1 (MQ=255)
                gCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                        <  1:6094551/150‑1 (MQ=255)
                  tGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                      <  2:1078563/150‑1 (MQ=255)
                   gTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTg                                                                                                                                     <  1:4177964/150‑1 (MQ=255)
                   gTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGGAAACCCCGCCAGCCAAGGTg                                                                                                                                     >  1:6327150/1‑150 (MQ=255)
                    tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGa                                                                                                                                    <  1:5851714/150‑1 (MQ=255)
                    tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGa                                                                                                                                    <  1:5851571/150‑1 (MQ=255)
                    tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCGCCAGCAAAGGtta                                                                                                                                    >  1:6233328/1‑148 (MQ=255)
                    tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGTAACCCCCGCCGGGCAAGGTgg                                                                                                                                    >  1:2523229/1‑149 (MQ=255)
                      tGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACg                                                                                                                                  >  1:5323884/1‑150 (MQ=255)
                      tGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACg                                                                                                                                  >  2:7455841/1‑150 (MQ=255)
                         cGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAg                                                                                                                               >  1:6781478/1‑150 (MQ=255)
                                 tGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGc                                                                                                                       <  1:7473190/150‑1 (MQ=255)
                                  gACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGCCCGTCGcc                                                                                                                      >  2:7152095/1‑149 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGccg                                                                                                                     >  2:6039593/1‑148 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                     <  1:2904340/150‑1 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                     <  2:1961649/150‑1 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                     <  1:7785882/150‑1 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                     <  1:3096266/150‑1 (MQ=255)
                                   aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTCCCTCTTGTGCCGGCCCCTGCGCGGGTAAACCCCCGCCAGCAAAGGTGACGTAGCCCCGCGCAg                                                                                                                     >  2:7205098/1‑150 (MQ=255)
                                     cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                  <  1:869141/150‑1 (MQ=255)
                                     cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                  <  2:4541227/150‑1 (MQ=255)
                                      aGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAg                                                                                                                 <  2:3957325/150‑1 (MQ=255)
                                       gTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGCCCCTCGCAGGAGc                                                                                                               >  2:3862951/1‑150 (MQ=255)
                                        tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCg                                                                                                              <  1:9035164/150‑1 (MQ=255)
                                        tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCGGCAAAGGTGACGTACCCCGTCGCAGGAGCg                                                                                                              >  1:9051518/1‑150 (MQ=255)
                                          ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCggg                                                                                                            <  1:955701/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                                                                                                             <  2:1982903/146‑1 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCgg                                                                                                             >  1:1982903/1‑146 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGttt                                                                                                         <  1:1563920/150‑1 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGttt                                                                                                         <  2:2994604/150‑1 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCCAAGGGGACGGAGCCCGGCGCACGAGCGGGgtt                                                                                                         >  1:5234190/1‑150 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAAGGTGACGTAGCCCCTCGCAGGAGCGGGttt                                                                                                         >  1:6789662/1‑150 (MQ=255)
                                             ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGttt                                                                                                         >  2:6881679/1‑150 (MQ=255)
                                              cAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTa                                                                                                        >  2:9015636/1‑150 (MQ=255)
                                              cAATGTGAAACGCTTGAGCTCGACCTAGGTGATGAAGGAAGTGAACCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCCCCAGCAAAGGTGACGTAGCACGTCCCCGGAGCGGGTTTa                                                                                                        >  2:6382727/1‑150 (MQ=255)
                                               aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTAc                                                                                                       <  1:4989152/150‑1 (MQ=255)
                                               aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTAc                                                                                                       <  1:4989377/150‑1 (MQ=255)
                                                   tgAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCcgc                                                                                                   <  2:4474694/150‑1 (MQ=255)
                                                   tgAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAAGGGGACGTAGCCCCTCGCAGGAGCGGGTTTAccccc                                                                                                   >  2:4857653/1‑148 (MQ=255)
                                                    gAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAAAACCCGCCAGCCAAGGTTTCGTAGCCCCGCGCCGGAGCCGGGGTTAccccg                                                                                                   >  1:2345495/1‑148 (MQ=255)
                                                     aaaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGTAAACCCCCCCGGCCAAAGTGGCGTAGCCCCTCGCAGGAGCGGGGTTAGCCCCGa                                                                                                 >  1:1002742/1‑150 (MQ=255)
                                                     aaaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGa                                                                                                 >  2:6624366/1‑150 (MQ=255)
                                                      aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCAGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCCaa                                                                                                >  2:698333/1‑150 (MQ=255)
                                                      aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGGAAACCCCCCCGGCAAAAGTGTCGGAGCCCCGCGCAGGCGGGGGGTTTAccccgac                                                                                                 >  1:2445018/1‑146 (MQ=255)
                                                      aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTATCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTTACGTAGCCCCTCGCAGGAGCGGGTTTACCCGCGac                                                                                                >  2:6109632/1‑149 (MQ=255)
                                                           ttGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCGCCGGCAAAGGTTACCTAGCCCGGCCCAGGAGCGGGTTTAACCGCGAAAAggc                                                                                           >  2:430225/1‑150 (MQ=255)
                                                               gCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCAGCAAACGTGACGTAGCCCCTCGCAGGAGCCGGTTTACCCCCCACGAggccggc                                                                                       >  2:2713396/1‑150 (MQ=255)
                                                                 tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCAGCAAACGTGGCCTAGCCCCTCGCACGAGCGGGGTTACCCCCCAACAGGCCGGccc                                                                                     >  1:7371767/1‑150 (MQ=255)
                                                                   gACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCat                                                                                   >  2:2605584/1‑149 (MQ=255)
                                                                     cgttGGTGATGAAGTAAGTGATCCAGTATGGGTTGCTCCCCATGTAGCTGTTTCCTCTTGGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCaaga                                                                                 <  1:7520687/148‑4 (MQ=38)
                                                                           tgatgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCCTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAAccat                                                                           >  2:6588390/1‑148 (MQ=255)
                                                                            gatgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAAGGTGAAGTAGCCCGTCGCCGGAGCGGGTTTACCCGCGAAGCGGCCGGCCCAGGCAACCatt                                                                          >  2:4244147/1‑150 (MQ=39)
                                                                              tgaCGGATGTGATCCCGGATCGGTTGCTCCCCATTTTGCTGTTTCCTCTTGGGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAACGGTGACGTAGCCGGTCGAAGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCa                                                                        <  2:1150216/150‑1 (MQ=37)
                                                                               gaaggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTGAACCCCCCCCAGCAAACGTTAACTAGCCCGGCCCAGGAGCGGGTTTAACCCCCAAGAGGCCGGCCCCGGCAAccaataaa                                                                       >  2:848820/1‑144 (MQ=255)
                                                                                   gcagTGACCCAGGATCGGTTGCGCCCCATTTAGCTGTTTCCTGTTTTGCCGGCCCCTTCGCGGGTAACCCCGCGCCTGCAAAGTTGACGTAGCCGGTCGCTGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGc                                                                   <  2:124105/148‑1 (MQ=255)
                                                                                   gaagTGATCCAGGATGGTTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGAAACGGTTACGTAGCCGGTCGCAGGAGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGc                                                                   <  2:124090/150‑1 (MQ=38)
                                                                                     agTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGccc                                                                 <  2:8933203/150‑1 (MQ=255)
                                                                                      gTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCt                                                                <  2:3631111/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)
                                                                                                         tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGTCCTAGCCCGTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCAGGCAACCATTCATCTGCCCGGAAAAACCGCACAtaccac                                             >  2:7828780/1‑145 (MQ=255)
                                                                                                             cATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCCTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCACGCAACCATTCAAATGCCCAGAAAAACCGCACAACGCCCCAAc                                        >  2:3464141/1‑150 (MQ=255)
                                                                                                              aTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCGGCCAAGGTGAACGACCCCGTCGCCGGCGCGGGGTTAACCGCCAAGAGGCCGGCCCAGGCAAACGATCAACTGCCCAGAAAACGCGCCCAACGcacccacc                                       >  2:3225333/1‑150 (MQ=255)
                                                                                                                 tAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCGCCAGCAAAGGTTAAGTAGCCCGGCGCCGGGGCGGGTTTACCCGCGTAGGGGCCGGCCCAGGCAACCATTCAACTGGCCTGAAAACGCGCACATAGCCCCAACCCac                                    >  2:4591751/1‑149 (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)
                                                                                                                                tgggcgggCCCCTTCGCGGTTAACCCCGCGCCTGCAAGGGTGACGTCGCCGGTGGCAGGCGCGGTTTTCCCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCg                      <  1:4697144/144‑1 (MQ=255)
                                                                                                                                   gccgccCCCTGCGGGGTAAACCCCGCGCCTCAACTGGTGCGGCCGCCGGCCGCTGGCGGGGTTTAACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                   <  1:800585/150‑1 (MQ=255)
                                                                                                                                    ccGGCCTCTTCGCGGGTAAAACCCCAACTGCGCCGGGCGACGTCCCCGTTGCAGGCGCGGGGTTACCCGCGGAGGGGCCCGCCCAGGCAACCTTTCATCTGCCCTGCAAACGCGCCCCACGCACCAACCCAGCCCACCGGGATCCGGTTg                  >  2:7716469/1‑150 (MQ=12)
                                                                                                                                     cGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGc                 <  2:7371767/150‑1 (MQ=255)
                                                                                                                                      ggCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCATTCATCTGCCCTGAAAACGCGCACATAGCACCAACCCACCCCCCCCTGATCCGGTTGcc                >  1:8399215/1‑150 (MQ=38)
                                                                                                                                               gcgGGTAAACCCGCTCCTGCGACGGGCTACGTCACCTTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAcccccccc       >  1:3066249/1‑150 (MQ=18)
                                                                                                                                                  gggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCcgccgc     <  2:5121867/150‑1 (MQ=255)
                                                                                                                                                   ggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc    >  1:2221671/1‑150 (MQ=255)
                                                                                                                                                    gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc   >  1:2604612/1‑150 (MQ=255)
                                                                                                                                                    gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc   <  2:6233328/150‑1 (MQ=255)
                                                                                                                                                     tAACCCCGCGCCTTAAATGGTGAGGTCGCCGGGCGCTGGAGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGcccc  <  1:120241/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                         
TGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCC  >  NC_002947/4151852‑4152148

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: