Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,910 A→G 7.8% 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 7.8% 223.2 / 6.6 78intergenic (‑56/+23)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (2/4);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459768‑460053
                                                                                                                                                     |                                                                                                                                               
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)
                      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:262099/149‑1 (MQ=255)
                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                                >  5:85230/1‑149 (MQ=255)
                                                                                cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                                <  5:82322/149‑1 (MQ=255)
                                                                                   gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                                      >  5:77128/1‑140 (MQ=255)
                                                                                   gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                                      <  6:77128/140‑1 (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                                                                                          >  5:50326/1‑111 (MQ=255)
                                                                                             ggCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGaa                                                                                          <  6:50326/111‑1 (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                        >  5:9873/1‑149 (MQ=255)
                                                                                                                         aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg                        >  2:6727/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)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg   <  5:16599/149‑1 (MQ=255)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg   <  5:11711/149‑1 (MQ=255)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg   <  6:9873/149‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTggg  >  4:289124/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGGG  >  NC_002947/459768‑460053

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: