Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,804,956 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,804,9560AC 5.4% 606.1 / 17.4 178K459Q (AAG→CAG) rbsA‑Iribose ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (74/94);  new base C (6/4);  total (80/98)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 212 total aligned reads displayed.

TGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAT  >  NC_002947/2804807‑2805105
                                                                                                                                                     |                                                                                                                                                     
tGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggca                                                                                                                                                       >  2:1134474/1‑150 (MQ=255)
tGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggca                                                                                                                                                       >  2:683068/1‑150 (MQ=255)
 gTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaa                                                                                                                                                      >  1:2223049/1‑150 (MQ=255)
   cTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaagg                                                                                                                                                    >  2:4376577/1‑150 (MQ=255)
    tGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCaaggcaaggc                                                                                                                                                   >  2:3695823/1‑150 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               >  2:2684285/1‑150 (MQ=255)
        ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtg                                                                                                                                               >  1:3228330/1‑150 (MQ=255)
           aaCCAGCAGACGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTAGGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggcg                                                                                                                                            <  1:1299248/150‑3 (MQ=255)
           aaCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtg                                                                                                                                            >  2:1730055/1‑150 (MQ=255)
           aaCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggcg                                                                                                                                            <  1:2346738/150‑3 (MQ=255)
            aCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggcgc                                                                                                                                           <  1:1605078/150‑4 (MQ=255)
             ccAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggt                                                                                                                                          <  2:722191/150‑1 (MQ=255)
             ccAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggt                                                                                                                                          <  1:2132257/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  2:2353436/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         <  2:1264492/150‑1 (MQ=255)
              cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                         >  1:1160241/1‑150 (MQ=255)
                gcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTggtggtggt                                                                                                                                       >  2:79289/1‑150 (MQ=255)
                 cagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGgtggtggtg                                                                                                                                      >  2:4064512/1‑150 (MQ=255)
                   gAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTc                                                                                                                                    <  2:1678799/150‑1 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  2:39226/150‑1 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:393757/150‑1 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:393572/150‑1 (MQ=255)
                     aGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                                  <  1:3020128/150‑1 (MQ=255)
                      ggtggtGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAg                                                                                                                                 >  2:3966185/1‑150 (MQ=255)
                       gtggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGt                                                                                                                                <  1:583203/150‑1 (MQ=255)
                        tggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTg                                                                                                                               >  1:937195/1‑150 (MQ=255)
                         ggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGa                                                                                                                              <  1:2412716/150‑1 (MQ=255)
                             aTCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTg                                                                                                                          >  1:79899/1‑150 (MQ=255)
                              tCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgc                                                                                                                         >  1:2019620/1‑150 (MQ=255)
                                ggCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       >  2:2376704/1‑150 (MQ=255)
                                ggCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                       <  1:864090/150‑1 (MQ=255)
                                   cGCTGGCTGGAGCGCGATTGCCTGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGaa                                                                                                                    <  2:3291692/150‑1 (MQ=255)
                                          tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                             <  2:1386896/150‑1 (MQ=255)
                                            gAGCGCTATTGCCAGGTGATGCTGTTCGACGAACCCACGCGAGCCATCGATGTGGGCACCAAGTTCGACATTTATGGCCTGCTGGACTAGCTGGCGCTCAAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTg                                                                                                           >  2:2302563/1‑150 (MQ=255)
                                            gAGCGCGATCGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTg                                                                                                           >  1:2934635/1‑150 (MQ=255)
                                              gcgcGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAt                                                                                                         <  2:2740662/150‑1 (MQ=255)
                                                gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                       <  2:527749/150‑1 (MQ=255)
                                                  gATTGCCAGGTGCTGCTGTTCGACGACCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGACCTGATGCTGATTTGc                                                                                                     >  2:1879918/1‑150 (MQ=255)
                                                    ttGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGa                                                                                                   <  2:168306/150‑1 (MQ=255)
                                                      gCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGAcc                                                                                                 <  1:1268867/150‑1 (MQ=255)
                                                      gCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGAcc                                                                                                 <  2:139572/150‑1 (MQ=255)
                                                       ccAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATTGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCg                                                                                                <  1:4038123/150‑1 (MQ=255)
                                                         aGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCa                                                                                              >  2:2723296/1‑150 (MQ=255)
                                                                gctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGt                                                                                       >  2:2690779/1‑150 (MQ=255)
                                                                 ctgTTCGACGCACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                      <  2:1793739/150‑1 (MQ=255)
                                                                 ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                      <  2:1408794/150‑1 (MQ=255)
                                                                 ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                      <  2:50858/150‑1 (MQ=255)
                                                                 ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                      <  2:3603150/150‑1 (MQ=255)
                                                                    ttCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGAc                                                                                                                             <  1:1035877/108‑1 (MQ=255)
                                                                    ttCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGAc                                                                                                                             >  2:1035877/1‑108 (MQ=255)
                                                                       gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                >  2:902260/1‑150 (MQ=255)
                                                                       gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                >  1:4481456/1‑150 (MQ=255)
                                                                       gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                >  1:2584173/1‑150 (MQ=255)
                                                                        acgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCg                                                                               <  2:1707642/150‑1 (MQ=255)
                                                                          gaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCt                                                                             >  2:2985749/1‑150 (MQ=255)
                                                                           aaCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTg                                                                            <  1:695541/150‑1 (MQ=255)
                                                                            aCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACGTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTgg                                                                           >  1:4011196/1‑150 (MQ=255)
                                                                            aCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTgg                                                                           >  2:3564923/1‑150 (MQ=255)
                                                                              ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgcc                                                                         >  2:2861639/1‑150 (MQ=255)
                                                                               cACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccg                                                                        >  1:3822679/1‑150 (MQ=255)
                                                                                  gcgTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGAc                                                                                                                             <  2:651543/94‑1 (MQ=255)
                                                                                    gTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGACCTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGa                                                                   >  2:3681641/1‑150 (MQ=255)
                                                                                           cGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGacac                                                            <  2:2681592/150‑1 (MQ=255)
                                                                                             aTGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCt                                                          >  2:2572282/1‑150 (MQ=255)
                                                                                                tgGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcg                                                       <  1:2350852/150‑1 (MQ=255)
                                                                                                 gggCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcgc                                                      <  1:2427831/150‑1 (MQ=255)
                                                                                                            tCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAtt                                           <  1:3478618/150‑1 (MQ=255)
                                                                                                             cGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTg                                          >  2:1430792/1‑150 (MQ=255)
                                                                                                             cGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTg                                          <  2:1674201/150‑1 (MQ=255)
                                                                                                               aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                        <  1:1199995/150‑1 (MQ=255)
                                                                                                               aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                        <  1:106900/150‑1 (MQ=255)
                                                                                                               aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                        <  1:1002249/150‑1 (MQ=255)
                                                                                                                 aTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGc                                      <  2:1945710/150‑1 (MQ=255)
                                                                                                                 aTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGc                                      <  2:3407681/150‑1 (MQ=255)
                                                                                                                    tATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAg                                   <  1:773367/150‑1 (MQ=255)
                                                                                                                     aTGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAgg                                  <  2:1222479/150‑1 (MQ=255)
                                                                                                                     aTGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAgg                                  >  2:3443068/1‑150 (MQ=255)
                                                                                                                          ctgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAg                             >  2:3053052/1‑150 (MQ=255)
                                                                                                                             ctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttg                          <  2:3259524/150‑1 (MQ=255)
                                                                                                                                  cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttgc                     >  2:470347/1‑150 (MQ=255)
                                                                                                                                  cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttgc                     >  1:1512045/1‑150 (MQ=255)
                                                                                                                                       tGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg                >  1:1446565/1‑150 (MQ=255)
                                                                                                                                       accaccaccAAGGCCTTGCCTTGGCGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg                >  1:3405902/8‑150 (MQ=255)
                                                                                                                                           gcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCtt            >  1:3062424/1‑150 (MQ=255)
                                                                                                                                           ccaccAAGGCCTTGCCTTGGCGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCtt            >  1:1511745/4‑150 (MQ=255)
                                                                                                                                               caaggcaaggCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGcc        >  1:3718572/1‑150 (MQ=255)
                                                                                                                                               caaggcaaggCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGcc        >  1:2502088/1‑150 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       >  1:3270581/1‑150 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  2:1594828/150‑1 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  2:3385491/150‑1 (MQ=255)
                                                                                                                                                aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg       <  1:2731145/150‑1 (MQ=255)
                                                                                                                                                 aggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCgg      >  1:2073510/1‑150 (MQ=255)
                                                                                                                                                 aggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCgg      <  1:1621224/150‑1 (MQ=255)
                                                                                                                                                  ggcacggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGt     >  1:3675421/1‑150 (MQ=255)
                                                                                                                                                   gcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGtt    >  1:4303220/1‑150 (MQ=255)
                                                                                                                                                   gcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGtt    <  1:2282005/150‑1 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  >  1:1267492/1‑150 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  >  1:1141154/1‑150 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  >  1:2959479/1‑150 (MQ=255)
                                                                                                                                                     aaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAt  >  1:3891171/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTAT  >  NC_002947/2804807‑2805105

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: