Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,164,150 G→C 6.4% intergenic (+105/‑104) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,164,1500GC 6.4% 738.7 / 31.1 219intergenic (+105/‑104)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (101/104);  new base C (11/3);  total (112/107)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 237 total aligned reads displayed.

GGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCAAAGGCA  >  NC_002947/4164001‑4164299
                                                                                                                                                     |                                                                                                                                                     
ggggggTTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcagg                                                                                                                                                       <  2:6102979/150‑1 (MQ=255)
 gggggCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcaggc                                                                                                                                                      <  1:4491984/150‑1 (MQ=255)
  ggggCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcaggca                                                                                                                                                     <  2:7135747/150‑1 (MQ=255)
          cGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCGTACCGGGGGAAACCCCGCCGACAGGTGGGGGGGGCCCGGAAGAAGGCAACCCCCCCGGAGGCGCGGGGTCAACCCCCACCCAGGCCGCACgggc                                                                                                                                             >  2:6854914/1‑148 (MQ=255)
            aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtctc                                                                                                                                           <  1:5019116/150‑1 (MQ=255)
             ttACCCTGGGGTTGCAGATGCTGTTGGGCTGACTGGAGACCGAGTGGGCTGATTGGGGGGTAAACCCGCTCCTACTGGTGGGTGCTGCCCTGATGACGGCACGGCCCCTTTTGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                                                                                                          <  1:2291574/150‑1 (MQ=255)
               aCCCGGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAg                                                                                                                                        <  1:1682686/150‑1 (MQ=255)
                 ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAAACCGCGCCAACAGGTGGGGGGTGCCCGGGAGACGGCCAGGCCCCCGGAGGAGCGGGGTGACCCCCGAAGCCGGCAGCGCGGGCTCAAGGt                                                                                                                                      >  2:7294677/1‑150 (MQ=255)
                  gggggggTGCAGATGCTGTTGGGCTGACTGGCGACCGAGTCGGCTGCTTCGGGGGTAACCCCGCTCCCCCCGGGGGGTGCTGCCCTGAAGCCGGCACGGCCCCTGTTGGAGCGGGTTCACCCGCGACGCAGGCAGCTCGGTCTCAAGGtt                                                                                                                                     <  1:6301024/148‑1 (MQ=255)
                   tGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttt                                                                                                                                    <  2:6331064/150‑1 (MQ=255)
                                          tGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGCAGACGGCAAGGCCCCTGGAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGt                                                                                                             <  2:6513053/150‑1 (MQ=255)
                                           gACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTc                                                                                                            <  2:5419604/150‑1 (MQ=255)
                                           gACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTc                                                                                                            <  2:2176566/150‑1 (MQ=255)
                                                       aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCgg                                                                                                <  2:4496816/150‑1 (MQ=255)
                                                           ggCTGCTTCGCGGGTAAACCCGCTCCTACTGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtg                                                                                            <  1:1850954/150‑1 (MQ=255)
                                                                ctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt                                                                                       <  2:7286436/150‑1 (MQ=255)
                                                                  gagcgggttTACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                                                                                                 >  1:5329136‑M1/10‑138 (MQ=255)
                                                                   agcgggttcacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTc                                                                                    >  1:3067400‑M2/69‑150 (MQ=255)
                                                                     cgGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaa                                                                                  <  1:493544/150‑1 (MQ=255)
                                                                      gggTAACCCCGCTCCTACCGGGGGGTGCGGCCCTGAAGACGGCCAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaaa                                                                                 <  1:6141944/150‑1 (MQ=255)
                                                                      gggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaaa                                                                                 >  1:6705222/1‑150 (MQ=255)
                                                                        gTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaaga                                                                               <  2:6418076/150‑1 (MQ=255)
                                                                         tAAACCGGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCTAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTcttagta                                                                              <  2:1513639/150‑7 (MQ=255)
                                                                         tAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                                                                              <  1:6782544/150‑1 (MQ=255)
                                                                         tAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                                                                              <  2:1988343/150‑1 (MQ=255)
                                                                          gaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaag                                                                             >  2:6255335‑M2/62‑150 (MQ=255)
                                                                          gaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaag                                                                             >  2:676156‑M2/62‑150 (MQ=255)
                                                                            ACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGa                                                                           >  1:4050481‑M1/1‑150 (MQ=255)
                                                                             cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt                                                                          >  2:2974489‑M2/59‑150 (MQ=255)
                                                                             cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt                                                                          >  1:5854324‑M2/59‑150 (MQ=255)
                                                                             cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt                                                                          >  2:4901716‑M2/59‑150 (MQ=255)
                                                                             CCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt                                                                          >  1:1220767‑M1/1‑150 (MQ=255)
                                                                              ccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         >  1:5287180‑M2/58‑150 (MQ=255)
                                                                              ccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         >  1:1549118‑M2/58‑150 (MQ=255)
                                                                              CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         <  2:2042287‑M1/150‑1 (MQ=255)
                                                                              CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         <  1:2094631‑M1/150‑1 (MQ=255)
                                                                              CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         <  2:1815800‑M1/150‑1 (MQ=255)
                                                                              CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                                         <  2:532108‑M1/150‑1 (MQ=255)
                                                                              CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGATGCCGTCCACATATCTCGGTATGGTGGGGTCAAAGAAGGATg                                                                         <  2:5572145‑M1/150‑1 (MQ=255)
                                                                                 ctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGttt                                                                      >  1:6853302‑M2/55‑150 (MQ=255)
                                                                                 CTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGttt                                                                      >  1:3717823‑M1/1‑150 (MQ=255)
                                                                                  TCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGtttt                                                                     <  1:2848903‑M1/150‑1 (MQ=255)
                                                                                   CCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTa                                                                    <  1:2728918‑M1/150‑1 (MQ=255)
                                                                                    ctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTAc                                                                   >  1:5602305‑M2/52‑150 (MQ=255)
                                                                                    CTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTAc                                                                   <  2:4421555‑M1/150‑1 (MQ=255)
                                                                                     tacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTAcc                                                                  >  1:1362278‑M2/51‑150 (MQ=255)
                                                                                     TACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGAt                                                                          >  1:846273‑M1/1‑142 (MQ=255)
                                                                                     TACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTAcc                                                                  <  2:5852207‑M1/150‑1 (MQ=255)
                                                                                       CAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCg                                                                <  1:5447688‑M1/150‑1 (MQ=255)
                                                                                       CAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCg                                                                <  2:4203614‑M1/150‑1 (MQ=255)
                                                                                       CAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCg                                                                >  1:4566735‑M1/1‑150 (MQ=255)
                                                                                        AGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcagg                                                                                                                                                       >  2:3709955‑M1/1‑62 (MQ=255)
                                                                                          gggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac                                                             >  2:1361394‑M2/46‑150 (MQ=255)
                                                                                          GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaaa                                                                                 >  1:5736099‑M1/1‑130 (MQ=255)
                                                                                          GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCaaa                                                                                 <  2:5736099‑M1/130‑1 (MQ=255)
                                                                                          GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAGAGAAGGATGTTTTACCCGcac                                                             <  1:2234994‑M1/150‑1 (MQ=255)
                                                                                          GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac                                                             <  2:3193877‑M1/150‑1 (MQ=255)
                                                                                          GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac                                                             <  1:2887849‑M1/150‑1 (MQ=255)
                                                                                           tgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcaca                                                            <  2:5162978/150‑1 (MQ=255)
                                                                                             ggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAgg                                                          <  2:936612/150‑1 (MQ=255)
                                                                                               ttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGaa                                                        >  2:3454179‑M2/41‑150 (MQ=255)
                                                                                                 ctgcCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAg                                                      >  2:182974/1‑150 (MQ=255)
                                                                                                  tgcCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGa                                                     >  1:1517005/1‑150 (MQ=255)
                                                                                                    cccTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcaca                                                                                                                              >  1:6912106/1‑75 (MQ=255)
                                                                                                    cccTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACt                                                   >  2:4227517/1‑150 (MQ=255)
                                                                                                      cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTa                                                 >  1:1612338/1‑150 (MQ=255)
                                                                                                      cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTa                                                 >  2:488060/1‑150 (MQ=255)
                                                                                                      cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTa                                                 >  2:3763748/1‑150 (MQ=255)
                                                                                                      cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTATTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTa                                                 <  2:2418966/150‑1 (MQ=255)
                                                                                                           gACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGAGGc                                            <  2:777365/150‑1 (MQ=255)
                                                                                                              ggCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt                                                                                       >  1:1050701/1‑104 (MQ=255)
                                                                                                              ggCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt                                                                                       <  2:1050701/104‑1 (MQ=255)
                                                                                                                 aaGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAAc                                      >  1:6347315/1‑150 (MQ=255)
                                                                                                                   ggCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGc                                    <  1:2650689/150‑1 (MQ=255)
                                                                                                                    gCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGcc                                   >  2:1111252/1‑150 (MQ=255)
                                                                                                                    ccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcaca                                                                                                                              <  2:6553811‑M2/40‑1 (MQ=255)
                                                                                                                    ccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcaca                                                                                                                              >  1:6553811‑M2/20‑59 (MQ=255)
                                                                                                                     ccccTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCa                                  >  2:3406809/1‑150 (MQ=255)
                                                                                                                     ccccTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCa                                  <  1:637275/150‑1 (MQ=255)
                                                                                                                      cccTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt                                 >  1:4014125/1‑150 (MQ=255)
                                                                                                                        ctgttggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGt                               <  1:1745294‑M2/135‑1 (MQ=255)
                                                                                                                        ctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGt                               >  1:185685‑M2/16‑150 (MQ=255)
                                                                                                                         tgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGtc                              <  2:5476812‑M2/136‑1 (MQ=255)
                                                                                                                           taggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGtcat                            <  1:3702610‑M2/138‑1 (MQ=255)
                                                                                                                             ggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGtcatca                          >  1:3821120‑M2/11‑150 (MQ=255)
                                                                                                                              gagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCaa                         >  1:7017393‑M2/10‑150 (MQ=255)
                                                                                                                                gcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATg                       >  2:39330‑M2/8‑150 (MQ=255)
                                                                                                                                   ggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGt                               >  1:7332525‑M2/5‑139 (MQ=255)
                                                                                                                                    gttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATc                   >  2:1853454‑M2/4‑150 (MQ=255)
                                                                                                                                     ttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATcc                  <  2:2906651‑M2/148‑1 (MQ=255)
                                                                                                                                        aCCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTc               >  1:4538822/1‑150 (MQ=255)
                                                                                                                                        aCCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTc               >  2:1616814/1‑150 (MQ=255)
                                                                                                                                        aCCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTc               >  2:4292296/1‑150 (MQ=255)
                                                                                                                                         cccGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCa              >  2:5376005/1‑150 (MQ=255)
                                                                                                                                            gcgAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGcc           >  1:4146101/1‑150 (MQ=255)
                                                                                                                                                 gcaggcagCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCaaa      <  2:697401/150‑1 (MQ=255)
                                                                                                                                                 gcaggcagCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCaaa      <  1:6498812/150‑1 (MQ=255)
                                                                                                                                                  caggcagCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCAAAg     <  1:6711512/150‑1 (MQ=255)
                                                                                                                                                  caggcagCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCAAAg     >  2:2328530/1‑150 (MQ=255)
                                                                                                                                                     gcagCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCAAAGGCa  >  1:1436712/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAATGCATCCGTCAGCCCCAAAGGCA  >  NC_002947/4164001‑4164299

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: