Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,908,812 T→G 8.9% intergenic (+63/+44) PP_5176 → / ← spuH hypothetical protein/spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,908,8120TG 8.9% 792.6 / 4.8 302intergenic (+63/+44)PP_5176/spuHhypothetical protein/spermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (143/118);  new base G (35/3);  total (181/121)
Fisher's exact test for biased strand distribution p-value = 4.36e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

GCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5908663‑5908961
                                                                                                                                                     |                                                                                                                                                      
gCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCt                                                                                                                                                        <  1:2285516/150‑1 (MQ=255)
 cccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCtg                                                                                                                                                       >  2:3822062/1‑150 (MQ=255)
  ccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgt                                                                                                                                                      <  2:6557854/150‑1 (MQ=255)
  ccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgt                                                                                                                                                      <  1:1919147/150‑1 (MQ=255)
    ggTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGa                                                                                                                                                   <  1:2843610/150‑1 (MQ=255)
    cgattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGACCTCACAACCCTGAAGGCTGTGAg                                                                                                                                                  >  1:751083‑M1/110‑150 (MQ=255)
       aCATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGt                                                                                                                                                <  2:1526375/150‑1 (MQ=255)
          ttgttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccc                                                                                                                                             <  1:3791870/150‑1 (MQ=255)
          ttgttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACAGGGCGCACCACAACCCAGAAGGCGGTGAGGTccc                                                                                                                                             >  1:6558788/1‑150 (MQ=255)
            gttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtg                                                                                                                                           <  1:606606/150‑1 (MQ=255)
             ttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          >  2:3881502/1‑150 (MQ=255)
             ttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCTGAAGGCTGCGAGGTCCCggt                                                                                                                                          >  1:3352144/1‑150 (MQ=255)
             aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          >  2:3174967‑M1/101‑150 (MQ=255)
              ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         >  1:4437408‑M1/100‑150 (MQ=255)
              ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGACCCCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  2:575569‑M1/100‑150 (MQ=255)
               gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  2:5546219‑M1/99‑150 (MQ=255)
               gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  1:3361578/150‑1 (MQ=255)
               gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  1:2126650/150‑1 (MQ=255)
               gCCAATCCATTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  2:5011156/150‑1 (MQ=255)
                 gatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGa                                                                                                                                      >  1:3633903‑M1/97‑150 (MQ=255)
                 gatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGa                                                                                                                                      >  1:6867952‑M1/97‑150 (MQ=255)
                 cAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCACCCACAGGGCTCACCACAACCCTGAAGGCTTTGAGGTCGCGGTgggg                                                                                                                                      >  1:4926013/1‑149 (MQ=255)
                   atcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggcccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGc                                                                                                                                    >  1:5965112‑M1/95‑150 (MQ=255)
                     gaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCgg                                                                                                                                  >  1:4708020‑M1/93‑150 (MQ=255)
                     cccTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCgg                                                                                                                                  <  1:1312893/150‑1 (MQ=255)
                      ccTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAGCCTTCAGGGTTGTGGTGAGCCCTGTGGGAGCgg                                                                                                                                  <  2:3027801/150‑1 (MQ=25)
                      ccTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATGGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGATCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCggg                                                                                                                                 <  2:1893265/150‑1 (MQ=255)
                       agagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGc                                                                                                                                >  1:1837519‑M1/91‑150 (MQ=255)
                        ttGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAAGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGTGGTCCCGGTGGGAGCGGGCg                                                                                                                               >  2:1292224/1‑150 (MQ=255)
                        ttGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCCGGCCTTCCCCCAACCACACGGCTCAACACAACCCAAGAAGCTTTGTGGTCCCCGTGGGAGCGGGCg                                                                                                                               >  1:4099940/1‑150 (MQ=255)
                           aGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGc                                                                                                                            >  1:4975603/1‑150 (MQ=255)
                             tgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgccgggggctgcgtcggcctgacctggccccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGccc                                                                                                                          >  1:3192936‑M1/85‑150 (MQ=255)
                                    ttGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCATCCCGCGCCCCCACGGCTCACCACAACCCAGAAGGCGGTGAGGTCCCGGGGGGTGCGGGCGTGCCCGCCAAAg                                                                                                                   >  1:4654894/1‑150 (MQ=255)
                                     tGGTCAGATGTATCAGGGACTATCCACATGGGGATGCGGTATGGCTTTGATTCGGGGCTGGCGCCTTCGCGGGCTTGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAgg                                                                                                                  <  2:1137320/150‑1 (MQ=38)
                                       gTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGc                                                                                                                <  1:5716594/150‑1 (MQ=255)
                                       gTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCGGTGGGCGCGGGTGTGCCCGCGAAAGGGc                                                                                                                >  2:3770786/1‑150 (MQ=255)
                                          aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                             >  2:6708705/1‑150 (MQ=255)
                                          aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCGCACCACAACCCGGAAGGCTGGGGGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                             >  1:6219399/1‑150 (MQ=39)
                                          aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCGGGAGCCTTCGCGGGCAGGCCCGCTCCCCCCGGGCGCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                             <  2:227588/150‑1 (MQ=39)
                                           gAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGCAAGGGCCAgg                                                                                                            >  2:6859303/1‑150 (MQ=255)
                                            aaGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGt                                                                                                           >  1:2188100/1‑150 (MQ=255)
                                              gTATCTGGCACTACCCACTTGGGTTTGGGTTAGGGCTTTGAATCGGTGCTGGCGCCTTCGCGGGCTTGCCCGCTCCCCCAGGGCTCACCACAACCCTGAAGGCTGTGTGGTCCCTGGGGGTGCGGGCGTGCCCGCGAAAGGGCCAGGTCa                                                                                                         <  1:5308001/150‑1 (MQ=255)
                                               tATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAAACCAGAAGGCGGTGAGGTCCCGGTGGGGGCGGGCGTGCCCGCGAAAGGGCCAgggaag                                                                                                        >  1:6396044/1‑146 (MQ=39)
                                                 tCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCAGCCCCCGCCCACCCGGCGCAACACAACCCAGAAGGCGGTGAGGTCCCCGGGGGAGCGGGCGTGCCCCCCAAAAGGCCCGGGCAGGc                                                                                                      >  1:2382139/1‑150 (MQ=255)
                                                     ggACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAATCCCGCGCCCCCACGGCTCAACACAAACCAGAAGGCGGGGTGGTGCCGGTGGGCGCGGGCGTGCCCCCGCAAGGGGCACGGCACGCCGAc                                                                                                  >  1:1167351/1‑150 (MQ=255)
                                                        cTATCCACATTGGGTTGCGTTTTGGCTTTGAATCGGTGCTGGAGCCTGCGCGGGCAGGCCCGCTCCCCCCGGGCACACCACAACCCTGAAGCCTGTGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCa                                                                                               <  1:1365300/150‑1 (MQ=38)
                                                            ccACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGGTCAGGCCGACGCAGCCa                                                                                           >  2:4882000/1‑150 (MQ=255)
                                                            ccACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGccc                                                                                           >  2:3557036/1‑149 (MQ=255)
                                                             cacaTTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCGCGCCCACACGGCTCAACACAAACCAGAAGGCGGTGAGGTCCCGCTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCAGGCCCACGCAGCCAc                                                                                          >  1:2685525/1‑150 (MQ=255)
                                                                gcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         >  2:5382983‑M1/50‑100 (MQ=255)
                                                                gcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         <  1:5382983‑M1/51‑1 (MQ=255)
                                                                     gATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCg                                                                                  >  2:3356853/1‑150 (MQ=255)
                                                                           gTAGGGCTTTGACTCGGTGCGGGCGCCTTCGCGGGCAGGCCCGCTCCCCCAGGGCTCACCACACCCCTGAAGGCGGTGAGGGCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCt                                                                            <  1:6745075/150‑1 (MQ=38)
                                                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCGCTCCCACAGGGCTCAACACAAACCCGAAGGCGGTGCGGTCCCTGGGGGCGCGGGCGGGCCCGCGGAAGGGCCAGGGCAGGCCGACGCAGCCACCGGCGCCGGGCGCTTGAcca                                                                      >  2:5807477/1‑149 (MQ=255)
                                                                                   ttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGc                                                                    <  1:5600993/150‑1 (MQ=255)
                                                                                     gAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCgg                                                                  <  2:4323838/150‑1 (MQ=255)
                                                                                     gAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCgg                                                                  >  2:2732055/1‑150 (MQ=255)
                                                                                       aTCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCg                                                                                                                                   >  1:2487623/1‑83 (MQ=255)
                                                                                         cGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACACCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGc                                                              <  2:2161657/150‑1 (MQ=255)
                                                                                         cGGGGCTGGAGCCTTCGCGGGCAGGCCCGCCCCCCCCGGGCACACCACACCCCTGAAGGCTGGGAGGTCCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGc                                                              <  1:390437/150‑1 (MQ=255)
                                                                                            tGCGGGAGCCTTCGCGGGCAGGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                           <  1:5256049/150‑1 (MQ=255)
                                                                                                ggAGCCTTCGCGGGCATGCCCGCGCCCACACGGCTCACCACAACCCAGAAGGCGGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCGGCGCCCGGCGCGTGACCAGCGGCTGCGAAGCGt                                                       >  2:3699604/1‑150 (MQ=38)
                                                                                                    ccTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGc                                                   >  1:432057/1‑150 (MQ=255)
                                                                                                     cTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCt                                                  >  2:2550830/1‑150 (MQ=255)
                                                                                                      ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTGGGCTg                                                 >  1:5181604/1‑150 (MQ=255)
                                                                                                       tCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgc                                                <  1:1034681/150‑1 (MQ=255)
                                                                                                       tCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgc                                                <  1:6266316/150‑1 (MQ=255)
                                                                                                            ggCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCAc                                           >  2:3297945/1‑150 (MQ=255)
                                                                                                              cATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTgcgc                                                                                    <  1:2770296/107‑1 (MQ=255)
                                                                                                              cATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCAccc                                         >  1:6337927/1‑150 (MQ=255)
                                                                                                               aTGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCct                                        >  2:4172327/1‑150 (MQ=255)
                                                                                                                tGCCCGCGCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctc                                       <  2:4625521/150‑1 (MQ=255)
                                                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                    <  2:2773386‑M1/150‑1 (MQ=255)
                                                                                                                     GCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGa                                  >  2:4740462‑M1/1‑150 (MQ=255)
                                                                                                                       TCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                <  1:3440621‑M1/150‑1 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTg                                                     >  2:774892‑M1/1‑126 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             <  2:2905179‑M1/150‑1 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             <  1:1628981‑M1/150‑1 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGACCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             <  2:2472569‑M1/150‑1 (MQ=255)
                                                                                                                             AGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                          <  2:810635‑M1/150‑1 (MQ=255)
                                                                                                                             AGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                          >  2:4661599‑M1/1‑150 (MQ=255)
                                                                                                                             AGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                          <  2:529143‑M1/150‑1 (MQ=255)
                                                                                                                                GCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGa                       >  2:2842294‑M1/1‑150 (MQ=255)
                                                                                                                                   caccaccacCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATc                    <  2:3907714/150‑1 (MQ=255)
                                                                                                                                   caccacAACCCGGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATc                    <  2:4926013/150‑1 (MQ=255)
                                                                                                                                     ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                  >  2:5170199/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  2:4841717/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGCGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:6629592/1‑150 (MQ=255)
                                                                                                                                       acaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCt                >  1:5759085/1‑150 (MQ=255)
                                                                                                                                       acaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCt                >  2:6071252/1‑150 (MQ=255)
                                                                                                                                       acaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCt                <  1:3248282/150‑1 (MQ=255)
                                                                                                                                           cccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcg            >  2:1029422/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  2:4374616/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  1:484549/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  1:4338860/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           <  1:1515482/150‑1 (MQ=255)
                                                                                                                                             cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac          >  1:1089169/1‑150 (MQ=255)
                                                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg         >  2:6362829/1‑150 (MQ=255)
                                                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg         >  1:3516324/1‑150 (MQ=255)
                                                                                                                                                     tgtgAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGctct  <  1:3822062/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5908663‑5908961

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