Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,804,957 A→T 5.8% K459M (AAG→ATG) ‡ 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,9570AT 5.8% 567.3 / 17.9 172K459M (AAG→ATG) ‡rbsA‑Iribose ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (81/81);  new base T (5/5);  total (86/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 188 total aligned reads displayed.

GTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTATC  >  NC_002947/2804808‑2805106
                                                                                                                                                     |                                                                                                                                                     
 tCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaag                                                                                                                                                      <  1:683895/150‑1 (MQ=255)
 tCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaag                                                                                                                                                      <  1:4464362/150‑1 (MQ=255)
   tGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCaaggcaaggc                                                                                                                                                    >  2:4751078/1‑150 (MQ=255)
      cGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCtt                                                                                                                                                 <  1:1660918/150‑1 (MQ=255)
        gTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtgg                                                                                                                                               <  1:3136204/150‑1 (MQ=255)
        gTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtgg                                                                                                                                               <  1:3136016/150‑1 (MQ=255)
        gTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtgg                                                                                                                                               <  2:1246599/150‑1 (MQ=255)
         tAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtggc                                                                                                                                              <  1:2034724/150‑2 (MQ=255)
         tAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCCTTGCCTtggc                                                                                                                                              <  1:2117599/150‑2 (MQ=255)
            ccAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGTGCTGGCGCGCCAAGGCAAGGCCTtggtggt                                                                                                                                           <  1:4294050/150‑1 (MQ=255)
             cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                          >  2:3264499/1‑150 (MQ=255)
             cagcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                          <  2:3893734/150‑1 (MQ=255)
              agcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtgg                                                                                                                                         <  2:3125414/150‑1 (MQ=255)
               gcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTggtggtggt                                                                                                                                        >  1:1597390/1‑150 (MQ=255)
                      gtggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGt                                                                                                                                 <  2:1330462/150‑1 (MQ=255)
                      gtggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGt                                                                                                                                 <  2:3910669/150‑1 (MQ=255)
                         gtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGAc                                                                                                                              >  2:4605823/1‑150 (MQ=255)
                           gATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCt                                                                                                                            <  1:920581/150‑1 (MQ=255)
                           gATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCt                                                                                                                            >  1:1504090/1‑150 (MQ=255)
                           gATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCt                                                                                                                            <  2:874520/150‑1 (MQ=255)
                            aTCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTg                                                                                                                           >  2:3191139/1‑150 (MQ=255)
                            aTCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTg                                                                                                                           >  1:561981/1‑150 (MQ=255)
                               ggCCGCTGGCTGGAGCGCGTTTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTgcgc                                                                                                                        <  1:2724260/150‑1 (MQ=255)
                                  cGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGaa                                                                                                                     <  2:2831955/150‑1 (MQ=255)
                                   gctggctgGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCACGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAAc                                                                                                                    <  2:3803266/150‑1 (MQ=255)
                                       gctggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGAt                                                                                                                >  2:1374393/1‑150 (MQ=255)
                                         tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                              >  2:2739105/1‑150 (MQ=255)
                                         tggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGc                                                                                                              <  1:1382475/150‑1 (MQ=255)
                                            aGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGa                                                                                                           <  1:2324190/150‑1 (MQ=255)
                                             gcgcGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAt                                                                                                          >  1:2720873/1‑150 (MQ=255)
                                               gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                        >  2:1723412/1‑150 (MQ=255)
                                               gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                        <  1:3264499/150‑1 (MQ=255)
                                               gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                        <  1:3255446/150‑1 (MQ=255)
                                               gcgATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                        <  2:4434864/150‑1 (MQ=255)
                                                 gATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGc                                                                                                      >  1:4605647/1‑150 (MQ=255)
                                                    tGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGAc                                                                                                   >  1:4136432/1‑150 (MQ=255)
                                                       cAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGc                                                                                                >  1:2246865/1‑150 (MQ=255)
                                                       cAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGCCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGc                                                                                                >  2:4795421/1‑150 (MQ=255)
                                                          gTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATc                                                                                             >  1:2774319/1‑150 (MQ=255)
                                                             ctgctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGcc                                                                                          <  2:4091822/150‑1 (MQ=255)
                                                              tgctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCg                                                                                         >  2:3188157/1‑150 (MQ=255)
                                                               gctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGt                                                                                        <  2:1416464/150‑1 (MQ=255)
                                                               gctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGt                                                                                        <  2:921456/150‑1 (MQ=255)
                                                                ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                       <  1:3211028/150‑1 (MQ=255)
                                                                ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTg                                                                                       <  1:1177267/150‑1 (MQ=255)
                                                                 tgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGc                                                                                      <  1:5049879/150‑1 (MQ=255)
                                                                  gTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGTTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCt                                                                                     >  2:3199013/1‑150 (MQ=255)
                                                                  gTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCt                                                                                     <  1:2585789/150‑1 (MQ=255)
                                                                   ttCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTg                                                                                    >  2:3374809/1‑150 (MQ=255)
                                                                      gacgaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTcc                                                                                 >  1:1096534/1‑150 (MQ=255)
                                                                         gaACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCt                                                                              >  2:430901/1‑150 (MQ=255)
                                                                             ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgcc                                                                          >  1:2714259/1‑150 (MQ=255)
                                                                               aCGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccgc                                                                        >  1:2434217/1‑150 (MQ=255)
                                                                                cgcgTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccgcc                                                                       >  2:4929840/1‑150 (MQ=255)
                                                                                    tGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGAt                                                                   >  2:2024480/1‑150 (MQ=255)
                                                                                        aTCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGac                                                               <  2:1193766/150‑1 (MQ=255)
                                                                                              gtgGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTc                                                         <  2:561981/150‑1 (MQ=255)
                                                                                               tgGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTTACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcg                                                        >  1:4006362/1‑150 (MQ=255)
                                                                                                gggCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcgc                                                       <  1:2699041/150‑1 (MQ=255)
                                                                                                    gcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGt                                                   >  2:1580963/1‑150 (MQ=255)
                                                                                                    gcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGt                                                   >  2:1718902/1‑150 (MQ=255)
                                                                                                      cAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGa                                                 >  1:4696752/1‑150 (MQ=255)
                                                                                                      cAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGa                                                 >  1:4696488/1‑150 (MQ=255)
                                                                                                      cAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGa                                                 >  1:2052253/1‑150 (MQ=255)
                                                                                                         gTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCa                                              <  2:4136432/150‑1 (MQ=255)
                                                                                                          ttCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAt                                             <  1:3068180/150‑1 (MQ=255)
                                                                                                           tCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAtt                                            >  1:383457/1‑150 (MQ=255)
                                                                                                           tCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAtt                                            <  2:480721/150‑1 (MQ=255)
                                                                                                           tCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCAtt                                            <  2:513112/150‑1 (MQ=255)
                                                                                                            cGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTg                                           <  1:1374519/150‑1 (MQ=255)
                                                                                                              aCATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGa                                         <  1:3437679/150‑1 (MQ=255)
                                                                                                                aTTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGc                                       >  1:3377481/1‑150 (MQ=255)
                                                                                                                  ttATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCa                                     >  1:3945820/1‑150 (MQ=255)
                                                                                                                    aTGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                        <  2:137947/81‑1 (MQ=255)
                                                                                                                    aTGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAgg                                   <  1:2738221/150‑1 (MQ=255)
                                                                                                                         tcgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAg                              >  2:2645236/3‑150 (MQ=255)
                                                                                                                         ctgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAg                              >  2:3699945/1‑150 (MQ=255)
                                                                                                                          tgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGt                             <  2:518033/150‑1 (MQ=255)
                                                                                                                          tgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGt                             <  2:3470779/150‑1 (MQ=255)
                                                                                                                           gctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGtt                            <  2:1086036/150‑1 (MQ=255)
                                                                                                                             tgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgc                          >  2:3174847/1‑150 (MQ=255)
                                                                                                                             tgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgc                          <  1:789529/150‑1 (MQ=255)
                                                                                                                             tgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgc                          >  2:240638/1‑150 (MQ=255)
                                                                                                                               gCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgctt                        <  1:3004582/150‑1 (MQ=255)
                                                                                                                               gCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgctt                        <  2:3998064/150‑1 (MQ=255)
                                                                                                                               gCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgctt                        <  2:3041494/150‑1 (MQ=255)
                                                                                                                                 cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttgcttgc                      >  1:1078820/1‑150 (MQ=255)
                                                                                                                                   aGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccg                    <  2:618452/150‑1 (MQ=255)
                                                                                                                                        gcgcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccgcc               >  2:3407633/1‑150 (MQ=255)
                                                                                                                                        gcgcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccgcc               >  1:4753227/1‑150 (MQ=255)
                                                                                                                                        gcgcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccgcc               >  2:1803588/1‑150 (MQ=255)
                                                                                                                                        caccaccAAGGCCTTGCCTTGGCGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTCCTGTCCGCTGGCCGCCAGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGATGCTTgccgccgcc               >  1:2198642/6‑150 (MQ=255)
                                                                                                                                           caccAAGGCCTTGCCTTGGCGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTc            >  2:4887907/3‑150 (MQ=255)
                                                                                                                                            gcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGACGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCg           >  2:4777203/1‑150 (MQ=255)
                                                                                                                                              caaggcaaggCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGcc         >  2:2249683/1‑150 (MQ=255)
                                                                                                                                              caaggcaaggCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGcc         >  1:4100326/1‑150 (MQ=255)
                                                                                                                                               aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg        <  1:4789468/150‑1 (MQ=255)
                                                                                                                                               aaggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg        >  2:2806190/1‑150 (MQ=255)
                                                                                                                                                  gcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGtt     >  1:4735379/1‑150 (MQ=255)
                                                                                                                                                     aggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTATc  <  2:1504090/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGTTATC  >  NC_002947/2804808‑2805106

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: