Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,432 A→C 10.0% R726R (AGG→CGG)  PP_1897 → DNA internalization‑related competence protein ComEC/Rec2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,147,4320AC10.0% 486.0 / 2.0 198R726R (AGG→CGG) PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base A (96/69);  new base C (2/22);  total (98/100)
Fisher's exact test for biased strand distribution p-value = 3.26e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

TGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTATT  >  NC_002947/2147283‑2147580
                                                                                                                                                     |                                                                                                                                                    
tGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTa                                                                                                                                                      >  2:506044/1‑150 (MQ=255)
 gATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAg                                                                                                                                                     >  1:1807956/1‑150 (MQ=255)
  aTTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgg                                                                                                                                                    >  2:1265473/1‑150 (MQ=255)
  aTTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgg                                                                                                                                                    >  1:592671/1‑150 (MQ=255)
   tttCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAggg                                                                                                                                                   >  1:791848/1‑150 (MQ=255)
    ttCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCTGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGTGGGGGCTTTGGGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgggg                                                                                                                                                  <  2:2721702/150‑1 (MQ=255)
      cgcgGGGGCGCAACACCAGCTACGGGCACCCGCACGTGCAGGTGGTCGCGCGGTATGGGCGGCATGGGGGGGTGAGGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGc                                                                                                                                                <  2:2980966/150‑1 (MQ=255)
             gcgcAACAACTGCTTCGGGCACCCGCAGGTGCTGGTGGTCGAGCGTTTTGGGCGGCTTGGGGTGGTTATGAATGATACGGGGGAGGAGGGGGCATTGGGGTTGGTGCTGGGGCGGCTTGGGGTAGTCGAGGGGGTTAGGGGGCAgcggcg                                                                                                                                         <  2:82924/150‑1 (MQ=255)
             gcgcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAgcggcg                                                                                                                                         >  1:4603543/1‑150 (MQ=255)
             gcgcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGTGCAGcggcgg                                                                                                                                        >  2:2075534/1‑150 (MQ=255)
                 aacaacAGCTTCGGGCACCCGCACGTGCTGGTGGTCTAGCGTTTTGGGCGGCTTGGGGTTGTTATGATTGATACGGCGGTGGTGGGGGCATTGCGGTTGGTGCTGGGGCGGCTTGGGTTATTCGTGGGGTTTAGGGGGCAGCGGCGGTTc                                                                                                                                     <  2:2882207/150‑1 (MQ=255)
                  aaaactGCTTGGGGCACCGGCAGGTGCTGGGGGTGTCGGGTTTTGGGGGGATGGGGGTGTTTTTGATTGATAGGGGGGTGGTGGGGGATTTGGGGTTGGTGGTGGGGCGGCTTGGGGTTGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCt                                                                                                                                    <  1:4833060/148‑1 (MQ=255)
                   caacaGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTg                                                                                                                                   >  2:5506740/1‑150 (MQ=255)
                     actgCTTGGGGCACCCGCACGTGCTGGTGGTCGAGCGGTTTGGGCGGCTTGGGGTGGTTATGCTTGATACGGCGGTGGAGGGGGCATTGGGGTTGGTGCTGGGGCGGCTTGGGTTAGTCGTGGGGGTTAGGGGGCAGCGGCGGTTCTGGc                                                                                                                                 <  2:2529082/147‑1 (MQ=255)
                      caGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCg                                                                                                                                >  2:4374387/1‑150 (MQ=255)
                       cgCTTCGGGCACCCGCACGTGCCGGTGGTCGAGCGGTATGGGCGGCATGGGGGGGTGAGGCATGAGAGGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCgg                                                                                                                               <  2:3323508/149‑1 (MQ=255)
                          ctcgaccACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGGTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGtt                                                                                                                            >  2:5106820/7‑150 (MQ=255)
                           tCGGGCCCCCGCAGGTGCTGGTGGTCGCGGGTTTTGGGGGGCTTGGGGTGGTTATGCTTGATAGGGCGGTGGTGGGGGCTTTGGGGTTGGTGGTGGGGGGGCATGGGGTAGCCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTc                                                                                                                           <  1:1570645/150‑1 (MQ=255)
                             gggCCCCCGCCGGTGCTGGGGGTCGCGGGGTTGGGGGGGCTGGGGGTGGTGTTGATTTATCGGGGGGTGGTGGGGGATTGGCGGTGGGTGCTGGGGCGGCTTGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                                                                                                                         <  1:4162819/150‑1 (MQ=255)
                             gggCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                                                                                                                         >  2:4992666/1‑150 (MQ=255)
                             gaccACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                                                                                                                         >  2:1026576/4‑150 (MQ=255)
                               gCACCCGCACGTGCTGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGAGGCATGATAGGGCGGTGGCGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTCGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAgg                                                                                                                       <  2:584401/150‑1 (MQ=255)
                               gCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCAGTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAgg                                                                                                                       >  2:3250175/1‑150 (MQ=255)
                               ccACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAgg                                                                                                                       >  1:3448531/2‑150 (MQ=255)
                                cccccGCACGTGCCGGGGGTCGCGCGGGTGGGGCGGCATGGGGGGGTGATGCCTGAGCCGGCGGAGGAGGGGGCATTGGGGTTGGTGCGGGGGCGGCAGGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg                                                                                                                      <  2:5553623/148‑1 (MQ=255)
                                cACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg                                                                                                                      >  1:4646648/1‑150 (MQ=255)
                                cACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg                                                                                                                      >  1:2928625/1‑150 (MQ=255)
                                cACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg                                                                                                                      >  1:865894/1‑150 (MQ=255)
                                  cccGCACGTGCCGGTGGGCGCGGGGTTTGGGCGGCATGGGGTGGTTTTGCATGATACGGGGGAGGAGGGGGCATTGGGGTTGGTGCTGGGGCGGCTTGGGGTTGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa                                                                                                                    <  2:1924095/150‑1 (MQ=255)
                                  cccGCACGTGCAGGTGGTCGACCGGTATGGCCGGAATGGGTTGGTTATGAATGATACGGCGGAGTAGGGGGCATTGCGTTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCAGGTTCGAGGGGa                                                                                                                    <  2:4293520/150‑1 (MQ=255)
                                   ccGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAg                                                                                                                   >  1:984799/1‑150 (MQ=255)
                                   ccGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGGGAGGGGGCAGCGGCGGTTCTGGCGGGTGCGAGGGGAg                                                                                                                   >  2:4718830/1‑150 (MQ=255)
                                     gCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTg                                                                                                                 >  1:3395000/1‑150 (MQ=255)
                                      caggTGCTGGTGGTCGCGCGGTATGGGCGGCATGGGGTGGTGATGCATGATCGGGGGGTGGTGGGGGCATTGGGGTTGGTGGTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTga                                                                                                                <  1:4199556/147‑1 (MQ=255)
                                      cACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTga                                                                                                                >  2:2455946/1‑150 (MQ=255)
                                       aggTGCTGGTGGTGGTGGGTTTTGGGGGGCTTGGGGTGGTGATGCTTTTTCGGGGGGTGGCGGGGGCTTTGGGGTTGGTGGTGGGGGGGCTTGGGGTAGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgag                                                                                                               <  1:3244162/148‑1 (MQ=255)
                                       aCGTGCTGGTGGTCGAGCGGTCTGGGCGGCTGGGGGGGGTGATGCTGGAGCGGGCGGAGGAGGGGGCTTTGCGGTTGGTGCTGGGGCGGCAGGGGGTCGTCGCGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgag                                                                                                               <  2:70242/150‑1 (MQ=255)
                                               gtggtCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgaggcagaggc                                                                                                       >  1:3724291/1‑150 (MQ=255)
                                               gtggtCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgaggcagaggc                                                                                                       >  1:384861/1‑150 (MQ=255)
                                                 ggtCGCGCGGTAGGGGCGGCTGGGGGTGGTTTTGCTTGATCGGGGGGCGGTGGGGGCATGGCGGTTGGTGCGGGGGCGGCTGGGGGTCGGCGCGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGaggcagaggcag                                                                                                     <  1:4461358/150‑1 (MQ=255)
                                                  gtCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGaggcagg                                                                                                    >  2:3900797/1‑150 (MQ=255)
                                                  gtCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGaggcagg                                                                                                    >  1:1454028/1‑150 (MQ=255)
                                                  gtCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGaggcagg                                                                                                    >  2:1553689/1‑150 (MQ=255)
                                                     caGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCAt                                                                                                 >  2:387708/2‑150 (MQ=255)
                                                       gggtttgggggggCTTGGGGTGTTGTTGATTGATCGGGGGGTGGTGGGGGATTTGGGGTTGGTGGTGGGGGGGATTGGGGTTGTGGTGGGGTTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCt                                                                                               <  1:2574010/143‑1 (MQ=255)
                                                       gCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCt                                                                                               >  1:5317882/1‑150 (MQ=255)
                                                         ggttggggCGGCCTGGGGTGGTGCTGCTTGATCCGGGGGCGGCGGGGGCTTTGGGGTTGGGGGTGGGGCGGCTTGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGc                                                                                             <  2:4856830/145‑1 (MQ=255)
                                                               ggcggcAGGGGGGGGTGATGCATGATCCGGCGGTGGCGGGGGCATTGCGGTTGGTGCGGGGGCGGCAGGGGGTAGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGccc                                                                                       <  2:5401428/150‑1 (MQ=255)
                                                                gcggcCTGGGGTGGTGATGCTTGATACGGCGGAGGCGGGGGCATTGCGGTTGGTGCTGGGGCGGCTTGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCa                                                                                      <  2:4127648/150‑1 (MQ=255)
                                                                     aTGGGGTGTTGATGCTTTATCGGGGGGTGTTGGGGGCATGGGGTTTGGGGGGGGGGCGGATGGGGGTCGCGGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGAccc                                                                                 <  1:5066229/150‑1 (MQ=255)
                                                                        gggTGGTGTTGAATGTTACGGGGGTGGAGGGGGATTTGGGGTTGGTGCTGGGGCGGCTTGGGGTTGTCTCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGAGCCCaa                                                                              <  2:3618194/150‑1 (MQ=255)
                                                                        gggTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCaa                                                                              >  2:5248497/1‑150 (MQ=255)
                                                                          gtggtgATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAAtt                                                                            >  1:5010937/1‑150 (MQ=255)
                                                                             gtgATGCATGATACGGCGGAGGCGGGGGCTTTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTccc                                                                         <  2:4170147/150‑1 (MQ=255)
                                                                              tgATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCt                                                                        >  2:2735747/1‑150 (MQ=255)
                                                                              ggcagCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGTGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCt                                                                        <  2:1718741/146‑1 (MQ=255)
                                                                                  gCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGAcc                                                                    >  1:1877538/1‑150 (MQ=255)
                                                                                  gCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGAcc                                                                    >  2:532783/1‑150 (MQ=255)
                                                                                    aTGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCtt                                                                  >  2:4960862/1‑150 (MQ=255)
                                                                                     tGATACGGCGGAGGAGGGGGCTTTGCGGTTGGTGCGGGGGCGGCATGGGGTAGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTc                                                                 <  2:399538/150‑1 (MQ=255)
                                                                                       aTACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCgg                                                               <  2:5010937/150‑1 (MQ=255)
                                                                                            gcggAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGACGAGGGGGTTAGGGGGCAGCGGCGGTTATGGCGGGTTCGAGGGGAGTGAGGAAGAGGCAGGCATCTGCCAGCCCAGACCACAATTAACTGACATACGGCAGAt                                                          >  1:1375280/1‑150 (MQ=255)
                                                                                               gtggggggggCTTGGGGGTTGGTGGGGGGGCGGCTTGGGGTCGGCGCGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAg                                                       <  1:3822673/145‑1 (MQ=255)
                                                                                               gtggaggGGGCTTTGCGGTGGGTGCTGGGGGGGATGGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAg                                                       <  2:5167026/148‑1 (MQ=255)
                                                                                               gaggcgGGGGCTTTGGGGTTGGTGCTGGGGCGGCATGGGGTTGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAg                                                       <  2:984799/150‑1 (MQ=255)
                                                                                                 ggaggGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGcc                                                     >  1:3352823/1‑150 (MQ=255)
                                                                                                  gaggGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGccc                                                    >  1:2533786/1‑150 (MQ=255)
                                                                                                    gggggCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTc                                                  >  2:2403532/1‑150 (MQ=255)
                                                                                                           tggggtgtggtggGGGGGCGGCTGGGGGTTGCGGCGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGgcgcgc                                           <  1:4473834/143‑1 (MQ=255)
                                                                                                                gtgggTGCGGGGGCGGCTTGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATg                                      <  2:4876672/147‑1 (MQ=255)
                                                                                                                gTTGGTGCGGGGGGGGCTGGGGGTCGGGGGGGGGTTTGGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATg                                      <  1:5269998/150‑1 (MQ=255)
                                                                                                                     tGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTaga                                 <  2:1207456/150‑1 (MQ=255)
                                                                                                                       ggggggggGCTGGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGt                               <  1:4374419/148‑1 (MQ=255)
                                                                                                                          gggCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTggc                            >  1:3317329/1‑150 (MQ=255)
                                                                                                                          gggCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTggc                            >  2:1946550/1‑150 (MQ=255)
                                                                                                                            gcggcATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTggcgg                          >  1:1041068/1‑150 (MQ=255)
                                                                                                                             gggcGTGGGGTCGGCGTGGGGGTTTGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTggcggc                         <  1:1274488/149‑1 (MQ=255)
                                                                                                                             gggcATGGGTCCGTCGCGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTggcggc                         <  1:4794056/149‑1 (MQ=255)
                                                                                                                               gcATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCt                       >  2:2008867/1‑150 (MQ=255)
                                                                                                                               gcATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCt                       >  1:3612366/1‑150 (MQ=255)
                                                                                                                                cATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCtt                      >  1:1906418/1‑150 (MQ=255)
                                                                                                                                 aTGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCttt                     <  2:3395000/150‑1 (MQ=255)
                                                                                                                                 aTGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCttt                     <  2:612891/150‑1 (MQ=255)
                                                                                                                                    gggTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTc                  <  2:3451391/150‑1 (MQ=255)
                                                                                                                                    gggTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTc                  <  2:4415995/150‑1 (MQ=255)
                                                                                                                                     ggTCGTCGCGGGGGTTTGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTcc                 <  1:1770187/150‑1 (MQ=255)
                                                                                                                                     ggTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTcc                 >  2:1474271/1‑150 (MQ=255)
                                                                                                                                       tAGTCGAGGGGGTTAGGTGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGaa              <  1:2075534/150‑1 (MQ=255)
                                                                                                                                        tgtgggggggtTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGaa              <  1:5170999/144‑1 (MQ=255)
                                                                                                                                         ggggggggggTTGGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAg             <  1:1242220/145‑1 (MQ=255)
                                                                                                                                         gTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAg             >  2:5342163/1‑150 (MQ=255)
                                                                                                                                          tCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAgg            >  2:5569820/1‑150 (MQ=255)
                                                                                                                                            gAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGAACTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAgggg          >  1:4617090/1‑150 (MQ=255)
                                                                                                                                               ggggTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGt       >  1:2710445/1‑150 (MQ=255)
                                                                                                                                               ggggTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGt       >  1:5532709/1‑150 (MQ=255)
                                                                                                                                                 ggTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGttt     <  2:4603543/150‑1 (MQ=255)
                                                                                                                                                  gTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTa    >  1:2012176/1‑150 (MQ=255)
                                                                                                                                                  gTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTa    <  2:3448531/150‑1 (MQ=255)
                                                                                                                                                    tAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTAtt  <  2:2180630/150‑1 (MQ=255)
                                                                                                                                                    tAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTAtt  <  2:3799993/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAGGCAGGCATCTGCCAGCCCAGACCCCAATTCCCTGACCTTCGGCAGATGAGCCCTCGGCGCGCCTATGGTAGAGTGGCGGCCTTTTTCCGAAGGGGCGTTTATT  >  NC_002947/2147283‑2147580

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