Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,189,035 A→G 14.7% intergenic (+25/+17) PP_1038 → / ← tatC‑I hypothetical protein/Sec‑independent protein translocase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,189,0350AG14.7% 160.3 / 18.5 68intergenic (+25/+17)PP_1038/tatC‑Ihypothetical protein/Sec‑independent protein translocase protein
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base G (2/8);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGGT  >  NC_002947/1188894‑1189178
                                                                                                                                              |                                                                                                                                               
tGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCcgc                                                                                                                                          >  4:158250/1‑149 (MQ=255)
tGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCcgc                                                                                                                                          <  5:190664/149‑1 (MQ=255)
   aaGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaa                                                                                                                                       <  2:80305/149‑1 (MQ=255)
    aGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaaa                                                                                                                                      <  6:5232/149‑1 (MQ=255)
     ggTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAg                                                                                                                                     >  3:113195/1‑149 (MQ=255)
     ggTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGACTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaat                                                                                                                                      <  5:281387/148‑2 (MQ=255)
      gTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatag                                                                                                                                    <  4:315091/149‑4 (MQ=255)
      gTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatag                                                                                                                                    <  2:295417/149‑4 (MQ=255)
       tGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatagc                                                                                                                                   <  5:221930/149‑5 (MQ=255)
        ggAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatagcc                                                                                                                                  >  2:221157/1‑144 (MQ=255)
         gAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatagcc                                                                                                                                  <  1:221157/148‑6 (MQ=255)
          aGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaatagccgt                                                                                                                                <  6:79802/149‑8 (MQ=255)
                     gCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAg                                                                                                                      <  8:330830/148‑1 (MQ=255)
                     gCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGc                                                                                                                     <  4:280986/149‑1 (MQ=255)
                         aCGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCcac                                                                                                                  <  3:158250/148‑1 (MQ=255)
                          cGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACaa                                                                                                                <  2:264602/149‑1 (MQ=255)
                           gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAt                                                                                                               >  4:245728/1‑149 (MQ=255)
                           gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAt                                                                                                               >  4:217203/1‑149 (MQ=255)
                           gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAt                                                                                                               >  7:363481/1‑149 (MQ=255)
                                  ttGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTg                                                                                                        <  5:133306/149‑1 (MQ=255)
                                  ttGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTg                                                                                                        <  8:363481/149‑1 (MQ=255)
                                      gCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTg                                                                                                    >  4:251673/1‑149 (MQ=255)
                                       cAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                   >  6:31305/1‑149 (MQ=255)
                                       cAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                   <  7:160572/149‑1 (MQ=255)
                                       cAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                   <  1:290963/149‑1 (MQ=255)
                                          gtcgtcGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGc                                                                                                                     <  7:249452/128‑1 (MQ=255)
                                          gtcgtcGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGc                                                                                                                     >  8:249452/1‑128 (MQ=255)
                                              tcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                   >  3:135756/1‑142 (MQ=255)
                                              tcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                   <  4:135756/142‑1 (MQ=255)
                                               cgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGcc                                                                                           >  4:297443/1‑149 (MQ=255)
                                                gCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCg                                                                                          >  8:249122/1‑149 (MQ=255)
                                                    cgcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCtt                                                                                       <  3:1778/148‑1 (MQ=255)
                                                      cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCaa                                                                                    <  4:113195/149‑1 (MQ=255)
                                                      cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGACTCGCGGCTCGCCCAGTTCTAGAACGCAACACGGCTATTCGCGGGCACGCCCGAGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCaa                                                                                    >  2:102977/1‑149 (MQ=255)
                                                        cAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCt                                                                                                     >  1:83865/1‑131 (MQ=255)
                                                        cAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCt                                                                                                     <  2:83865/130‑1 (MQ=255)
                                                         aaGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                  <  1:31921/148‑1 (MQ=255)
                                                                gCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                                          <  5:31305/149‑1 (MQ=255)
                                                                     ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                                     >  7:10794/1‑149 (MQ=255)
                                                                            agaCGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCCCAAGCCAGCAAccc                                                               >  7:249032/1‑148 (MQ=255)
                                                                                gCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCACAtct                                                          >  8:78249/1‑149 (MQ=255)
                                                                                  cAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCca                                                                             <  7:130322/128‑1 (MQ=255)
                                                                                  cAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCca                                                                             >  8:130322/1‑128 (MQ=255)
                                                                                    gCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                                      <  2:326960/149‑1 (MQ=255)
                                                                                        aTGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAg                                                   >  6:51843/1‑148 (MQ=255)
                                                                                           aaGTCTGGCGGCTCGCCGAGTTCTAGCACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAcc                                               <  8:178482/149‑1 (MQ=255)
                                                                                                         cGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCAc                                    <  2:362957/147‑1 (MQ=255)
                                                                                                         cGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACta                                  >  1:362957/1‑148 (MQ=255)
                                                                                                          gCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                                 >  3:265031/1‑149 (MQ=255)
                                                                                                           ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAc                                                <  4:247316/133‑1 (MQ=255)
                                                                                                           ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAc                                                >  3:247316/1‑133 (MQ=255)
                                                                                                           ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                                 >  3:302936/1‑148 (MQ=255)
                                                                                                            cGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCa                               <  5:51843/149‑1 (MQ=255)
                                                                                                              aGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGc                             <  7:249122/149‑1 (MQ=255)
                                                                                                                ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTACGGCCGCATCAACATTACACAAGCCAGCAACCCCAAAACCAACAGCACCCACATACGCCCAGCCAGCa                            >  2:377859/1‑148 (MQ=255)
                                                                                                                   tAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                  <  7:264361/91‑1 (MQ=255)
                                                                                                                   tAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                  >  8:264361/1‑91 (MQ=255)
                                                                                                                      aaCGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTg                     <  4:265031/149‑1 (MQ=255)
                                                                                                                        cGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGc                             <  1:95312/139‑1 (MQ=255)
                                                                                                                        cGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGc                             >  2:95312/1‑139 (MQ=255)
                                                                                                                           aTCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCTGGCCGCTTCAACTTTCCACAAGCCAGCAACCCCAACTCTAACAGCACCAACCTCGGCACTGACAGCATCGTCTGCGaa                 >  7:302259/1‑148 (MQ=255)
                                                                                                                            tCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa               <  8:10794/149‑1 (MQ=255)
                                                                                                                            tCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa               <  4:302936/149‑1 (MQ=255)
                                                                                                                              aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAac              <  7:78249/148‑1 (MQ=255)
                                                                                                                                   ctttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATc         <  8:150695/146‑1 (MQ=255)
                                                                                                                                   ctttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATc         >  7:150695/3‑148 (MQ=255)
                                                                                                                                       cgcgGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcggc     <  3:245728/148‑1 (MQ=255)
                                                                                                                                         cgGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGgt  <  1:377859/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCGGT  >  NC_002947/1188894‑1189178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: