Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,899 C→T 6.4% intergenic (‑45/+34) 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,8990CT 6.4% 277.4 / 6.3 93intergenic (‑45/+34)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base T (2/4);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GTGGCTTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTAC  >  NC_002947/459751‑460045
                                                                                                                                                    |                                                                                                                                                  
gTGGCTTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGccc                                                                                                                                                    <  7:82198/149‑1 (MQ=255)
 tGGCTTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGcccc                                                                                                                                                   <  4:180326/149‑1 (MQ=255)
     tttCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCgg                                                                                                                                                 >  6:264304/1‑147 (MQ=255)
        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                                                                                                                   <  3:137698/142‑1 (MQ=255)
                                        cGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGc                                                                                                                   >  4:137698/1‑142 (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                                                                                 <  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)
                                                                                                      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                                 >  8:305099/1‑147 (MQ=255)
                                                                                                                     gCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTc                                 <  7:305099/147‑1 (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:80203/1‑140 (MQ=255)
                                                                                                                                                cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggt             >  6:85533/1‑140 (MQ=255)
                                                                                                                                                 gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc     >  3:212650/1‑147 (MQ=255)
                                                                                                                                                   cccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGctac  >  8:184517/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTGGCTTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTAC  >  NC_002947/459751‑460045

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: