Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,396:1 +T 7.9% coding (1542/2016 nt) uvrB → uvrABC excinuclease protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,238,3961.T 7.9% 662.3 / 52.7 164coding (1542/2016 nt)uvrBuvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base . (76/75);  new base T (8/5);  total (84/80)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 166 total aligned reads displayed.

GCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAGCGAGCCATCGA  >  NC_002947/2238247‑2238545
                                                                                                                                                      |                                                                                                                                                     
  cGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCt                                                                                                                                                     <  2:4686106/150‑1 (MQ=255)
  cGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCt                                                                                                                                                     <  2:4686170/150‑1 (MQ=255)
   gACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtg                                                                                                                                                    >  1:1678431/1‑150 (MQ=255)
       acgacgTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtggtgg                                                                                                                                                >  2:4441148/1‑150 (MQ=255)
         gacgTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGcc                                                                                                                                              <  2:3160316/150‑1 (MQ=255)
          acgTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACGTGCCCGAGGTGTCGCTGGTGGCCa                                                                                                                                             <  1:3996156/150‑1 (MQ=255)
            gTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATc                                                                                                                                           >  1:1380601/1‑150 (MQ=255)
                     cGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCg                                                                                                                                  <  1:4959182/150‑1 (MQ=255)
                       cTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGa                                                                                                                                >  2:1130851/1‑150 (MQ=255)
                         aCCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACa                                                                                                                              >  1:2155104/1‑150 (MQ=255)
                          ccTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACac                                                                                                                             >  2:3527577/1‑149 (MQ=255)
                           cTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaag                                                                                                                            >  1:1136212/1‑150 (MQ=255)
                              cACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaa                                                                                                                         <  1:2050265/150‑1 (MQ=255)
                              cACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaa                                                                                                                         <  2:3077540/150‑1 (MQ=255)
                                cTCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaagg                                                                                                                       >  1:885340/1‑150 (MQ=255)
                                 tCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaaggc                                                                                                                      >  1:2400998/1‑149 (MQ=255)
                                 tCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaaggc                                                                                                                      >  1:276246/1‑149 (MQ=255)
                                 tCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGAACTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaaggc                                                                                                                      >  1:2172267/1‑149 (MQ=255)
                                     aCATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGACGGTTTTc                                                                                                                  >  2:4023470/1‑150 (MQ=255)
                                     aCATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGcttcc                                                                                                                  <  2:1832817/150‑3 (MQ=255)
                                      cATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCt                                                                                                                 <  2:1835908/150‑1 (MQ=255)
                                      cATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCt                                                                                                                 <  2:1836036/150‑1 (MQ=255)
                                       aTCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTg                                                                                                                >  2:1740552/1‑150 (MQ=255)
                                         cGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCg                                                                                                              <  2:5280452/150‑1 (MQ=255)
                                          gACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGt                                                                                                             <  1:4039396/150‑1 (MQ=255)
                                           acacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGtt                                                                                                            >  1:77871/1‑150 (MQ=255)
                                           acacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGtt                                                                                                            <  1:2905871/150‑1 (MQ=255)
                                           acacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGtt                                                                                                            >  1:55163/1‑150 (MQ=255)
                                             acGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTcc                                                                                                          <  1:5344909/150‑1 (MQ=255)
                                               ggTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtggtgg                                                                                                                                                <  1:4230232/110‑1 (MQ=255)
                                                    aGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCt                                                                                                   <  1:5014710/150‑1 (MQ=255)
                                                      cGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTcc                                                                                                 <  2:2053346/150‑1 (MQ=255)
                                                      cGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTcc                                                                                                 <  1:1740552/150‑1 (MQ=255)
                                                            gagaTGATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATc                                                                                           <  2:1370567/150‑1 (MQ=255)
                                                              gaTTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCa                                                                                         >  2:575744/1‑150 (MQ=255)
                                                              gaTTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCa                                                                                         >  2:555570/1‑150 (MQ=255)
                                                                ttATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGa                                                                                       <  2:1413600/150‑1 (MQ=255)
                                                                ttATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGa                                                                                       >  2:2502757/1‑150 (MQ=255)
                                                                    ccGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATCTTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCag                                                                                   >  2:4525924/1‑149 (MQ=255)
                                                                    ccGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCAt                                                                                   >  1:2335702/1‑150 (MQ=255)
                                                                    ccGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCAt                                                                                   <  1:3212174/150‑1 (MQ=255)
                                                                      gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                                 >  2:3048459/1‑150 (MQ=255)
                                                                      gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                                 >  1:2289647/1‑150 (MQ=255)
                                                                         atcTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCATCggtc                                                                              <  2:3572417/148‑3 (MQ=255)
                                                                           cTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCcgc                                                                            >  1:4545213/1‑150 (MQ=255)
                                                                            tGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCcgcg                                                                           >  2:3184700/1‑150 (MQ=255)
                                                                            tGCGCCTGGGCACCTACGACGTGCTGGTGGGCATCAACCTGCAGCGCGAGGGCCTGGACATGACCGAGGTGTCGCTGGTTGACATCCTCGATGAGGACAAGGAAGGTTTTCAGCGTTCCGAGCGCTCCCTCATCCAGACCAAAGGCcgcg                                                                           >  2:5002928/1‑150 (MQ=255)
                                                                               gcCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCgg                                                                        >  1:3092937/1‑150 (MQ=255)
                                                                                ccTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGc                                                                       <  1:3859265/150‑1 (MQ=255)
                                                                                ccTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGc                                                                       >  2:3296965/1‑150 (MQ=255)
                                                                                ccTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGc                                                                       >  1:4659281/1‑150 (MQ=255)
                                                                                  tGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGTGCTCCCTGATCCAGACCATTGGCCGCGCGGCAc                                                                     >  1:5449068/1‑150 (MQ=255)
                                                                                  tGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCAc                                                                     <  2:2271155/150‑1 (MQ=255)
                                                                                   gggCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACg                                                                    >  1:1999016/1‑150 (MQ=255)
                                                                                    ggCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGt                                                                   <  2:2289647/150‑1 (MQ=255)
                                                                                    ggCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGt                                                                   >  1:1080905/1‑150 (MQ=255)
                                                                                         cTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCt                                                              >  1:2564559/1‑150 (MQ=255)
                                                                                         cTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGAAACCt                                                              >  1:3036921/1‑150 (MQ=255)
                                                                                            cGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCaa                                                           <  2:328332/150‑1 (MQ=255)
                                                                                            cGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCaa                                                           <  2:2564559/150‑1 (MQ=255)
                                                                                              aCGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATg                                                         <  2:3341429/150‑1 (MQ=255)
                                                                                              aCGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATg                                                         <  1:1087583/150‑1 (MQ=255)
                                                                                                    tggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCggg                                                   >  2:3697420/1‑150 (MQ=255)
                                                                                                    tggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCggg                                                   >  2:4292418/1‑150 (MQ=255)
                                                                                                     ggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGc                                                  <  2:3782419/150‑1 (MQ=255)
                                                                                                     ggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGc                                                  <  2:3782844/150‑1 (MQ=255)
                                                                                                     ggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGc                                                  <  1:2407460/150‑1 (MQ=255)
                                                                                                      gtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCAtt                                                                                  >  2:1456639/1‑117 (MQ=255)
                                                                                                      gtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCAtt                                                                                  <  1:1456639/117‑1 (MQ=255)
                                                                                                      gtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGcc                                                 >  1:51334/1‑150 (MQ=255)
                                                                                                       tggGCATCAACCTGCGGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCa                                                >  1:2404443/1‑150 (MQ=255)
                                                                                                          gCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTc                                             >  2:3517007/1‑150 (MQ=255)
                                                                                                          gCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTc                                             >  2:2012575/1‑150 (MQ=255)
                                                                                                            aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                                           >  2:4604611/1‑150 (MQ=255)
                                                                                                            aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                                           >  1:2220342/1‑150 (MQ=255)
                                                                                                            aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                                           <  1:1226507/150‑1 (MQ=255)
                                                                                                                 ccTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGc                                      <  2:1080905/150‑1 (MQ=255)
                                                                                                                   tgctgcGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCg                                    >  2:4837967/1‑150 (MQ=255)
                                                                                                                      tgcGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATa                                 <  1:335794/150‑1 (MQ=255)
                                                                                                                      tgcGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATa                                 <  1:335446/150‑1 (MQ=255)
                                                                                                                      tgcGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATa                                 <  1:2672127/150‑1 (MQ=255)
                                                                                                                       gcgcGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATac                                >  2:5391514/1‑149 (MQ=255)
                                                                                                                       gcgcGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATaa                                >  2:799961/1‑150 (MQ=255)
                                                                                                                        cgcgAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAAc                               <  2:2751686/150‑1 (MQ=255)
                                                                                                                         gcgAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATaa                                >  2:2455019/1‑148 (MQ=255)
                                                                                                                         gcgAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATaa                                <  1:2455019/148‑1 (MQ=255)
                                                                                                                            aGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTaa                           <  2:2220342/150‑1 (MQ=255)
                                                                                                                             gggCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAAc                          <  2:905564/150‑1 (MQ=255)
                                                                                                                                ccTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAgg                       >  2:4203129/1‑150 (MQ=255)
                                                                                                                                 cTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCt                                            >  2:1789084/1‑128 (MQ=255)
                                                                                                                                    gACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTcc                   >  2:3221215/1‑150 (MQ=255)
                                                                                                                                     aCATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCa                  >  1:788203/1‑150 (MQ=255)
                                                                                                                                      cATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCAt                 >  2:2456736/1‑150 (MQ=255)
                                                                                                                                       aTGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATg                <  1:4441148/150‑1 (MQ=255)
                                                                                                                                        tGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGc               <  2:3138698/150‑1 (MQ=255)
                                                                                                                                          cccGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAg             >  2:4698379/1‑150 (MQ=255)
                                                                                                                                          cccGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAg             >  2:255861/1‑150 (MQ=255)
                                                                                                                                             gAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAGCGa          <  2:822570/150‑1 (MQ=255)
                                                                                                                                              aGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAGCGAg         <  2:1966215/150‑1 (MQ=255)
                                                                                                                                                     gCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAGCGAGCCATcga  <  2:2167162/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATGCAGCGAGCCATCGA  >  NC_002947/2238247‑2238545

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: