Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,807,983 A→C 5.4% K459Q (AAG→CAG) ‡ rbsA‑I → ribose ABC transporter ‑ ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,807,9830AC 5.4% 934.4 / 28.4 278K459Q (AAG→CAG) ‡rbsA‑Iribose ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (105/156);  new base C (9/7);  total (114/164)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAT  >  NC_002947/2807834‑2808132
                                                                                                                                                     |                                                                                                                                                     
 gTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaa                                                                                                                                                      <  1:6536960/150‑1 (MQ=255)
 gTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaa                                                                                                                                                      <  2:350371/150‑1 (MQ=255)
 gTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGAGATTGCCAGGTGCTGCTGCTCGACGAACCCACTCGTGGCATCGAGGTGGGCGCCAAGTTCCACATTTATGGCCTGCTGGCAGAGCTGGCGCgccaagccaa                                                                                                                                                      <  2:5617383/150‑1 (MQ=255)
  tCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaag                                                                                                                                                     <  2:4028727/150‑1 (MQ=255)
    tGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCaaggcaaggc                                                                                                                                                   >  1:1757770/1‑150 (MQ=255)
      gCGGTAACCAGCAGACGGGGGTGATCGGCCGCTGGCTGGAGCGCGATGGCCCGGTGCGGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACAGTTATGGCCTGCTGGCCGAGCTGGCGCGCCCAGGCCAGGCCt                                                                                                                                                 >  2:2981407/1‑150 (MQ=255)
      gCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCt                                                                                                                                                 >  1:4053580/1‑150 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               >  1:1351768/1‑150 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               <  2:6525782/150‑1 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               <  2:3296848/150‑1 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               >  1:4706760/1‑150 (MQ=255)
          tAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtggc                                                                                                                                             <  2:5503850/150‑2 (MQ=255)
          tAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggt                                                                                                                                             <  1:4983740/150‑1 (MQ=255)
            aCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtggcgc                                                                                                                                           <  1:4318265/150‑4 (MQ=255)
            aCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtgg                                                                                                                                           >  1:2801915/1‑150 (MQ=255)
             ccAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggt                                                                                                                                          <  2:3856448/150‑1 (MQ=255)
             ccAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggt                                                                                                                                          <  2:6214659/150‑1 (MQ=255)
              cagcagACGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGCACCCACGCGTGGCATCGATGTGGGCGCCTAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  2:5967276/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  1:4139303/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         >  1:3032317/1‑150 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  2:4558717/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  1:5483318/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         >  1:6729209/1‑150 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCGGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  2:5035504/150‑1 (MQ=255)
                gcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTggtggtggt                                                                                                                                       >  1:4614382/1‑150 (MQ=255)
                gcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTggtggtggt                                                                                                                                       >  1:3173202/1‑150 (MQ=255)
                gcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGATGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTggtggtggt                                                                                                                                       >  2:5585911/1‑150 (MQ=255)
                   gAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTc                                                                                                                                    >  1:729229/1‑150 (MQ=255)
                   gAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTc                                                                                                                                    >  1:5856635/1‑150 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:783769/150‑1 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:2452159/150‑1 (MQ=255)
                       gtggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGt                                                                                                                                <  2:5886115/150‑1 (MQ=255)
                        tggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTg                                                                                                                               >  2:3016777/1‑150 (MQ=255)
                        tggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTg                                                                                                                               <  2:6472741/150‑1 (MQ=255)
                          gtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGAc                                                                                                                             <  1:5015221/150‑1 (MQ=255)
                              tCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgc                                                                                                                         <  1:1660765/150‑1 (MQ=255)
                                ggCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACGTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       <  2:6008135/150‑1 (MQ=255)
                                ggCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       <  2:1745699/150‑1 (MQ=255)
                                 gCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGcgcg                                                                                                                      >  1:6401606/1‑150 (MQ=255)
                                 gCCGCTGGCTGGAGCGCGAGTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCATGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGcgcg                                                                                                                      >  1:972450/1‑150 (MQ=255)
                                   cGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGaa                                                                                                                    <  1:6804293/150‑1 (MQ=255)
                                        gctggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGAt                                                                                                               <  1:4550837/150‑1 (MQ=255)
                                          tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                             >  2:774863/1‑150 (MQ=255)
                                          tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                             >  1:6316354/1‑150 (MQ=255)
                                          tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                             <  1:4655102/150‑1 (MQ=255)
                                                gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                       <  1:4426611/150‑1 (MQ=255)
                                                gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                       >  1:939759/1‑150 (MQ=255)
                                                          ggTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCAt                                                                                             >  2:6599378/1‑150 (MQ=255)
                                                                gctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGt                                                                                       >  1:4258095/1‑150 (MQ=255)
                                                                   gTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCt                                                                                    <  1:3994643/150‑1 (MQ=255)
                                                                       gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                >  1:2510675/1‑150 (MQ=255)
                                                                       gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                >  1:864314/1‑150 (MQ=255)
                                                                          gaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCt                                                                             >  2:836466/1‑150 (MQ=255)
                                                                            aCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTgg                                                                           >  1:1972112/1‑150 (MQ=255)
                                                                             cccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgc                                                                          >  1:4914940/1‑150 (MQ=255)
                                                                              ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgcc                                                                         <  2:2195962/150‑1 (MQ=255)
                                                                                aCGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccgc                                                                       <  2:1455080/150‑1 (MQ=255)
                                                                                aCGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccgc                                                                       >  2:6070592/1‑150 (MQ=255)
                                                                                      ggCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:6509466/85‑1 (MQ=255)
                                                                                        cATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGa                                                               <  1:4045420/150‑1 (MQ=255)
                                                                                        cATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGa                                                               <  2:5771774/150‑1 (MQ=255)
                                                                                        cATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGa                                                               <  1:6029572/150‑1 (MQ=255)
                                                                                         aTCGATGTGGGCGCCAAGTTCGACATCTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGac                                                              >  2:5202667/1‑150 (MQ=255)
                                                                                           cGAGGTGGGCGCCAAGTTCGACTTTTATGTCCTGCTGGCCGAGCGGGCGCGCCAAGGAAAGGCATTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGAATCGCCGTGCTGTCCGCTGGCCGCCTTATCGacac                                                            <  1:2642022/150‑1 (MQ=255)
                                                                                               gtgGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGa                                                                                                   >  2:6370694/1‑107 (MQ=255)
                                                                                               gtgGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTc                                                        <  1:5180553/150‑1 (MQ=255)
                                                                                                   gcgcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcgcgc                                                    <  2:1661111/150‑1 (MQ=255)
                                                                                                   gcgcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcgcgc                                                    >  1:1156698/1‑150 (MQ=255)
                                                                                                       cAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGa                                                >  1:2270661/1‑150 (MQ=255)
                                                                                                       cAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGa                                                <  1:5630511/150‑1 (MQ=255)
                                                                                                        aaGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGAt                                               <  1:6350583/150‑1 (MQ=255)
                                                                                                          gTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCa                                             <  1:1520859/150‑1 (MQ=255)
                                                                                                           ttCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAt                                            >  2:2658272/1‑150 (MQ=255)
                                                                                                              gACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       <  1:4773320/72‑1 (MQ=255)
                                                                                                              gACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       >  2:4773320/1‑72 (MQ=255)
                                                                                                               aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                        <  1:1004875/150‑1 (MQ=255)
                                                                                                               aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                        <  2:2519426/150‑1 (MQ=255)
                                                                                                                 aTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGc                                      >  1:3346412/1‑150 (MQ=255)
                                                                                                                    tATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAg                                   <  2:6316354/150‑1 (MQ=255)
                                                                                                                          ctgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAg                             <  1:3852219/150‑1 (MQ=255)
                                                                                                                           tgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGt                            <  1:3456731/150‑1 (MQ=255)
                                                                                                                           tgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGt                            <  2:343967/150‑1 (MQ=255)
                                                                                                                                gCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAAttgctt                       >  2:1886273/1‑150 (MQ=255)
                                                                                                                                 gcGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCGGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttg                      <  1:2981407/149‑1 (MQ=255)
                                                                                                                                  cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttgc                     >  1:2713457/1‑150 (MQ=255)
                                                                                                                                  cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttgc                     >  2:1974229/1‑150 (MQ=255)
                                                                                                                                       tGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg                <  2:3653879/150‑1 (MQ=255)
                                                                                                                                       tGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg                >  2:543359/1‑150 (MQ=255)
                                                                                                                                       accaccaccAAGGCCTTGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg                >  1:2058924/8‑150 (MQ=255)
                                                                                                                                          cgcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCt             <  2:1757888/150‑1 (MQ=255)
                                                                                                                                            caccAAGGCCTTGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTc           >  2:1054735/3‑150 (MQ=255)
                                                                                                                                            caccAAGGCCTTGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTc           >  2:128151/3‑150 (MQ=255)
                                                                                                                                              ccAAGGCCTTGCCTTGGCGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGc         >  1:3765698/1‑150 (MQ=255)
                                                                                                                                               caaggcaaggCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGcc        <  2:2764783/150‑1 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  1:2834431/150‑1 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  2:4614398/150‑1 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  2:4334920/150‑1 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  >  1:2244608/1‑150 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  <  2:2560779/150‑1 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGACCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGCTCGACACCTTCGCGCGTGAACATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  <  2:3921330/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAT  >  NC_002947/2807834‑2808132

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