Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,350,678 Δ1 bp 15.6% intergenic (‑219/+11) elbB ← / ← arcB isoprenoid biosynthesis protein with amidotransferase‑like domain/aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,6780T.15.6% 356.4 / 62.9 128intergenic (‑219/+11)elbB/arcBisoprenoid biosynthesis protein with amidotransferase‑like domain/aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base T (48/59);  new base . (10/10);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 128 total aligned reads displayed.

CGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACC  >  NC_000913/3350580‑3350775
                                                                                                  |                                                                                                  
cgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTc                                                                                                   <  2:2236279/100‑1 (MQ=255)
   ccTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAg                                                                                                >  2:1967469/1‑100 (MQ=255)
      gAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCCGTATTCGCGCACCCCGGTCTAG‑cc                                                                                             >  1:2375849/1‑99 (MQ=255)
      gAAAAACAGTGCTGTGACCATGTAACTCATCATAACCATTTACCGCGCAGCCAAGATTTCCCTGGTGATGACGCACAATTCGCGCACGCCGGTCTAG‑cc                                                                                             >  1:831897/1‑99 (MQ=255)
       aaaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCgg                                                                                           <  2:1568690/100‑1 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCggg                                                                                          >  1:1936085/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCggg                                                                                          >  1:1353601/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCggg                                                                                          >  1:408771/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCggg                                                                                          >  1:426085/1‑100 (MQ=255)
           aCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtg                                                                                       <  1:1159743/100‑2 (MQ=255)
            cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgc                                                                                      <  2:1864915/100‑3 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                     <  2:444466/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                     <  1:2532105/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                     <  2:1591413/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCAGCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGGTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                     <  1:2464555/100‑4 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtt                                                                                    >  2:852515/1‑100 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtt                                                                                    <  1:156224/100‑1 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcgc                                                                                    <  2:2078070/100‑5 (MQ=255)
               tGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttt                                                                                   <  1:934766/100‑1 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                  >  2:1249635/1‑100 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                  >  1:1279553/1‑100 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                  <  2:299032/100‑1 (MQ=255)
                 cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcgcgaa                                                                                 <  1:1636760/100‑8 (MQ=38)
                 cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttttt                                                                                 <  1:2062265/100‑1 (MQ=255)
                    tgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAg                                                                              <  1:1064918/100‑1 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTg                                                                            <  2:2350329/100‑1 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTg                                                                            >  2:2555464/1‑100 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTg                                                                            >  2:2555823/1‑100 (MQ=255)
                        cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGc                                                                          <  2:1443073/100‑1 (MQ=255)
                        cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGc                                                                          <  2:1409213/100‑1 (MQ=255)
                          tGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCtt                                                                        >  2:2211467/1‑100 (MQ=255)
                             aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTg                                                                     <  1:2037564/100‑1 (MQ=255)
                             aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTg                                                                     >  2:1605247/1‑100 (MQ=255)
                              aCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGc                                                                    <  1:2319515/100‑1 (MQ=255)
                               ctaatcatAATAATTGACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCATTATTCGCGCATCCCGGTCTAG‑CCGGGGTCATTTTTTAGt                                                                             <  1:899670/87‑1 (MQ=255)
                               cTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTCTTAGTGGCTTTTGcc                                                                   >  2:2549995/1‑100 (MQ=255)
                                   tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                               >  1:2012967/1‑100 (MQ=255)
                                   tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                               >  1:1303153/1‑100 (MQ=255)
                                    catAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccaccca                                                              >  2:1692560/1‑100 (MQ=255)
                                     ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             <  1:139885/100‑1 (MQ=255)
                                     ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             >  1:240426/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  1:189616/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            <  2:575862/100‑1 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  2:240091/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  2:1473462/1‑100 (MQ=255)
                                       aataatTTACGGCGCAGCAAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             <  1:2295812/98‑1 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       >  1:1820424/1‑100 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       <  1:95065/100‑1 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       <  1:2552649/100‑1 (MQ=255)
                                            tttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGTGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                      <  1:174449/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1462633/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:2138472/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1844508/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1844507/100‑1 (MQ=255)
                                              tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCCCCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                    <  2:1674940/100‑1 (MQ=255)
                                               aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                   <  2:605386/100‑1 (MQ=255)
                                               aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                   >  1:87304/1‑100 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTCG‑CCGGGGTCATCTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 <  1:1808454/100‑1 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 >  2:46548/1‑100 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 <  1:424559/100‑1 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  2:1311246/1‑100 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  1:1129631/1‑100 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  1:2216510/1‑100 (MQ=255)
                                                   cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  2:1528489/1‑99 (MQ=255)
                                                       gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg                                           >  2:446642/1‑100 (MQ=255)
                                                         cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                         >  2:2251178/1‑100 (MQ=255)
                                                         cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                         >  1:1982812/1‑100 (MQ=255)
                                                          aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                        <  2:1226621/100‑1 (MQ=255)
                                                          aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                        <  2:2077678/100‑1 (MQ=255)
                                                              tttCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGAc                                    <  2:60601/100‑1 (MQ=255)
                                                                     ggTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCActc                             >  1:837064/1‑100 (MQ=255)
                                                                     ggTGTGGGCGCAGTATTCGCGCACCCCGGTGTGT‑CCGGGGTCATTTTTTAGTGGCTTGTGCCGCCCACGCTTTCAGCACTTCTACGTCGTGACGCCAc                               <  1:1985441/98‑1 (MQ=255)
                                                                          tGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTcttctt                        >  2:1774638/1‑100 (MQ=255)
                                                                           ggCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCttt                       <  1:2591319/100‑1 (MQ=255)
                                                                            gcgcAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTc                      >  2:777044/1‑100 (MQ=255)
                                                                             cgcAGTATTCGCGCATCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCa                     <  1:1009869/100‑1 (MQ=255)
                                                                             cgcAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCa                     <  1:1080903/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  2:527219/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  1:1914033/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  2:335749/100‑1 (MQ=255)
                                                                                 ctaaaaaatgaCCCCGG‑CTAGACCGGGGTCATTTTTTATTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                 >  1:1759189/11‑100 (MQ=37)
                                                                                 ctaaaaaatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                 >  2:2153342/11‑100 (MQ=37)
                                                                                 ataaaaaatgaCCCCGG‑CTAGACCGGGGCCATTTGTTAGTGGCTTTTGCCACCCACGCTTTCAGCAATTATACGTCGTAACGCCACTCTTCTTTCAtctc                 >  1:246713/11‑100 (MQ=25)
                                                                                  taaaaaatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATctct                >  1:1349059/10‑100 (MQ=38)
                                                                                   aaaaaatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt               >  2:1837895/9‑100 (MQ=38)
                                                                                   aTTCGCGCACCACGGGCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACCCTTCTTTCAtct                  >  1:560909/1‑97 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:822654/100‑1 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:1544794/100‑1 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:333054/100‑1 (MQ=255)
                                                                                      cgcgCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCaa            >  2:1610667/1‑100 (MQ=255)
                                                                                       gcgcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           <  2:2194044/100‑1 (MQ=255)
                                                                                       aatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           >  2:1392731/5‑100 (MQ=255)
                                                                                       aatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           >  1:2093232/5‑100 (MQ=255)
                                                                                        cgcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATc          <  1:1473453/100‑1 (MQ=255)
                                                                                        atgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACCCTTTCAGCACTTCTCCGTCGTTACGCCACTCTTCTTTCATCTCTTCAGTc          >  1:211742/4‑100 (MQ=255)
                                                                                         tgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATcc         >  2:1076110/3‑100 (MQ=255)
                                                                                         gcACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGGCGCCACTCTTCTTTCATCTCTTCAATcc         >  1:1930542/1‑100 (MQ=255)
                                                                                          cACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCa        >  1:2476790/1‑100 (MQ=255)
                                                                                                ggTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa    <  1:1116645/98‑1 (MQ=255)
                                                                                                ggTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc  <  2:400422/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACC  >  NC_000913/3350580‑3350775

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