Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,893 T→A 5.6% intergenic (+28/‑55) benR → / → PP_3160 BenABC operon transcriptional activator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,580,8930TA 5.6% 279.0 / 2.6 84intergenic (+28/‑55)benR/PP_3160BenABC operon transcriptional activator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base A (3/2);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAGAAGCTTCGGC  >  NC_002947/3580750‑3581042
                                                                                                                                                |                                                                                                                                                     
gcAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCa                                                                                                                                                 >  2:1463697/1‑150 (MQ=255)
  aGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGt                                                                                                                                               >  2:377361/1‑150 (MQ=255)
   gtgtGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGtt                                                                                                                                              >  1:1988185/1‑150 (MQ=255)
   gtgtGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggt                                                                                                                                              <  1:2273342/150‑3 (MQ=255)
      tgACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCg                                                                                                                                           >  1:368128/1‑150 (MQ=255)
       gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCAttt                                                                                                                                                       <  1:33531/137‑1 (MQ=255)
       gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCAttt                                                                                                                                                       >  2:33531/1‑137 (MQ=255)
       gACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGa                                                                                                                                          <  1:316004/150‑1 (MQ=255)
               tGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTg                                                                                                                                  <  1:1410729/150‑1 (MQ=255)
               tGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTg                                                                                                                                  <  1:2070/150‑1 (MQ=255)
                 gCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑TTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCt                                                                                                                                <  1:432259/150‑1 (MQ=255)
                 gCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCGGGCATTTATGCCAGTTGCGAGCGCTTTGCt                                                                                                                                <  2:885282/150‑1 (MQ=255)
                  cccTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTg                                                                                                                               >  1:1788968/1‑150 (MQ=255)
                    cTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGtt                                                                                                                             >  1:293870/1‑150 (MQ=255)
                      ggATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTaa                                                                                                                           <  1:377361/150‑1 (MQ=255)
                            cGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAaccacc                                                                                                                     >  2:41825/1‑150 (MQ=255)
                             ggTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAAccacca                                                                                                                    >  1:2401267/1‑150 (MQ=255)
                               tttCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCAtt                                                                                                                  >  2:830561/1‑150 (MQ=255)
                                 tCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAg                                                                                                                <  2:2401267/150‑1 (MQ=255)
                                    ttCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTc                                                                                                             >  1:283181/1‑150 (MQ=255)
                                    ttCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTc                                                                                                             >  1:751164/1‑150 (MQ=255)
                                       aTCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTAt                                                                                                          >  2:2390564/1‑150 (MQ=255)
                                       aTCTGGGGCGGTTTTACGAGGTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTAt                                                                                                          >  2:2395190/1‑150 (MQ=255)
                                         cTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAac                                                                                                                         >  2:1367535/1‑133 (MQ=255)
                                         cTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAac                                                                                                                         <  1:1367535/133‑1 (MQ=255)
                                         cTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGc                                                                                                        >  1:959217/1‑150 (MQ=255)
                                           ggggCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCg                                                                                                      <  1:830561/150‑1 (MQ=255)
                                            gggCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGAACTGGCATTTATGCCa                                                                                                                                                 <  1:1574431/106‑1 (MQ=255)
                                            gggCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGAACTGGCATTTATGCCa                                                                                                                                                 >  2:1574431/1‑106 (MQ=255)
                                                 gTTTTCCGATATTTATCGGCAGCAGTTTGGCGAGTTTCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCTTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGaa                                                                                                <  1:1139320/150‑1 (MQ=255)
                                                    ttCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg                                                                                             <  1:41825/150‑1 (MQ=255)
                                                    ttCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATgg                                                                                             <  1:446847/150‑1 (MQ=255)
                                                       cGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggt                                                                                                                                              >  1:2031694/1‑96 (MQ=255)
                                                       cGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATAAATGCCAggt                                                                                                                                              <  2:2031694/98‑3 (MQ=255)
                                                       cGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGACTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTca                                                                                          >  2:196749/1‑150 (MQ=255)
                                                         agaTTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTcaca                                                                                        >  1:1513652/1‑150 (MQ=255)
                                                            tttATCGGCAGCAGTTTGGCGAGTTGCCTTTACCGCCCTT‑AACCGCATCAAGGCGACTTCGTTTGGCTTGTGGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGcc                                                                                     >  1:2223682/1‑150 (MQ=255)
                                                              tATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCAt                                                                                   >  1:1841038/1‑150 (MQ=255)
                                                                tCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt                                                                                                                              <  1:2384696/105‑1 (MQ=255)
                                                                tCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGt                                                                                                                              >  2:2384696/1‑105 (MQ=255)
                                                                tCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATgc                                                                                 <  2:1988185/150‑1 (MQ=255)
                                                                 cGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATgca                                                                                <  1:2390564/150‑1 (MQ=255)
                                                                   gcagcaGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTTCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAAt                                                                                               >  1:549871/1‑133 (MQ=255)
                                                                   gcagcaGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAAt                                                                                               <  2:549871/133‑1 (MQ=255)
                                                                   gcagcaGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATgcagc                                                                              >  2:740094/1‑150 (MQ=255)
                                                                    cagcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATgcagca                                                                             >  2:1548402/1‑150 (MQ=255)
                                                                    cagcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATgcagca                                                                             <  2:293870/150‑1 (MQ=255)
                                                                     agcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCAt                                                                                                                   <  2:633216/111‑1 (MQ=255)
                                                                     agcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCAt                                                                                                                   >  1:633216/1‑111 (MQ=255)
                                                                      gcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAg                                                                                                                <  2:1042424/113‑1 (MQ=255)
                                                                      gcagTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAg                                                                                                                >  1:1042424/1‑113 (MQ=255)
                                                                          tttGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCa                                                                       >  2:1515805/1‑150 (MQ=255)
                                                                                gAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACa                                                                 >  1:210832/1‑150 (MQ=255)
                                                                                   ttGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAat                                                              >  1:1033753/1‑150 (MQ=255)
                                                                                      ccTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAattta                                                           >  1:1779260/1‑150 (MQ=255)
                                                                                       cTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAatttat                                                          <  1:1515805/150‑1 (MQ=255)
                                                                                        ttCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAatttatt                                                         >  2:170206/1‑150 (MQ=255)
                                                                                           acaGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCt                                                      >  2:981837/1‑150 (MQ=255)
                                                                                           acaGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCt                                                      <  2:1301574/150‑1 (MQ=255)
                                                                                             aGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGc                                                                                      >  1:2343371/1‑116 (MQ=255)
                                                                                             aGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGc                                                                                      <  2:2343371/116‑1 (MQ=255)
                                                                                             aGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTAc                                                    <  2:283181/150‑1 (MQ=255)
                                                                                                 aTTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTc                                                <  2:368128/150‑1 (MQ=255)
                                                                                                  ttCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGGATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTcc                                               >  2:1311940/1‑150 (MQ=255)
                                                                                                   tCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTccc                                              <  2:1841038/150‑1 (MQ=255)
                                                                                                   tCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTccc                                              <  2:1788968/150‑1 (MQ=255)
                                                                                                    cAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGt                                                                           <  2:154078/120‑1 (MQ=255)
                                                                                                    cAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGt                                                                           >  1:154078/1‑120 (MQ=255)
                                                                                                    cAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCg                                             >  1:717051/1‑150 (MQ=255)
                                                                                                     aaGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGa                                            >  2:580797/1‑150 (MQ=255)
                                                                                                       gcggcgCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAg                                          <  1:170206/150‑1 (MQ=255)
                                                                                                        ccgatctc‑gTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAgg                                         <  2:2351525/143‑1 (MQ=255)
                                                                                                          gcgcgc‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGAc                                       >  2:2215014/1‑150 (MQ=255)
                                                                                                               c‑gTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCg                                  <  2:717051/150‑1 (MQ=255)
                                                                                                               c‑gTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGACAGTTGCGAACGCTTTGCTGATAACCAACCATTAGTTCTATGCCGGACGAATTGTTACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCg                                  >  2:275004/1‑150 (MQ=255)
                                                                                                                 gTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGc                                 <  1:1311940/150‑1 (MQ=255)
                                                                                                                  tGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCt                                <  2:1513652/150‑1 (MQ=255)
                                                                                                                  tGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCt                                <  2:959217/150‑1 (MQ=255)
                                                                                                                   gACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGAc                                                                                                   <  1:1214438/82‑1 (MQ=255)
                                                                                                                   gACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGAc                                                                                                   >  2:1214438/1‑82 (MQ=255)
                                                                                                                           ttGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAata                       >  2:109274/1‑150 (MQ=255)
                                                                                                                                ttGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTg                  >  1:1796746/1‑150 (MQ=255)
                                                                                                                                 tGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTgg                 <  2:2223682/150‑1 (MQ=255)
                                                                                                                                      aacTGGCATAAATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAg            >  2:403879/3‑150 (MQ=255)
                                                                                                                                             aTTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATAAGTTGGCGCAGAAGCTTc     >  1:568705/1‑150 (MQ=255)
                                                                                                                                             aTTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAAGATGTTGGCGCAGAAGCTTc     >  2:1587435/1‑150 (MQ=255)
                                                                                                                                                tATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTAGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAGAAGCTTCGGc  <  2:1180955/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAGTTTGGCGAGTTGCCTTCACAGACATTCAAGCGGCGCGC‑GTGACGGCGTTTGGCTTGTAGACCTGGCATTTATGCCAGTTGCGAGCGCTTTGCTGTTAAACCACCATTAGTTCTATGCCGGACGAATGGTCACAGCCATGCAGCAGTCGCAGAAACAAATTTATTTCTACGCTCCCGAAGGACTTTCGCTGGTCAAATATGTTGGCGCAGAAGCTTCGGC  >  NC_002947/3580750‑3581042

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