Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,921,923 T→A 6.9% intergenic (‑114/+26) acs ← / ← peaD acyl‑CoA synthetase/quinohemoprotein amine dehydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,921,9230TA 6.9% 339.0 / 15.3 116intergenic (‑114/+26)acs/peaDacyl‑CoA synthetase/quinohemoprotein amine dehydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (54/53);  new base A (3/5);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 120 total aligned reads displayed.

GTCCTGGGTGCCCTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTGCCGCCCAGGTA  >  NC_002947/3921774‑3922071
                                                                                                                                                     |                                                                                                                                                    
gTCCTGGGTGCCCTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGtt                                                                                                                                                      >  2:2801166/1‑150 (MQ=255)
         gCCCTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTATGTGCTGGCCTCTTCTCGTGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCga                                                                                                                                             >  1:1149547/1‑150 (MQ=255)
          cccTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgcc                                                                                                                                            <  2:232561/150‑3 (MQ=255)
            cTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCgaaga                                                                                                                                          >  1:1904991/1‑150 (MQ=255)
             tGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCgaagaa                                                                                                                                         <  1:1575059/150‑1 (MQ=255)
             tGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgccgac                                                                                                                                         <  1:1874038/150‑6 (MQ=255)
                gtgtAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGTGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGca                                                                                                                                      >  2:777244/1‑150 (MQ=255)
                  gtAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGcaca                                                                                                                                    <  2:2318539/150‑1 (MQ=255)
                   tAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGcacac                                                                                                                                   >  2:2240470/1‑150 (MQ=255)
                       ttGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGGGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGAGGAGGATTACGTCGGCGCAGGTTTACTCGCGGAGAAGGACGCACCg                                                                                                                               >  1:1445793/1‑150 (MQ=255)
                            tGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTa                                                                                                                          >  1:787216/1‑150 (MQ=255)
                            tGGGACATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTa                                                                                                                          <  2:1701220/150‑1 (MQ=255)
                             ggggCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGGGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAg                                                                                                                         <  1:997956/150‑1 (MQ=255)
                                gCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGCAGAAGGCGCCACATCTGGGGCCTGTGGAGGCTTTCGTCGGCCCGGGTTTACCCGCGAAGAAGCACACACCGACTTGGcga                                                                                                                      >  2:2099194/1‑150 (MQ=38)
                                gCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGcga                                                                                                                      >  2:1892054/1‑150 (MQ=255)
                                  aTGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGcgaa                                                                                                                     <  2:2712889/149‑1 (MQ=255)
                                  aTGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGcgaa                                                                                                                     >  1:2712889/1‑149 (MQ=255)
                                     ggCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGcgaacgaa                                                                                                                 <  2:2673400/150‑1 (MQ=255)
                                      gCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTCCACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCgaacgaac                                                                                                                >  1:2090919/1‑150 (MQ=255)
                                          tctcTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCt                                                                                                            <  1:315707/150‑1 (MQ=255)
                                            tctTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTgc                                                                                                          <  2:1220157/150‑1 (MQ=255)
                                              ttGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTgcgg                                                                                                        <  1:2753825/150‑1 (MQ=255)
                                                  tGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGAGAAGACGAACACTCCGGTTTAGCGAACGAACACCTGCGGcgtc                                                                                                    >  2:2254116/1‑150 (MQ=255)
                                                  tGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCACGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtc                                                                                                    >  2:2219562/1‑150 (MQ=255)
                                                   gTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCgaagaa                                                                                                                                         <  1:1742262/112‑1 (MQ=255)
                                                   gTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCgaagaa                                                                                                                                         >  2:1742262/1‑112 (MQ=255)
                                                   gTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtcg                                                                                                   >  2:3020274/1‑150 (MQ=255)
                                                    ttCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtcgt                                                                                                  <  1:2240470/150‑1 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTAccc                                                                                                                                                 <  2:2502564/100‑1 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTAccc                                                                                                                                                 >  1:2502564/1‑100 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtcgt                                                                                                  >  1:1403830/1‑147 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtcgt                                                                                                  <  2:1403830/147‑1 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGt                                                                                               >  1:913043/1‑150 (MQ=255)
                                                       ttGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGt                                                                                               >  2:20477/1‑150 (MQ=255)
                                                        tgtgTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTa                                                                                              <  1:1179838/150‑1 (MQ=255)
                                                          tgtTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGa                                                                                            <  1:399209/150‑1 (MQ=255)
                                                           gtTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGAc                                                                                           <  1:2254116/150‑1 (MQ=255)
                                                             tttCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACAt                                                                                         <  1:1892054/150‑1 (MQ=255)
                                                             tttCTGTGCTGGCCTCTTCGCGGGTGAACCCCCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACAt                                                                                         <  1:2099194/150‑1 (MQ=255)
                                                                       ggCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGa                                                                                            <  2:1158563/137‑1 (MQ=255)
                                                                       ggCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGa                                                                                            >  1:1158563/1‑137 (MQ=255)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTAcc                                                                                                                                                  >  1:2664293/1‑82 (MQ=255)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTAcc                                                                                                                                                  <  2:2664293/82‑1 (MQ=255)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCgaag                                                                                                                                           <  2:203145/89‑1 (MQ=255)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTGTACCCGCgaag                                                                                                                                           >  1:203145/1‑89 (MQ=255)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgccga                                                                                                                                          <  1:1540506/90‑5 (MQ=37)
                                                                        gCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgccga                                                                                                                                          >  2:1540506/1‑86 (MQ=37)
                                                                          ctctTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGc                                                                            >  1:1293497/1‑150 (MQ=255)
                                                                           tctTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCa                                                                           >  2:2368795/1‑150 (MQ=255)
                                                                                 cgGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgc                                                                                                                                             >  1:1780642/1‑77 (MQ=38)
                                                                                 cgGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGCGCCTGTGGAGGTTTTCGTCGGCGCGGGTAAACCcgcgc                                                                                                                                             <  2:1780642/78‑2 (MQ=38)
                                                                                 cgGGTGAACACGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTCTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGa                                                                     <  1:1079480/150‑1 (MQ=255)
                                                                                  gggTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGAt                                                                    <  1:20477/150‑1 (MQ=255)
                                                                                        aCCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGttctt                                                              >  2:2752383/1‑150 (MQ=255)
                                                                                         cccGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTca                                                                                     <  1:525516/126‑1 (MQ=255)
                                                                                         cccGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTca                                                                                     >  2:525516/1‑126 (MQ=255)
                                                                                          ccGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCa                                                            >  2:2732750/1‑150 (MQ=255)
                                                                                          ccGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCa                                                            <  1:2561344/150‑1 (MQ=255)
                                                                                            gCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTa                                                                                              >  1:1938975/1‑114 (MQ=255)
                                                                                            gCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTa                                                                                              <  2:1938975/114‑1 (MQ=255)
                                                                                            gCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCAcc                                                          <  2:322629/150‑1 (MQ=255)
                                                                                            gCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCAcc                                                          >  2:1758916/1‑150 (MQ=255)
                                                                                                cTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCtttt                                                      >  2:2672522/1‑150 (MQ=255)
                                                                                                 tACACAGGAGGCGCCAACCTTGGGGACTGTGGAGGCTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGAACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTc                                                     >  2:2163494/1‑150 (MQ=255)
                                                                                                  acacAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtc                                                                                                    >  2:1472997/1‑102 (MQ=255)
                                                                                                  acacAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtc                                                                                                    <  1:1472997/102‑1 (MQ=255)
                                                                                                  acacAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGcgtc                                                                                                    >  2:1473363/1‑102 (MQ=255)
                                                                                                   cacaGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCa                                                   <  2:1904991/150‑1 (MQ=255)
                                                                                                     aaggaggCTGCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACATCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGc                                                 <  2:1445793/149‑1 (MQ=255)
                                                                                                        gaggCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTa                                              <  1:1098428/150‑1 (MQ=255)
                                                                                                           gcgcCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAGGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCg                                                >  1:219310/1‑145 (MQ=255)
                                                                                                             gcCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGAt                                                                    >  2:2265904/1‑123 (MQ=255)
                                                                                                             gcCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCgg                                         >  1:1761357/1‑150 (MQ=255)
                                                                                                               cacaTTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCAc                                                           <  1:1221694/130‑1 (MQ=255)
                                                                                                               cacaTTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCAc                                                           >  2:1221694/1‑130 (MQ=255)
                                                                                                                 aaTTTGGGGCCCGTGGATGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTg                                     <  2:2522380/149‑1 (MQ=255)
                                                                                                                  aTTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTcc                                                    <  2:2146217/134‑1 (MQ=255)
                                                                                                                  aTTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTcc                                                    >  1:2146217/1‑134 (MQ=255)
                                                                                                                  aTTTGGGGCCTGTGGAGGTTTCCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGa                                    >  2:74237/1‑150 (MQ=255)
                                                                                                                   tttGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGaa                                   >  2:261009/1‑150 (MQ=255)
                                                                                                                   tttGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGaa                                   <  1:2466818/150‑1 (MQ=255)
                                                                                                                   tttGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGaa                                   <  2:1293497/150‑1 (MQ=255)
                                                                                                                       gggCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACAcc                               >  2:1940266/1‑150 (MQ=255)
                                                                                                                            tgtgGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAg                          <  1:1758916/150‑1 (MQ=255)
                                                                                                                             gtgGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAgg                         <  1:1940266/150‑1 (MQ=255)
                                                                                                                                gAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCg                      >  2:2434589/1‑150 (MQ=255)
                                                                                                                                 aGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGATATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTAGCCGGGTTGAACACCGCCAGGTCGt                     <  2:1149547/150‑1 (MQ=255)
                                                                                                                                 aGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGt                     <  1:2368795/150‑1 (MQ=255)
                                                                                                                                   gTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTg                   <  1:2434589/150‑1 (MQ=255)
                                                                                                                                      ttCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAg                <  1:3020274/150‑1 (MQ=255)
                                                                                                                                      ttCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAg                <  1:855316/150‑1 (MQ=255)
                                                                                                                                        cgtcgGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGACAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGt              <  1:2719748/150‑1 (MQ=255)
                                                                                                                                         ttcttcgcgGGTAAACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGa                  <  1:3008131/140‑1 (MQ=255)
                                                                                                                                         ttcttcgcgGGTAAACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGa                  >  2:3008131/6‑145 (MQ=255)
                                                                                                                                          tcgGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCACCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTgc            <  2:1778479/150‑1 (MQ=255)
                                                                                                                                            ggCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGtt                                      <  2:559481/122‑1 (MQ=255)
                                                                                                                                                cgGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTGCCGCCCa      <  1:2752383/150‑1 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTGCCGCCCAGGt   >  1:983673/1‑150 (MQ=255)
                                                                                                                                                    tttACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGttc                                                                <  1:1808326/88‑1 (MQ=255)
                                                                                                                                                    tttACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTGCCGCCCAGGTa  >  2:2166042/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTCCTGGGTGCCCTGGGTGTAGCTTGGCTGGGGCATGGGCGTTCTCTTGCTGTTCTTGTGTTTTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACACAGGAGGCGCCACATTTGGGGCCTGTGGAGGTTTTCGTCGGCGCGGGTTTACCCGCGAAGAAGCACACACCGGTTTAGCGAACGAACACCTGCGGCGTCGTGGTAGACATGTCACCACCCGGCAACTTGATGTTCTTCACCTTTTCCAGCGTATCCGGGTTGAACACCGCCAGGTCGTTGAAGGTGCCGCCCAGGTA  >  NC_002947/3921774‑3922071

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