Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,152,000 A→T 6.6% 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.6% 444.0 / 7.8 143intergenic (+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 (61/70);  new base T (6/4);  total (68/75)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

CGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCT  >  NC_002947/4151855‑4152149
                                                                                                                                                 |                                                                                                                                                     
cGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACC                                                                                                                                                   <  2:3341078‑M1/150‑1 (MQ=255)
 ggCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCC                                                                                                                                                  >  1:3689101‑M1/1‑150 (MQ=255)
   cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                                <  1:3323776‑M1/150‑1 (MQ=255)
   cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                                <  1:4402988‑M1/150‑1 (MQ=255)
   agcacggagagtgccagctggctgaaggccttggggtcgaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCcgc                                                                                                                                                >  2:746990‑M1/133‑150 (MQ=255)
   agcaaggagagtgccagctggctgaaggccttggggtagaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCcgc                                                                                                                                                >  1:1308621‑M1/133‑150 (MQ=255)
   agcaaggagagtgccagctggctgaaggccttggggtagaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCcgc                                                                                                                                                >  1:1351654‑M1/133‑150 (MQ=255)
   agcaaggagagggccagctggctgacggccttggggtagcggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCcgc                                                                                                                                                >  2:747165‑M1/133‑150 (MQ=255)
     tttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAACCCcccc                                                                                                                                              >  1:2451335‑M1/1‑146 (MQ=255)
     tttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCgcgc                                                                                                                                              >  2:2885374/1‑150 (MQ=255)
     tttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCgcgc                                                                                                                                              <  1:2893590/150‑1 (MQ=255)
       tCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCt                                                                                                                                            <  2:4394501/150‑1 (MQ=255)
       tCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCt                                                                                                                                            <  1:4511769/150‑1 (MQ=255)
        cAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaagg                                                                                                                                           <  1:4082548/150‑5 (MQ=255)
         gagagtgccagctggctgaaggccttggggtagaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCCGCGCCTGc                                                                                                                                          <  2:562416‑M1/24‑1 (MQ=255)
         aGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGc                                                                                                                                          >  2:1968866/1‑150 (MQ=255)
          gtaggtGTATGAGGAGCACCTGACCAGTATTCCCAATGTGAAACGGTTGAGCTCGAGGTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGGTCCCCTTGTAGCTGTTTCCTCTTGTGCGGGCCCCTTCGCGGGTTTAGCCGCGCCTggg                                                                                                                                         >  2:1908920/6‑148 (MQ=255)
          gAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAACCCGCGCCTGCa                                                                                                                                         >  2:146134/1‑150 (MQ=255)
             gCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                      <  2:4490179/150‑1 (MQ=255)
               tGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCTGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                    <  2:3811535/150‑1 (MQ=255)
                gTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTg                                                                                                                                   <  2:3172261/150‑1 (MQ=255)
                 tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGa                                                                                                                                  <  2:662065/150‑1 (MQ=255)
                 tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGa                                                                                                                                  <  1:756651/150‑1 (MQ=255)
                 tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACACGCGCCTGGAAATGTGa                                                                                                                                  >  1:729279/1‑150 (MQ=255)
                  aTGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGAc                                                                                                                                 <  2:2313508/150‑1 (MQ=255)
                  aTGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGAc                                                                                                                                 <  1:3167912/150‑1 (MQ=255)
                  aTGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGAc                                                                                                                                 <  1:229467/150‑1 (MQ=255)
                       gAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGc                                                                                                                            <  2:71195/150‑1 (MQ=255)
                         gCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGccc                                                                                                                          <  1:453076/150‑1 (MQ=255)
                            ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTc                                                                                                                       <  2:111387/150‑1 (MQ=255)
                                aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                   <  1:3648228/150‑1 (MQ=255)
                                 ccAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCGCCAGCAAAGGTGACGTAGCCCCTCGCAgg                                                                                                                  >  1:2545470/1‑150 (MQ=255)
                                   aGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAg                                                                                                                >  2:3952723‑M2/1‑149 (MQ=255)
                                    gTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGAAAACCCCCCCCGCCAAGGGGCCGTAGCCCCGCGCCGGAgc                                                                                                               <  1:2446959‑M2/150‑3 (MQ=255)
                                      aTTCCCAATGTGCAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCgcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:1297490‑M1/150‑36 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:2502708‑M1/150‑36 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:754363‑M1/150‑36 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:2365346‑M1/150‑36 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgg                                                                                                             >  1:2796354‑M2/1‑146 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgg                                                                                                             >  2:2501423‑M2/1‑146 (MQ=255)
                                      aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCGTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:793480‑M1/150‑36 (MQ=255)
                                       ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcggg                                                                                                            >  1:1571007‑M2/1‑145 (MQ=255)
                                       ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCGTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcggg                                                                                                            >  2:2460246‑M2/1‑145 (MQ=255)
                                       tagaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcccCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCggg                                                                                                            >  2:1220042‑M1/97‑150 (MQ=255)
                                        agaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgtgccggcccCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGt                                                                                                           >  2:509942‑M1/96‑150 (MQ=255)
                                        agaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGt                                                                                                           >  2:494000‑M1/96‑150 (MQ=255)
                                        agaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGt                                                                                                           >  2:1611497‑M1/96‑150 (MQ=255)
                                        agaggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcctCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGt                                                                                                           >  1:3936305‑M1/96‑150 (MQ=255)
                                          ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggttt                                                                                                         >  1:1600493‑M2/1‑142 (MQ=255)
                                          aggggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcttgtgccggcccCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGttt                                                                                                         >  2:3631863‑M1/94‑150 (MQ=255)
                                            aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttac                                                                                                       <  2:2155313‑M1/150‑42 (MQ=255)
                                             ggcggcgggggtggcaaccggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaaaggttgcctgggccggcccCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTAcc                                                                                                      >  2:2633231‑M1/91‑150 (MQ=255)
                                                   aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGGAAACCCCCCCCGCCAAGGGGTCCGTGCCCCGCGCAGGAGCCGGGTTTCCcgccca                                                                                                >  1:1165125/1‑147 (MQ=255)
                                                     gggtggcaacgggatcacggtgggctgggttggtgcgatgtgcgcgttttcagggcagatgaatggttgcctgggccggcctCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAaga                                                                                              >  2:468498‑M1/83‑150 (MQ=255)
                                                     gggtggcaaccgtatcagggtgggctgggttggtgcgatgtgcgcgttttctgggcagatgaaaggttgccggggcgggcctCTTGGCGGGTAAACCCGCGCCTGCAATGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAaga                                                                                              >  2:2723833‑M1/83‑150 (MQ=255)
                                                      gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATGGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAagag                                                                                             <  2:4275831/150‑1 (MQ=255)
                                                          gAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGACGTTTCCCGTCGCAGGAGCGGGTTTTACCCCGAAGAggccg                                                                                         >  2:2786286/1‑150 (MQ=255)
                                                               cggatcccggtgggctgggttggtgcgaggtgcgcgttttccgggcagatcaaagggtgcctgtgccggcccCTTCGCGGGTAAACCCGCGCCGGCAACGGTGACGTAGCCCGTCGCCGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCa                                                                                    >  2:1568338‑M1/73‑150 (MQ=255)
                                                                 aCCTGGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                  <  2:2959035/150‑1 (MQ=255)
                                                                  ccTTGGTGTTGACGTACGTGACCCTGGAGGGGTTGCTCCCCATGTAGCTGTTTCCTCTTGGGCGGGCCTCTGCGCGGGTAACCCCGCtcctgcgacgggctacgtcacctttgctggcgcgggtttacccgcgaaggggccggcacaaga                                                                                 <  2:629522‑M1/150‑64 (MQ=255)
                                                                   cTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCCCCAGCAAAGGTGACGTATCCCCTCGCAGGAGCGGGGTTAACCCCGAAGAGGCCCGCCCACGCa                                                                                >  1:651943/1‑150 (MQ=255)
                                                                     tGGTGATGACGGAAGTGATCCAGGATCGGTTGCACCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTCCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAAc                                                                              <  2:2415119/150‑1 (MQ=255)
                                                                        tgatgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGGCGTAGCCCCGCGCCGGAGCGGGGTTAACCCCCAAGAGGCCGGCCCACGCAAccat                                                                           >  1:2226331/1‑148 (MQ=255)
                                                                         gatgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCttt                                                                          <  2:1095147/150‑1 (MQ=255)
                                                                          atgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                         <  2:2243915/150‑1 (MQ=255)
                                                                            gaaggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             >  2:3573279‑M1/1‑77 (MQ=255)
                                                                            gaaggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:3573279‑M1/112‑36 (MQ=255)
                                                                            gaaggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGGAGCCCCGCGCAGGAGCGGGTTTACCCGCGAACAGGCCGGCCCAGGCAACCATTCaa                                                                       >  2:3172885/1‑149 (MQ=38)
                                                                                 aagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGcc                                                                  <  2:3365738/150‑1 (MQ=255)
                                                                                        cctggATGGGTTGCTCCCCTTGTCGCTGTTTCCTCTTGTGCCGGCCCCTGCGGGGGAAAACCCGCGCCTGCAATGGTGAGGTAGCCGGTCGCTGGTGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGaaaa                                                           <  1:1004708/147‑1 (MQ=255)
                                                                                        ccgggATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGGCGCAGGAGCGGGTTTACCCGCCAAGAGGCCGGCCCAGGCAACCATTAATCAGCCCAGaaaa                                                           >  2:1725119/4‑150 (MQ=255)
                                                                                           ggATCGGTTGCTCCCCTTTTTGGTTTTTCCTCTTGTGCCGGCCCCTTCGGGGGAAAACCCGCGCCTCAAATGGTGACGACGCCCGTCGCTGTTGGGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgc                                                        <  1:2155952/150‑1 (MQ=255)
                                                                                            gATCGGTGGCTCCCCATGTAGCTGTTACCTCTTGGGCCGGCCCCTTCGCGGGTAACCCCGCGCCTCCAAAGGTGACGTAGCCCGTCGCAGTAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgcg                                                       <  1:1810631/150‑1 (MQ=39)
                                                                                                      tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             >  2:2259373‑M1/1‑51 (MQ=255)
                                                                                                      tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  1:2259373‑M1/86‑36 (MQ=255)
                                                                                                      tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAATGTGACGGAGCCCGGCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCAc                                             >  2:303399/1‑150 (MQ=255)
                                                                                                      tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCGGAAAACGCGCACATCGCAc                                             >  2:303606/1‑150 (MQ=255)
                                                                                                      tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGGCGTAGCCCGTAGCAGGAGCGGGTTTTGCCGCGAAGGGGCCGGGCCAGGCAACCTTTCAACTGGCCAGAAAAAGCGCACAAAGCAc                                             >  2:3425681/1‑150 (MQ=255)
                                                                                                       ccccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCAGCACAGGGGCCCGAGCCCCTCCCAGGAGCCGGGTGACCCCCGAAGAGGCCCGCCCCCGCCACCAGTCATCAGCCCAGAAAACGCCCACAACGCAcc                                            >  2:3841819/1‑150 (MQ=255)
                                                                                                          cATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAAACCCCGCCCGCCAAGGGGAAGTAGCCCGGCGCAGGAGCGGGGGTAACCCCGAACAGGCCGGCCCAGGCCCCCGTTAATCAGCCCGGCAAACGCGCACAACGCCCCAAc                                         >  2:2636383/1‑150 (MQ=255)
                                                                                                             tttgCTGTTTCCTCTTGTGCCGGCCCCTTCGGGGTTAAACCCGCGCCTGCAATGGTGAGGTAGCCGGTCGCAGTTGGGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCa                                      <  1:3671433/147‑1 (MQ=255)
                                                                                                               aGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAACCCCGCGCCTGCAACGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGc                                    <  2:2451335/150‑1 (MQ=255)
                                                                                                                gCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGCCCGGCGCCGGGGCGGGGTTACCCGCGCAGAGGCCGGCCCAGGCCACCATTCATCTGCCCTGAAAAAGCGCCCCTAGCACCAACCCCGcc                                   >  2:1940659/1‑150 (MQ=255)
                                                                                                                    tttCCTCTTGTGCCGGCCCCTTCGCGGGTAAGCCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCAcc                               <  2:4027960/150‑1 (MQ=255)
                                                                                                                         tctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAAGGTGACGTAGCCCGTCGCAGGAgcggggttacccgcgcagaggcccgcccaagaggaaacagctacatggggagcaacccctccagggtcacgtacgtcaacaccaagg                          <  2:629522‑M2/150‑88 (MQ=255)
                                                                                                                         tctTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGTCGTAGCCCGTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCAGGCAACCGTTAAACTGCCCTGAAAAAGCGCACAAAGCACCAAACCAGCCCCCCCTGAt                          >  1:3006182/1‑150 (MQ=255)
                                                                                                                          ctTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCAGCAAAGGTGGCGTAGCCCGTCGCCGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCACGCCACCTTTCATATGCCCAGAAAAAGCGCACAACGCACCAACCCACcccccccggata                         >  2:2666412/1‑144 (MQ=25)
                                                                                                                              ggcCGGCCCCGGCGCGGGTCACCCCGCgccggcgccgggccacggcacctttgccggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgag                     <  1:3677393‑M1/149‑124 (MQ=255)
                                                                                                                               gCCGGCCCCTTCGCGGTTAAACCCGCGCCTGCAATGTTTACTTCCCCCGGCGCTGGAGGGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                    <  1:1968866/150‑1 (MQ=255)
                                                                                                                                          cggggttaaccccGCtcctgcgccgggctacgtcacctttgcaggcgggggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttc         <  2:3786199‑M1/144‑136 (MQ=255)
                                                                                                                                          cgcgGGTAAACCCGCtcctgcgacgggctacgtcacttttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttc         <  2:3479090‑M1/150‑136 (MQ=255)
                                                                                                                                              GGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc     <  2:1803505‑M1/150‑1 (MQ=255)
                                                                                                                                               gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttacccgcgaaggggccggcacaagaggacacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacatg    <  2:3721181‑M2/150‑110 (MQ=255)
                                                                                                                                               gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacatt    <  2:2155313‑M2/150‑110 (MQ=255)
                                                                                                                                               GTAAACCCGCGCCTGCAACGGTGACGTAGCCCGTCGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc    >  2:776802‑M1/1‑150 (MQ=255)
                                                                                                                                               GTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCCTCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc    >  2:776788‑M1/1‑150 (MQ=255)
                                                                                                                                               GTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc    <  2:3458960‑M1/150‑1 (MQ=255)
                                                                                                                                               GTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc    >  2:2336296‑M1/1‑150 (MQ=255)
                                                                                                                                                 aaaCCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgg  >  1:1600493‑M1/1‑8 (MQ=255)
                                                                                                                                                 aaaCCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgg  >  1:1600479‑M1/1‑8 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCT  >  NC_002947/4151855‑4152149

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: