Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,899 C→T 7.6% 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 7.6% 224.8 / 5.5 78intergenic (‑45/+34)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/4);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTC  >  NC_002947/459756‑460039
                                                                                                                                                      |                                                                                                                                            
tttCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCgg                                                                                                                                           >  5:325212/1‑147 (MQ=255)
tttCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCgg                                                                                                                                           >  2:262203/1‑147 (MQ=255)
            aTCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGtt                                                                                                                             >  2:57181/1‑148 (MQ=255)
             tCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGa                                                                                                                            <  6:235103/149‑1 (MQ=255)
               aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                          <  3:138863/149‑5 (MQ=255)
               aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGaaa                                                                                                                          <  5:87881/149‑1 (MQ=255)
               aCGGTAGCCCGTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTTGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                          <  8:239358/149‑5 (MQ=255)
                   tAGCCCTTGAAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCCCGACCCTACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGtttt                                                                                                                           <  1:115525/149‑4 (MQ=255)
                       ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                  <  7:280146/149‑1 (MQ=255)
                       ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                  >  5:60316/1‑149 (MQ=255)
                       ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGACGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                  >  1:75531/1‑149 (MQ=255)
                           gCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTg                                                                                                               <  5:345680/148‑1 (MQ=255)
                              gTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCt                                                                                                            >  6:82322/1‑148 (MQ=255)
                                 gACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCa                                                                                                                                         >  7:284317/1‑116 (MQ=255)
                                 gACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCa                                                                                                                                         <  8:284317/116‑1 (MQ=255)
                                  aCGCATTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                       <  5:51578/149‑1 (MQ=255)
                                  aCGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                       <  5:46487/149‑1 (MQ=255)
                                    gcagttgcagttCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGCACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGCGGCGTGCACCCCTAAGCCTGGCTTCGCGGt                                                                                                     >  4:101164/1‑149 (MQ=255)
                                     cagttgcagttCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGt                                                                                                     >  8:235000/1‑148 (MQ=255)
                                            agttCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAg                                                                                             <  7:235000/149‑1 (MQ=255)
                                             gttCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGt                                                                                            <  6:60316/149‑1 (MQ=255)
                                                ccACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTAc                                                                                          <  2:183032/148‑1 (MQ=255)
                                                  aCTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGAGTGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                   >  1:240422/1‑138 (MQ=255)
                                                  aCTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                   <  2:240422/137‑1 (MQ=255)
                                                  aCTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                       >  2:115742/1‑149 (MQ=255)
                                                      cGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGa                                                                                   >  5:315830/1‑149 (MQ=255)
                                                      cGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGAACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                     <  1:115742/147‑1 (MQ=255)
                                                       gggAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGCTCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGAc                                                                                  <  6:262528/149‑1 (MQ=255)
                                                               ccG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                              >  6:262099/1‑147 (MQ=255)
                                                               ccG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGaa                                                                            <  2:41787/149‑1 (MQ=255)
                                                                cG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGGGGTTAGCGAAGTACAGGGTGACGTCGAAg                                                                           >  4:251389/1‑149 (MQ=255)
                                                                    ccACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGcc                                                                         >  8:199974/1‑149 (MQ=255)
                                                                        ccGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACTGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGc                                                                          <  8:86763/145‑1 (MQ=255)
                                                                        ccGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACGGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGc                                                                          >  7:86763/1‑145 (MQ=255)
                                                                          ggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATAc                                                                    >  2:116988/1‑149 (MQ=255)
                                                                                ggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGt                                                              >  3:64470/1‑149 (MQ=255)
                                                                                     cGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCa                                                            >  5:7603/1‑146 (MQ=255)
                                                                                     cGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCa                                                            <  6:7603/146‑1 (MQ=255)
                                                                                      gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                       >  5:145542/1‑118 (MQ=255)
                                                                                      gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                       <  6:145542/118‑1 (MQ=255)
                                                                                      gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                        <  1:116988/149‑1 (MQ=255)
                                                                                            cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                  <  5:82322/149‑1 (MQ=255)
                                                                                            cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                  >  5:85230/1‑149 (MQ=255)
                                                                                            cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                  <  5:262099/149‑1 (MQ=255)
                                                                                               gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                        <  6:77128/140‑1 (MQ=255)
                                                                                               gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                        >  5:77128/1‑140 (MQ=255)
                                                                                               gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGtt                                               <  3:251389/149‑1 (MQ=255)
                                                                                                aTGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGttt                                              <  6:325212/149‑1 (MQ=255)
                                                                                                  tGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTg                                             >  6:11711/1‑149 (MQ=255)
                                                                                                  tGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTg                                             >  6:16599/1‑149 (MQ=255)
                                                                                                     aCATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGTAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGAt                                                                      <  2:82875/121‑1 (MQ=255)
                                                                                                     aCATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGTAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGAt                                                                      >  1:82875/1‑121 (MQ=255)
                                                                                                       aTGGCCACTCTCTATGCAGAAACCGGATCGAAAAATACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                          >  1:225294/1‑132 (MQ=255)
                                                                                                       aTGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                          <  2:225294/131‑1 (MQ=255)
                                                                                                        tGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCca                                        >  5:250347/1‑148 (MQ=255)
                                                                                                         ggCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGaa                                                                            <  6:50326/111‑1 (MQ=255)
                                                                                                         ggCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGaa                                                                            >  5:50326/1‑111 (MQ=255)
                                                                                                                   tATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACt                            >  3:315182/1‑149 (MQ=255)
                                                                                                                     tGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGTCTTCGAAGCCGATACGCAGGTCTGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                          <  6:57804/149‑1 (MQ=255)
                                                                                                                     tGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                          >  7:106472/1‑149 (MQ=255)
                                                                                                                     tGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                          <  7:199974/149‑1 (MQ=255)
                                                                                                                        aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTATCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                       <  3:101164/149‑1 (MQ=255)
                                                                                                                        aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                       <  4:64470/149‑1 (MQ=255)
                                                                                                                        aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                       >  6:72336/1‑149 (MQ=255)
                                                                                                                         gAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGATGCCGATACGCAGGTCAGTGTATGCAGGTTTTGTCCACATGGAATTACTCCTTcc                      <  1:262203/149‑1 (MQ=255)
                                                                                                                         gAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTcc                      <  6:250347/149‑1 (MQ=255)
                                                                                                                               ggATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAAt                  >  1:153776/1‑147 (MQ=255)
                                                                                                                                gATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGt                                                                                            >  4:175273/1‑72 (MQ=255)
                                                                                                                                gATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGt                                                                                            <  3:175273/72‑1 (MQ=255)
                                                                                                                                     aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg          >  2:6727/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg          >  5:9873/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg          >  5:254988/1‑149 (MQ=255)
                                                                                                                                       aaaacgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg        <  6:85230/149‑1 (MQ=255)
                                                                                                                                              ttcacgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCt                         <  4:26204/123‑1 (MQ=255)
                                                                                                                                              ttcacgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCt                         >  3:26204/4‑126 (MQ=255)
                                                                                                                                              aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc  <  4:315182/149‑1 (MQ=255)
                                                                                                                                                   gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg              >  7:158273/1‑132 (MQ=255)
                                                                                                                                                   gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg              <  8:158273/132‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
TTTCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTC  >  NC_002947/459756‑460039

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