Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,350,682:1 +A 16.1% intergenic (‑223/+7) 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,6821.A16.1% 371.7 / 82.2 124intergenic (‑223/+7)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 . (46/58);  new base A (10/10);  total (56/68)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 128 total aligned reads displayed.

CCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGT  >  NC_000913/3350583‑3350780
                                                                                                    |                                                                                                  
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:408771/1‑100 (MQ=255)
     aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCggg                                                                                               >  1:1353601/1‑100 (MQ=255)
        aCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtg                                                                                            <  1:1159743/100‑2 (MQ=255)
         cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgc                                                                                           <  2:1864915/100‑3 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGTCGCAGTATTCGCGCACCCCGGTCTAG‑CCgggg                                                                                              >  1:2128529/1‑96 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                          <  2:1591413/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                          <  2:444466/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                          <  1:2532105/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCAGCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGGTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcg                                                                                          <  1:2464555/100‑4 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtt                                                                                         <  1:156224/100‑1 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtt                                                                                         >  2:852515/1‑100 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcgc                                                                                         <  2:2078070/100‑5 (MQ=255)
            tGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttt                                                                                        <  1:934766/100‑1 (MQ=255)
             gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                       >  1:1279553/1‑100 (MQ=255)
             gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAtttt                                                                                       >  2:1249635/1‑100 (MQ=255)
              cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcgcgaa                                                                                      <  1:1636760/100‑8 (MQ=38)
              cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCAttttt                                                                                      <  1:2062265/100‑1 (MQ=255)
               tgtgCCCTTGTAACTCATCATACTAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGG‑CTAGACCGGGGtgcgcgaat                                                                                     <  1:1751650/100‑9 (MQ=37)
                 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)
                         tAACTCATCATAATAGTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGAGCACCCCGATCTAG‑ACGGGGTCATTTTTTAGTgg                                                                                >  1:891024/1‑95 (MQ=255)
                          aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTg                                                                          >  2:1605247/1‑100 (MQ=255)
                           aCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGc                                                                         <  1:2319515/100‑1 (MQ=255)
                            ctcctcATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGcaaaagcc                                                                        <  1:1941937/100‑8 (MQ=255)
                            ctaatcatAATAATTGACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCATTATTCGCGCATCCCGGTCTAG‑CCGGGGTCATTTTTTAGt                                                                                  <  1:899670/87‑1 (MQ=255)
                               atcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacc                                                                     <  2:1051265/100‑1 (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                                                                 >  2:1473462/1‑100 (MQ=255)
                                   taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                 >  2:240091/1‑100 (MQ=255)
                                   taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                 <  2:575862/100‑1 (MQ=255)
                                    aataatTTACGGCGCAGCAAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                                  <  1:2295812/98‑1 (MQ=255)
                                        atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                            <  1:2552649/100‑1 (MQ=255)
                                        atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                            >  1:1820424/1‑100 (MQ=255)
                                         tttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGTGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                           <  1:174449/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:1844508/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:1844507/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:2138472/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCCCCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  2:1674940/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  2:2496497/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  1:341408/100‑1 (MQ=255)
                                            aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                        <  2:605386/100‑1 (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                                                     >  1:1129631/1‑100 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                     >  2:1528489/1‑99 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTc                                                    <  1:1824516/100‑1 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTc                                                    <  1:2469474/100‑1 (MQ=255)
                                                   aGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                 >  2:1344804/1‑100 (MQ=255)
                                                    gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg                                                >  2:446642/1‑100 (MQ=255)
                                                      cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                              >  2:2251178/1‑100 (MQ=255)
                                                       aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                             <  2:2077678/100‑1 (MQ=255)
                                                       aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                             <  2:1226621/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)
                                                                         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                        <  2:335749/100‑1 (MQ=255)
                                                                            cAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                        <  1:1914033/100‑1 (MQ=255)
                                                                              ctaaaaaatgaCCCCGG‑CTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                      >  2:2153342/11‑100 (MQ=37)
                                                                               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:1544794/100‑1 (MQ=255)
                                                                                 ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                   <  2:333054/100‑1 (MQ=255)
                                                                                 ttCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                   <  2:822654/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‑CCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACGCTTTCCACCCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATc               <  1:514080/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)
                                                                                         ccccGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCAtc           <  1:2014607/100‑2 (MQ=255)
                                                                                           ccGGTCTAG‑CCAGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa         >  2:298885/1‑100 (MQ=255)
                                                                                             ggTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa         <  1:1116645/98‑1 (MQ=255)
                                                                                             ggTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc       <  2:400422/100‑1 (MQ=255)
                                                                                                cTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGa     >  1:2377711/1‑99 (MQ=255)
                                                                                                cTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGAc    <  1:2290651/100‑1 (MQ=255)
                                                                                                  aG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGt  >  2:300374/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAG‑CCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGT  >  NC_000913/3350583‑3350780

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: ‑