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.

CGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACC  >  NC_000913/3350580‑3350775
                                                                                                  |                                                                                                  
cgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTc                                                                                                   <  2:2236279/100‑1 (MQ=255)
   ccTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAg                                                                                                >  2:1967469/1‑100 (MQ=255)
      gAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCCGTATTCGCGCACCCCGGTCTAGcc                                                                                             >  1:2375849/1‑99 (MQ=255)
      gAAAAACAGTGCTGTGACCATGTAACTCATCATAACCATTTACCGCGCAGCCAAGATTTCCCTGGTGATGACGCACAATTCGCGCACGCCGGTCTAGcc                                                                                             >  1:831897/1‑99 (MQ=255)
       aaaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCgg                                                                                           <  2:1568690/100‑1 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                          >  1:1936085/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                          >  1:1353601/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                          >  1:408771/1‑100 (MQ=255)
        aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                          >  1:426085/1‑100 (MQ=255)
           aCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                       <  1:1159743/100‑2 (MQ=255)
            cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgc                                                                                      <  2:1864915/100‑3 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                     <  2:444466/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                     <  1:2532105/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                     <  2:1591413/100‑4 (MQ=255)
             aGTGCTGTGCCCTTGTAACTCAGCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGGTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                     <  1:2464555/100‑4 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtt                                                                                    >  2:852515/1‑100 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtt                                                                                    <  1:156224/100‑1 (MQ=255)
              gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcgc                                                                                    <  2:2078070/100‑5 (MQ=255)
               tGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAttt                                                                                   <  1:934766/100‑1 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                  >  2:1249635/1‑100 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                  >  1:1279553/1‑100 (MQ=255)
                gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                  <  2:299032/100‑1 (MQ=255)
                 cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcgcgaa                                                                                 <  1:1636760/100‑8 (MQ=38)
                 cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAttttt                                                                                 <  1:2062265/100‑1 (MQ=255)
                    tgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAg                                                                              <  1:1064918/100‑1 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTg                                                                            <  2:2350329/100‑1 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTg                                                                            >  2:2555464/1‑100 (MQ=255)
                      cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTg                                                                            >  2:2555823/1‑100 (MQ=255)
                        cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGc                                                                          <  2:1443073/100‑1 (MQ=255)
                        cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGc                                                                          <  2:1409213/100‑1 (MQ=255)
                          tGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCtt                                                                        >  2:2211467/1‑100 (MQ=255)
                             aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                     <  1:2037564/100‑1 (MQ=255)
                             aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                     >  2:1605247/1‑100 (MQ=255)
                              aCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGc                                                                    <  1:2319515/100‑1 (MQ=255)
                               ctaatcatAATAATTGACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCATTATTCGCGCATCCCGGTCTAGCCGGGGTCATTTTTTAGt                                                                             <  1:899670/87‑1 (MQ=255)
                               cTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTCTTAGTGGCTTTTGcc                                                                   >  2:2549995/1‑100 (MQ=255)
                                   tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                               >  1:2012967/1‑100 (MQ=255)
                                   tcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                               >  1:1303153/1‑100 (MQ=255)
                                    catAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccca                                                              >  2:1692560/1‑100 (MQ=255)
                                     ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             <  1:139885/100‑1 (MQ=255)
                                     ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             >  1:240426/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  1:189616/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            <  2:575862/100‑1 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  2:240091/1‑100 (MQ=255)
                                      taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                            >  2:1473462/1‑100 (MQ=255)
                                       aataatTTACGGCGCAGCAAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                             <  1:2295812/98‑1 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       >  1:1820424/1‑100 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       <  1:95065/100‑1 (MQ=255)
                                           atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                       <  1:2552649/100‑1 (MQ=255)
                                            tttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGTGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                      <  1:174449/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1462633/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:2138472/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1844508/100‑1 (MQ=255)
                                             ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                     <  2:1844507/100‑1 (MQ=255)
                                              tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCCCCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                    <  2:1674940/100‑1 (MQ=255)
                                               aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                   <  2:605386/100‑1 (MQ=255)
                                               aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                   >  1:87304/1‑100 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTCGCCGGGGTCATCTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 <  1:1808454/100‑1 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 >  2:46548/1‑100 (MQ=255)
                                                 ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                 <  1:424559/100‑1 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  2:1311246/1‑100 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  1:1129631/1‑100 (MQ=255)
                                                  gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  1:2216510/1‑100 (MQ=255)
                                                   cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                >  2:1528489/1‑99 (MQ=255)
                                                       gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg                                           >  2:446642/1‑100 (MQ=255)
                                                         cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                         >  2:2251178/1‑100 (MQ=255)
                                                         cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                         >  1:1982812/1‑100 (MQ=255)
                                                          aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                        <  2:1226621/100‑1 (MQ=255)
                                                          aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                        <  2:2077678/100‑1 (MQ=255)
                                                              tttCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGAc                                    <  2:60601/100‑1 (MQ=255)
                                                                     ggTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCActc                             >  1:837064/1‑100 (MQ=255)
                                                                     ggTGTGGGCGCAGTATTCGCGCACCCCGGTGTGTCCGGGGTCATTTTTTAGTGGCTTGTGCCGCCCACGCTTTCAGCACTTCTACGTCGTGACGCCAc                               <  1:1985441/98‑1 (MQ=255)
                                                                          tGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTcttctt                        >  2:1774638/1‑100 (MQ=255)
                                                                           ggCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCttt                       <  1:2591319/100‑1 (MQ=255)
                                                                            gcgcAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTc                      >  2:777044/1‑100 (MQ=255)
                                                                             cgcAGTATTCGCGCATCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCa                     <  1:1009869/100‑1 (MQ=255)
                                                                             cgcAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCa                     <  1:1080903/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  2:527219/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  1:1914033/100‑1 (MQ=255)
                                                                               cAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                   <  2:335749/100‑1 (MQ=255)
                                                                                 ctaaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTATTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                 >  1:1759189/11‑100 (MQ=37)
                                                                                 ctaaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                 >  2:2153342/11‑100 (MQ=37)
                                                                                 ataaaaaatgaCCCCGGCTAGACCGGGGCCATTTGTTAGTGGCTTTTGCCACCCACGCTTTCAGCAATTATACGTCGTAACGCCACTCTTCTTTCAtctc                 >  1:246713/11‑100 (MQ=25)
                                                                                  taaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATctct                >  1:1349059/10‑100 (MQ=38)
                                                                                   aaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt               >  2:1837895/9‑100 (MQ=38)
                                                                                   aTTCGCGCACCACGGGCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACCCTTCTTTCAtct                  >  1:560909/1‑97 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:822654/100‑1 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:1544794/100‑1 (MQ=255)
                                                                                    ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc              <  2:333054/100‑1 (MQ=255)
                                                                                      cgcgCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCaa            >  2:1610667/1‑100 (MQ=255)
                                                                                       gcgcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           <  2:2194044/100‑1 (MQ=255)
                                                                                       aatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           >  2:1392731/5‑100 (MQ=255)
                                                                                       aatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt           >  1:2093232/5‑100 (MQ=255)
                                                                                        cgcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATc          <  1:1473453/100‑1 (MQ=255)
                                                                                        atgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACCCTTTCAGCACTTCTCCGTCGTTACGCCACTCTTCTTTCATCTCTTCAGTc          >  1:211742/4‑100 (MQ=255)
                                                                                         tgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATcc         >  2:1076110/3‑100 (MQ=255)
                                                                                         gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGGCGCCACTCTTCTTTCATCTCTTCAATcc         >  1:1930542/1‑100 (MQ=255)
                                                                                          cACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCa        >  1:2476790/1‑100 (MQ=255)
                                                                                                ggTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa    <  1:1116645/98‑1 (MQ=255)
                                                                                                ggTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc  <  2:400422/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACC  >  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: