Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,910 A→G 6.7% intergenic (‑56/+23) 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,9100AG 6.7% 263.4 / 6.1 90intergenic (‑56/+23)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base G (2/4);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459763‑460053
                                                                                                                                                   |                                                                                                                                               
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)
                 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                                                                                                   <  1:29522/147‑1 (MQ=255)
                                               cGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                                   <  3:128006/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                                                                                                                 >  2:286331/1‑108 (MQ=255)
                                                                        ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCGGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                                 <  1:286331/108‑1 (MQ=255)
                                                                        ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGAACGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                        <  8:4729/149‑1 (MQ=255)
                                                                            gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTAt                                                                    <  3:202070/149‑1 (MQ=255)
                                                                            gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGAACGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTCGCTTAGCGGTTAGCGAAGTACATCGTGACTTCGAAGCCGATACTCAGCTAAGTCTAt                                                                    >  1:178806/1‑149 (MQ=255)
                                                                               ccTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                         <  2:364819/125‑1 (MQ=255)
                                                                               ccTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                         >  1:364819/1‑125 (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)
                                                                                          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)
                                                                                                                                                   |                                                                                                                                               
GTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459763‑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: ‑