Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,905 G→T 6.5% intergenic (‑51/+28) pqqB ← / ← pqqA coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947459,9050GT 6.5% 277.0 / 7.4 92intergenic (‑51/+28)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base T (2/4);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

CAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459759‑460053
                                                                                                                                                  |                                                                                                                                                    
cAGGGTGCCATCTCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAgg                                                                                                                                                     <  3:15268/148‑1 (MQ=255)
  gcgagCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACAGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGAAGAAACCGGATCGAAAAAACGAAAACCCCCCGGCAGGCCg                                                                                                                                                  >  8:139710/5‑149 (MQ=255)
    gTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAg                                                                                                                                                >  2:250272/1‑149 (MQ=255)
      gCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGc                                                                                                                                              <  5:264304/149‑1 (MQ=255)
         aTCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCAGGCAGGCCGAGGCGt                                                                                                                                            >  3:40627/1‑148 (MQ=255)
            aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACCAAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                                        <  7:300560/149‑5 (MQ=255)
             cGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAAc                                                                                                                                       <  1:310777/149‑1 (MQ=255)
              ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttttcgt                                                                                                                                      <  2:260005/149‑7 (MQ=255)
              ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAAcc                                                                                                                                      <  4:64763/149‑1 (MQ=255)
                   cccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGGACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAg                                                                                                                                 >  1:331180/1‑149 (MQ=255)
                     cTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGc                                                                                                                                                    >  5:250722/1‑128 (MQ=255)
                     cTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGc                                                                                                                                                    <  6:250722/128‑1 (MQ=255)
                     cTTGCAGTTGACGCAGTTGCAGTTCCACTGCGCGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGCATCGAAAAAACGAAAACGCCCCCGCAGGCCGACGCGTGAAAACCAAAGcc                                                                                                                               >  6:37124/1‑149 (MQ=255)
                        gCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                            >  6:40336/1‑149 (MQ=255)
                         ctgTTGACGCAGTTGCAGGTCCACTGCGGGAACCCGGCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCTGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGGCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                            <  2:278686/146‑1 (MQ=255)
                                cGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGc                                                                                                                           >  4:137698/1‑142 (MQ=255)
                                cGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGc                                                                                                                           <  3:137698/142‑1 (MQ=255)
                                  cagttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGtt                                                                                                                  <  5:40336/149‑1 (MQ=255)
                                   cgttgcaggtCCTCTGCGGGAACCCGCCACCGTCGGCGTAGCCGAGAACCTGGAGTTACATGGCCACTCTCTATGCAGAAACCTGATCCAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  8:121623/148‑1 (MQ=255)
                                   agttgcagttCCACTGCGGGGACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  2:331180/149‑1 (MQ=255)
                                   agttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGt                                                                                                                   >  1:14939/1‑147 (MQ=255)
                                   agttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  6:323356/149‑1 (MQ=255)
                                   agttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  7:248299/149‑1 (MQ=255)
                                    gttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAg                                                                                                                >  8:68430/1‑149 (MQ=255)
                                     ttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGc                                                                                                               <  1:250272/149‑1 (MQ=255)
                                          gttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 >  3:170125/1‑142 (MQ=255)
                                          gttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGt                                                                                                          <  6:52179/149‑1 (MQ=255)
                                          gttCCACTGCGGCAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  4:170125/142‑1 (MQ=255)
                                           ttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGCCAGGCCCAGGCGAGAAACCCTAAGCCTGGCTTAGAGGTTAGCGAAGTa                                                                                                         >  1:225576/1‑149 (MQ=255)
                                              cACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACa                                                                                                       <  7:68430/148‑1 (MQ=255)
                                              cACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACAt                                                                                                      >  4:202070/1‑149 (MQ=255)
                                               aCTTCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCCGATCGAAAAAACGAAAACGCCCCGACAGGCCGAGGCGTGAAACCCTAAGCCTGGCTAACCGGTTAGCGAAGTACATg                                                                                                     >  8:327483/1‑149 (MQ=255)
                                                cTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATgg                                                                                                    >  8:96195/1‑149 (MQ=255)
                                                   cGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                                   <  3:128006/147‑1 (MQ=255)
                                                   cGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                                   <  1:29522/147‑1 (MQ=255)
                                                     ggAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTACGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACt                                                                                               >  5:356236/1‑149 (MQ=255)
                                                     ggAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACt                                                                                               >  8:123958/1‑149 (MQ=255)
                                                       aaCCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTc                                                                                             <  2:12343/149‑1 (MQ=255)
                                                         cccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                            >  6:133854/1‑148 (MQ=255)
                                                          ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                            >  6:211592/1‑147 (MQ=255)
                                                           cGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCTCTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGaa                                                                                                            >  3:261710/1‑130 (MQ=255)
                                                           cGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTATCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                         <  4:40627/149‑1 (MQ=255)
                                                           cGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGaa                                                                                                            <  4:261710/130‑1 (MQ=255)
                                                           cGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACGCTCGCTGCAGAGGCCGGAGCGAAAAAACGAAAACGCCCCGGCAGGCCGCGGCGTGCACCCCTAAGCCTGGCGTAGCGGTTTGCGAAGTACATGGGGACTTCGacg                                                                                         >  5:374924/1‑147 (MQ=255)
                                                                           gCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAg                                                                         >  8:36465/1‑149 (MQ=255)
                                                                            ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAACACCAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                        <  2:361435/149‑1 (MQ=255)
                                                                            ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCGGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  1:286331/108‑1 (MQ=255)
                                                                            ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCGGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 >  2:286331/1‑108 (MQ=255)
                                                                            ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGAACGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                        <  8:4729/149‑1 (MQ=255)
                                                                                gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTAt                                                                    <  3:202070/149‑1 (MQ=255)
                                                                                gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGAACGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTCGCTTAGCGGTTAGCGAAGTACATCGTGACTTCGAAGCCGATACTCAGCTAAGTCTAt                                                                    >  1:178806/1‑149 (MQ=255)
                                                                                   ccTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                         >  1:364819/1‑125 (MQ=255)
                                                                                   ccTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                         <  2:364819/125‑1 (MQ=255)
                                                                                    gtGGATGAAAATGGCCACTCTGTATGCAGAAACGGGATCGAAAAAACCAAAACGCCCCGGCAGGCCGATGCGTGAAAGCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                                <  2:178806/148‑1 (MQ=255)
                                                                                    cTGGATGGACCTGGCCACGCTCTAGCCAGCACCCGCATCGAAACAACGAAACCGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCGTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGCCAGTGTATGCAg                                                                <  6:374924/149‑1 (MQ=255)
                                                                                       gATGTACATTGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCGGGCTGGCCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                            <  3:303988/86‑1 (MQ=255)
                                                                                       gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGTCGTGAAACCCTAAGCCTgg                                                                                                                            >  4:303988/1‑86 (MQ=255)
                                                                                            aCATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                                        >  4:175241/1‑149 (MQ=255)
                                                                                             cATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGAGTTCGAAGCCGATACGCAGGTCAGTCTATGCAGGTTTGGTc                                                        >  4:41634/1‑148 (MQ=255)
                                                                                             cATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                                        >  7:241433/1‑148 (MQ=255)
                                                                                              aTGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAgg                                                                                                                                                     >  1:214971/1‑54 (MQ=255)
                                                                                              aTGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAgg                                                                                                                                                     <  2:214971/54‑1 (MQ=255)
                                                                                              aTGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGAGTATGCAGGTTTGGTcc                                                       >  7:42431/1‑148 (MQ=255)
                                                                                                 gCCTCTCTCTATGCAGAAACCGGATCGAAAAAACGAAAAGGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGGCTGGCTTAGCGTTTAGCCAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                                    <  5:37124/148‑1 (MQ=255)
                                                                                                 gCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGAGTATGCAGGTTTGGTCCACAt                                                   <  8:42431/149‑1 (MQ=255)
                                                                                                   cACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATgg                                                 <  2:14939/149‑1 (MQ=255)
                                                                                                   cACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATg                                                  >  2:71707/1‑148 (MQ=255)
                                                                                                        tctATGGAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTa                                            >  6:4546/1‑149 (MQ=255)
                                                                                                         ctATGCAGAAACCGGATCCAAAAAACGAAAACGCACCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAACCCGATACGCAGGTCAGTGTATGCAGGTTAGGTCCAGATGCAATTa                                            >  8:223805/1‑148 (MQ=255)
                                                                                                             gCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTc                                         <  7:305099/147‑1 (MQ=255)
                                                                                                             gCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTc                                         >  8:305099/1‑147 (MQ=255)
                                                                                                               aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                                     >  5:125771/1‑149 (MQ=255)
                                                                                                               aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTCCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                                     >  1:13999/1‑149 (MQ=255)
                                                                                                                gAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTcc                                    <  7:96195/149‑1 (MQ=255)
                                                                                                                   aCCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGaa                                 <  3:175241/149‑1 (MQ=255)
                                                                                                                    ccGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAAt                                >  4:294494/1‑149 (MQ=255)
                                                                                                                      ggATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGa                              <  6:356236/149‑1 (MQ=255)
                                                                                                                           gAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggtt                         <  5:211592/149‑1 (MQ=255)
                                                                                                                            aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg                        >  3:212922/1‑149 (MQ=255)
                                                                                                                            aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg                        >  5:198270/1‑149 (MQ=255)
                                                                                                                              aaaacgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg                      >  7:21340/1‑149 (MQ=255)
                                                                                                                                aacgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTTTGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggtt                    <  7:139710/149‑1 (MQ=255)
                                                                                                                                   gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGgttgtt                 >  2:51051/1‑149 (MQ=255)
                                                                                                                                   gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGgttgtt                 >  5:327482/1‑149 (MQ=255)
                                                                                                                                    taaacgCCCCGGCAGGCCGAGGCGTTAAACCCTAAGCCGGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAATCCTATACGCAGGTCAGTGTATGCAGGTTTGGTCCTCATGGAATTACTCCTTCCGAATGAGTGTTGGGTTGTTc                <  6:337395/148‑1 (MQ=255)
                                                                                                                                     aaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGTAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc                <  6:134815/148‑1 (MQ=255)
                                                                                                                                       acgCCCCTGCAGGCCGAGGCGTGAAACCCTAAGCCTGGGTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc             <  7:327483/149‑1 (MQ=255)
                                                                                                                                        cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCTCATGGAATTACTCCTTCCGAATGAgggttgggt                     <  5:80203/140‑1 (MQ=255)
                                                                                                                                        cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggt                     <  5:85533/140‑1 (MQ=255)
                                                                                                                                        cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggt                     >  6:85533/1‑140 (MQ=255)
                                                                                                                                        cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggt                     >  6:80203/1‑140 (MQ=255)
                                                                                                                                         gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc             >  3:212650/1‑147 (MQ=255)
                                                                                                                                           cccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGctac          >  8:184517/1‑148 (MQ=255)
                                                                                                                                               tcAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTAGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTTGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACtt     <  3:41634/148‑1 (MQ=255)
                                                                                                                                                cAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTg    <  6:327482/149‑1 (MQ=255)
                                                                                                                                                  ggCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTggg  >  6:40301/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459759‑460053

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