Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 35,009 A→G 6.7% intergenic (+27/+31) czrSA → / ← PP_0031 metal‑sensor histidine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294735,0090AG 6.7% 615.0 / 21.6 194intergenic (+27/+31)czrSA/PP_0031metal‑sensor histidine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (90/91);  new base G (6/7);  total (96/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 200 total aligned reads displayed.

AGTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTCT  >  NC_002947/34861‑35158
                                                                                                                                                    |                                                                                                                                                     
aGTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCaa                                                                                                                                                      >  2:1200956/1‑150 (MQ=255)
 gTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAAt                                                                                                                                                     >  1:1415260/1‑150 (MQ=255)
     aTCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGAc                                                                                                                                                 >  1:2840379/1‑150 (MQ=255)
     aTCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGAc                                                                                                                                                 >  2:217277/1‑150 (MQ=255)
      tCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACg                                                                                                                                                >  1:5038551/1‑150 (MQ=255)
       cATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACgg                                                                                                                                               >  2:2307316/1‑150 (MQ=255)
       cATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                               <  1:4790668/150‑1 (MQ=255)
       cATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                               <  2:1076628/150‑1 (MQ=255)
         tGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                             <  2:1008219/150‑1 (MQ=255)
          gAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGcc                                                                                                                                            <  1:1355126/150‑1 (MQ=255)
            gTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGccgt                                                                                                                                          >  2:5480444/1‑148 (MQ=255)
            gTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCtt                                                                                                                                          >  1:2421810/1‑150 (MQ=255)
              cTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGCCtaat                                                                                                                                        <  1:1533434/150‑4 (MQ=255)
              cTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCtttt                                                                                                                                        <  1:3462715/150‑1 (MQ=255)
                gCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgt                                                                                                                                      <  2:107231/150‑1 (MQ=255)
                 cACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgtt                                                                                                                                     >  2:4604740/1‑150 (MQ=255)
                         aTCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCaa                                                                                                                             >  1:5499825/1‑150 (MQ=255)
                            ggCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGc                                                                                                                          <  1:4602645/150‑1 (MQ=255)
                            ggCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATACAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGc                                                                                                                          >  1:2827569/1‑150 (MQ=255)
                             gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCa                                                                                                                         >  1:397640/1‑150 (MQ=255)
                             gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCa                                                                                                                         >  1:2858232/1‑150 (MQ=255)
                              cGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCAt                                                                                                                        <  2:3988246/150‑1 (MQ=255)
                                aGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCa                                                                                                                      <  2:3121505/150‑1 (MQ=255)
                                  cGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGt                                                                                                                    <  2:5223311/150‑1 (MQ=255)
                                        ggTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa                                                                                                              <  1:1646497/150‑1 (MQ=255)
                                        ggTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa                                                                                                              <  1:5013670/150‑1 (MQ=255)
                                         gTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAg                                                                                                             <  1:230321/150‑1 (MQ=255)
                                             ccGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGAcc                                                                                                         >  2:1625677/1‑150 (MQ=255)
                                              cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttg                                                                                                                                       <  1:3426771/119‑1 (MQ=255)
                                              cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACca                                                                                                        >  1:3896113/1‑150 (MQ=255)
                                                 agaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAg                                                                                                     <  2:3501080/150‑1 (MQ=255)
                                                  gaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAgg                                                                                                    >  1:1593069/1‑150 (MQ=255)
                                                   aaaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAggg                                                                                                   <  2:1822565/150‑1 (MQ=255)
                                                   aaaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAggg                                                                                                   >  2:4743562/1‑150 (MQ=255)
                                                   aaaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAggg                                                                                                   >  2:5548439/1‑150 (MQ=255)
                                                    aaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGt                                                                                                  <  1:4753326/150‑1 (MQ=255)
                                                       cTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCATTGCCGAAGGACCACAGGGTCCa                                                                                               >  2:3436544/1‑150 (MQ=255)
                                                        ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGTACCACAGGGTCCAt                                                                                              >  2:3414855/1‑150 (MQ=255)
                                                         tCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATg                                                                                             >  1:896252/1‑150 (MQ=255)
                                                           aCCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATggg                                                                                           <  2:265661/150‑1 (MQ=255)
                                                           aCCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATggg                                                                                           <  2:5128048/150‑1 (MQ=255)
                                                            ccTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGc                                                                                          <  1:4066596/150‑1 (MQ=255)
                                                            ccTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGc                                                                                          <  2:3370427/150‑1 (MQ=255)
                                                                   ccACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCa                                                                                   >  1:2119251/1‑150 (MQ=255)
                                                                     aCCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCt                                                                                                                                <  2:2994148/103‑1 (MQ=255)
                                                                     aCCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGa                                                                                 <  1:5316144/150‑1 (MQ=255)
                                                                        ttCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTATGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACAACAGGGTCCATGGGCGACCGCAGAATc                                                                              >  2:3291902/1‑150 (MQ=255)
                                                                        ttCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATc                                                                              >  1:4489298/1‑150 (MQ=255)
                                                                           aCTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAcc                                                                           >  2:634129/1‑150 (MQ=255)
                                                                             tttCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGa                                                                         >  2:1447359/1‑150 (MQ=255)
                                                                             tttCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGa                                                                         >  2:769194/1‑150 (MQ=255)
                                                                             tttCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGa                                                                         <  1:426433/150‑1 (MQ=255)
                                                                                  ggAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTtaa                                                                    <  1:1024979/150‑1 (MQ=255)
                                                                                  ggAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTtaa                                                                    >  2:674941/1‑150 (MQ=255)
                                                                                         cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                             <  2:4420969/150‑1 (MQ=255)
                                                                                         cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                             >  1:4197469/1‑150 (MQ=255)
                                                                                          aGCAAATGCAGCCGGCAATGACAGCCGCTTAAATTGCGCAGACACGATTAGTCCCGTCAATCGACGGGCCTTTTGTTTGACTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCATAATCACCGAATTAATAAtctta                                                            >  1:520325/1‑149 (MQ=255)
                                                                                           ccAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATGAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtcttct                                                           <  2:4129837/149‑1 (MQ=255)
                                                                                               aTGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTc                                                       <  1:3006427/150‑1 (MQ=255)
                                                                                                 gCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGt                                                     <  2:328836/150‑1 (MQ=255)
                                                                                                  cAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTc                                                    <  2:144745/150‑1 (MQ=255)
                                                                                                   aGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCa                                                   >  2:5274655/1‑150 (MQ=255)
                                                                                                   aGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCa                                                   <  1:808396/150‑1 (MQ=255)
                                                                                                   aGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCa                                                   <  2:2737205/150‑1 (MQ=255)
                                                                                                   aGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGGTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCa                                                   >  1:7695/1‑150 (MQ=255)
                                                                                                      cGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGa                                                <  1:634129/150‑1 (MQ=255)
                                                                                                          aaTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTgg                                            >  1:2761730/1‑150 (MQ=255)
                                                                                                           aTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                           >  2:3914373/1‑150 (MQ=255)
                                                                                                           aTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                           <  2:967912/150‑1 (MQ=255)
                                                                                                           aTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                           >  1:5384025/1‑150 (MQ=255)
                                                                                                              aCCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGtt                                        <  2:825222/150‑1 (MQ=255)
                                                                                                               ccgcctCTTATTTTGCGCAGACACGTTTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGTCCACTGGGTCCATGGGCGACCTCAGAATCACCGACTTAATATTCTTCTGTTCGTCTGGATTGGCGttt                                       <  2:3893760/150‑1 (MQ=255)
                                                                                                                cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGtttt                                      >  2:2624294/1‑150 (MQ=255)
                                                                                                                cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGtttt                                      >  1:3029096/1‑150 (MQ=255)
                                                                                                                cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGtttt                                      <  1:4431251/150‑1 (MQ=255)
                                                                                                                  ccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTc                                    <  1:1200956/150‑1 (MQ=255)
                                                                                                                  ccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTc                                    <  1:872324/150‑1 (MQ=255)
                                                                                                                  ccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTc                                    <  1:5396343/150‑1 (MQ=255)
                                                                                                                         aaTTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGgt                             >  2:2888459/1‑150 (MQ=255)
                                                                                                                           ttGCGAAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCCACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGgtgt                           <  2:3806309/150‑1 (MQ=255)
                                                                                                                            tGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGtgtg                          >  2:5236509/1‑150 (MQ=255)
                                                                                                                               gcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTgaga                       <  2:117522/150‑1 (MQ=255)
                                                                                                                               gcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTACATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTgaga                       <  1:2049923/150‑1 (MQ=255)
                                                                                                                                 aGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAg                     <  2:5442171/150‑1 (MQ=255)
                                                                                                                                 aGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAg                     <  1:3117930/150‑1 (MQ=255)
                                                                                                                                 aGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAg                     <  2:5367226/150‑1 (MQ=255)
                                                                                                                                   acacGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTa                   >  1:1657843/1‑150 (MQ=255)
                                                                                                                                    cacGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTAc                  <  2:1175293/150‑1 (MQ=255)
                                                                                                                                       gATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGa               <  1:4001544/150‑1 (MQ=255)
                                                                                                                                        aaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAg              >  1:1789955/4‑150 (MQ=255)
                                                                                                                                        aaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAg              >  1:2578355/4‑150 (MQ=255)
                                                                                                                                         aaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGc             >  1:5614247/3‑150 (MQ=255)
                                                                                                                                           aGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCAt           <  1:4397948/150‑1 (MQ=255)
                                                                                                                                            ggCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCAtt          >  1:3565923/1‑150 (MQ=255)
                                                                                                                                             gCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTc         >  1:2138263/1‑150 (MQ=255)
                                                                                                                                               ccGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCAc       >  1:3632321/1‑150 (MQ=255)
                                                                                                                                                 gTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACtt     >  1:4737497/1‑150 (MQ=255)
                                                                                                                                                    aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTtct  <  1:217277/150‑1 (MQ=255)
                                                                                                                                                    aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTtct  >  2:1980881/1‑150 (MQ=255)
                                                                                                                                                    aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTtct  <  2:2421810/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTCT  >  NC_002947/34861‑35158

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: