Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,904 A→C 8.1% intergenic (‑50/+29) 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,9040AC 8.1% 225.1 / 8.4 74intergenic (‑50/+29)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (2/4);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGG  >  NC_002947/459768‑460052
                                                                                                                                               |                                                                                                                                                    
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                                                                                                                               >  5:60316/1‑149 (MQ=255)
           ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                               <  7:280146/149‑1 (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                                                                         <  6:7603/146‑1 (MQ=255)
                                                                         cGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCa                                                                         >  5:7603/1‑146 (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                                                                                         >  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:199974/149‑1 (MQ=255)
                                                                                                         tGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                                       >  7:106472/1‑149 (MQ=255)
                                                                                                            aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTATCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                                    <  3:101164/149‑1 (MQ=255)
                                                                                                            aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                                    >  6:72336/1‑149 (MQ=255)
                                                                                                            aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTc                                    <  4:64470/149‑1 (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                                      >  3:26204/4‑126 (MQ=255)
                                                                                                                                  ttcacgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCt                                      <  4:26204/123‑1 (MQ=255)
                                                                                                                                  aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc               <  4:315182/149‑1 (MQ=255)
                                                                                                                                       gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg                           >  7:158273/1‑132 (MQ=255)
                                                                                                                                       gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg                           <  8:158273/132‑1 (MQ=255)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg  <  6:9873/149‑1 (MQ=255)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg  <  5:16599/149‑1 (MQ=255)
                                                                                                                                               aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg  <  5:11711/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAA‑‑CCCG‑‑CCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGG  >  NC_002947/459768‑460052

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: