Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,585,620 A→G 5.4% I255V (ATC→GTC) ‡ benR → BenABC operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,585,6200AG 5.4% 1268.2 / 35.5 390I255V (ATC→GTC) ‡benRBenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (200/169);  new base G (8/13);  total (208/182)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 391 total aligned reads displayed.

CCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3585475‑3585769
                                                                                                                                                 |                                                                                                                                                     
ccTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCAc                                                                                                                                                   <  2:7541612/150‑1 (MQ=255)
 cTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACg                                                                                                                                                  >  1:5497715/1‑150 (MQ=255)
 cTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACg                                                                                                                                                  >  1:6111756/1‑150 (MQ=255)
  tGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGc                                                                                                                                                 >  1:10024/1‑150 (MQ=255)
  tGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGc                                                                                                                                                 >  1:1268077/1‑150 (MQ=255)
      tGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCaa                                                                                                                                             <  2:3532422/150‑1 (MQ=255)
      tGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCaa                                                                                                                                             <  1:4836717/150‑1 (MQ=255)
        ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgt                                                                                                                                           <  1:3970202/150‑2 (MQ=255)
         cccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCa                                                                                                                                          >  2:6141231/1‑150 (MQ=255)
          ccGCAGGCNGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCAt                                                                                                                                         >  1:4576146/1‑150 (MQ=255)
           cGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCAtt                                                                                                                                        >  1:7490469/1‑150 (MQ=255)
            gCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGACATTGGAATTGTCGGCTGAAGTGCGGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCAAgtgctg                                                                                                                                       <  2:5614722/150‑6 (MQ=255)
            gCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCAAGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgctg                                                                                                                                       <  1:7159906/150‑6 (MQ=255)
             cAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTAc                                                                                                                                      <  2:2267630/150‑1 (MQ=255)
             cAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTGTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTAc                                                                                                                                      <  1:5998858/150‑1 (MQ=255)
                 cTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                                  <  2:2957557/150‑1 (MQ=255)
                   gggCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgcgc                                                                                                                                >  1:5017133/1‑150 (MQ=255)
                    ggCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGcc                                                                                                                               >  2:6434441/1‑150 (MQ=255)
                     gCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCa                                                                                                                              >  1:2896728/1‑150 (MQ=255)
                       tAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                            >  2:742674/1‑150 (MQ=255)
                            cgcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                       >  2:1033376/1‑150 (MQ=255)
                                 tCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGa                                                                                                                  >  1:6473407/1‑150 (MQ=255)
                                  cTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                                 <  1:5012842/150‑1 (MQ=255)
                                    tGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACg                                                                                                               <  2:5789088/150‑1 (MQ=255)
                                     gATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                              <  1:916907/150‑1 (MQ=255)
                                       ttACATCGAACGCAACCTGAAATTGGAATAGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGt                                                                                                            >  1:8009189/1‑150 (MQ=255)
                                          cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTAGGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                         >  1:7655147/1‑150 (MQ=255)
                                           aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAt                                                                                                        <  2:656321/150‑1 (MQ=255)
                                           aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCCCAAACTGGAACGGGTGCAt                                                                                                        >  1:1117589/1‑150 (MQ=255)
                                              gAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcc                                                                                                     >  2:4232757/1‑150 (MQ=255)
                                               aaCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcct                                                                                                    <  1:539486/150‑1 (MQ=255)
                                                 cGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                                  <  1:528522/150‑1 (MQ=255)
                                                      ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                             <  2:5576703/150‑1 (MQ=255)
                                                      ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                             >  1:3698383/1‑150 (MQ=255)
                                                      ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                             <  1:4489009/150‑1 (MQ=255)
                                                       cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGc                                                                                            >  1:1391448/1‑150 (MQ=255)
                                                       cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGc                                                                                            >  1:1703898/1‑150 (MQ=255)
                                                       cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGc                                                                                            >  2:8502806/1‑150 (MQ=255)
                                                       cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGc                                                                                            >  2:7892061/1‑150 (MQ=255)
                                                         gAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGa                                                                                          >  1:2381017/1‑150 (MQ=255)
                                                          aaaTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAc                                                                                         <  2:5205529/150‑1 (MQ=255)
                                                           aaTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                        <  1:5251413/150‑1 (MQ=255)
                                                           aaTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                        <  2:4834914/150‑1 (MQ=255)
                                                            aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAccc                                                                                       >  2:5659920/1‑150 (MQ=255)
                                                             ttGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCt                                                                                      >  2:8534227/1‑150 (MQ=255)
                                                              tGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTa                                                                                     >  1:7010438/1‑150 (MQ=255)
                                                                 aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACtt                                                                                  <  1:6702877/150‑1 (MQ=255)
                                                                 aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACtt                                                                                  <  2:6346890/150‑1 (MQ=255)
                                                                   ttGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                                >  1:5784050/1‑150 (MQ=255)
                                                                           tGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTgcg                                                                        <  2:1801345/150‑1 (MQ=255)
                                                                             aaGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCa                                                                      <  1:3328363/150‑1 (MQ=255)
                                                                             aaGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCa                                                                      >  2:6653697/1‑150 (MQ=255)
                                                                              aGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAg                                                                     <  2:8103726/150‑1 (MQ=255)
                                                                               gTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgt                                                                    <  2:4420199/150‑1 (MQ=255)
                                                                                tGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgtg                                                                   <  1:2237203/150‑1 (MQ=255)
                                                                                   tGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGa                                                                >  1:4126947/1‑150 (MQ=255)
                                                                                    ggCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                               <  2:3426813/150‑1 (MQ=255)
                                                                                     gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACg                                                              >  2:7110483/1‑150 (MQ=255)
                                                                                     gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACg                                                              >  2:912816/1‑150 (MQ=255)
                                                                                       agagCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGa                                                            >  1:2567718/1‑150 (MQ=255)
                                                                                        gagCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCaa                                                                                                                                             <  2:7354587/68‑1 (MQ=38)
                                                                                        gagCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCaa                                                                                                                                             >  1:7354594/1‑68 (MQ=255)
                                                                                           cAGGCTTGCCTGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTg                                                        <  2:7735874/150‑1 (MQ=255)
                                                                                            aGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTgg                                                       >  2:3055056/1‑150 (MQ=255)
                                                                                                 tGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTg                                                  <  1:2457966/150‑1 (MQ=255)
                                                                                                  gCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTgg                                                 <  2:2674564/150‑1 (MQ=255)
                                                                                                  gCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTgg                                                 <  1:3569519/150‑1 (MQ=255)
                                                                                                   cATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                <  2:4106170/150‑1 (MQ=255)
                                                                                                   cATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                >  1:8416820/1‑150 (MQ=255)
                                                                                                   cATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                <  1:3622650/150‑1 (MQ=255)
                                                                                                        gTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACg                                           >  2:1671559/1‑150 (MQ=255)
                                                                                                        gTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCAGTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACg                                           <  1:6528455/150‑1 (MQ=255)
                                                                                                         tttGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACgg                                          >  2:3934226/1‑150 (MQ=255)
                                                                                                              gTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                     <  2:2057664/150‑1 (MQ=255)
                                                                                                              gTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                     >  1:7267679/1‑150 (MQ=255)
                                                                                                                tCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCtt                                   >  1:2192337/1‑150 (MQ=255)
                                                                                                                  gCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCa                                 <  1:7904336/150‑1 (MQ=255)
                                                                                                                  gCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCa                                 <  1:5661048/150‑1 (MQ=255)
                                                                                                                   cTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCAt                                <  2:6155699/150‑1 (MQ=255)
                                                                                                                    tGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATc                               >  1:6596676/1‑150 (MQ=255)
                                                                                                                     gTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCt                              <  2:8712443/150‑1 (MQ=255)
                                                                                                                       aTGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                            <  2:3267828/150‑1 (MQ=255)
                                                                                                                           gttgttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCg                        <  1:5578356/150‑1 (MQ=255)
                                                                                                                                 cGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc                  <  1:1483494/150‑1 (MQ=255)
                                                                                                                                 cGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc                  <  2:5443989/150‑1 (MQ=255)
                                                                                                                                  gACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                 <  1:6274218/150‑1 (MQ=255)
                                                                                                                                     cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              >  1:3767125/1‑150 (MQ=255)
                                                                                                                                     cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              >  2:7600918/1‑150 (MQ=255)
                                                                                                                                     cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              >  1:314013/1‑150 (MQ=255)
                                                                                                                                      aGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga             <  2:1399193/150‑1 (MQ=255)
                                                                                                                                      aGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga             >  2:8314509/1‑150 (MQ=255)
                                                                                                                                         gcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGAttt          >  2:4708858/2‑150 (MQ=255)
                                                                                                                                           ttGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACGGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAt        >  1:4196160/1‑150 (MQ=255)
                                                                                                                                           ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAt        >  2:114742/1‑150 (MQ=255)
                                                                                                                                           ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAt        >  1:3372262/1‑150 (MQ=255)
                                                                                                                                           ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTGGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAt        >  2:2169551/1‑150 (MQ=255)
                                                                                                                                              ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCgg     <  1:8565988/150‑1 (MQ=255)
                                                                                                                                               gCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgc    >  2:6865816/1‑150 (MQ=255)
                                                                                                                                                cATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGNTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgca   <  2:4576146/150‑1 (MQ=255)
                                                                                                                                                cATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgca   <  1:6853596/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3585475‑3585769

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: