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.

CCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGT  >  NC_000913/3350583‑3350780
                                                                                                    |                                                                                                  
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:408771/1‑100 (MQ=255)
     aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                               >  1:1353601/1‑100 (MQ=255)
        aCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                            <  1:1159743/100‑2 (MQ=255)
         cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgc                                                                                           <  2:1864915/100‑3 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGTCGCAGTATTCGCGCACCCCGGTCTAGCCgggg                                                                                              >  1:2128529/1‑96 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                          <  2:1591413/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                          <  2:444466/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                          <  1:2532105/100‑4 (MQ=255)
          aGTGCTGTGCCCTTGTAACTCAGCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGGTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcg                                                                                          <  1:2464555/100‑4 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtt                                                                                         <  1:156224/100‑1 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtt                                                                                         >  2:852515/1‑100 (MQ=255)
           gTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcgc                                                                                         <  2:2078070/100‑5 (MQ=255)
            tGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAttt                                                                                        <  1:934766/100‑1 (MQ=255)
             gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                       >  1:1279553/1‑100 (MQ=255)
             gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                       >  2:1249635/1‑100 (MQ=255)
              cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCCGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcgcgaa                                                                                      <  1:1636760/100‑8 (MQ=38)
              cTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAttttt                                                                                      <  1:2062265/100‑1 (MQ=255)
               tgtgCCCTTGTAACTCATCATACTAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtgcgcgaat                                                                                     <  1:1751650/100‑9 (MQ=37)
                 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)
                         tAACTCATCATAATAGTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGAGCACCCCGATCTAGACGGGGTCATTTTTTAGTgg                                                                                >  1:891024/1‑95 (MQ=255)
                          aaCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                          >  2:1605247/1‑100 (MQ=255)
                           aCTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGc                                                                         <  1:2319515/100‑1 (MQ=255)
                            ctcctcATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGcaaaagcc                                                                        <  1:1941937/100‑8 (MQ=255)
                            ctaatcatAATAATTGACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCATTATTCGCGCATCCCGGTCTAGCCGGGGTCATTTTTTAGt                                                                                  <  1:899670/87‑1 (MQ=255)
                               atcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacc                                                                     <  2:1051265/100‑1 (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                                                                 >  2:1473462/1‑100 (MQ=255)
                                   taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                 >  2:240091/1‑100 (MQ=255)
                                   taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                 <  2:575862/100‑1 (MQ=255)
                                    aataatTTACGGCGCAGCAAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                                                  <  1:2295812/98‑1 (MQ=255)
                                        atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                            <  1:2552649/100‑1 (MQ=255)
                                        atTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTc                                                            >  1:1820424/1‑100 (MQ=255)
                                         tttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGTGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                           <  1:174449/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:1844508/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:1844507/100‑1 (MQ=255)
                                          ttACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                          <  2:2138472/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCCCCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  2:1674940/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  2:2496497/100‑1 (MQ=255)
                                           tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                                                         <  1:341408/100‑1 (MQ=255)
                                            aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCa                                                        <  2:605386/100‑1 (MQ=255)
                                              ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTCGCCGGGGTCATCTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                      <  1:1808454/100‑1 (MQ=255)
                                              ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                      >  2:46548/1‑100 (MQ=255)
                                              ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                      <  1:424559/100‑1 (MQ=255)
                                               gcgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                     >  1:1129631/1‑100 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                     >  2:1528489/1‑99 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTc                                                    <  1:1824516/100‑1 (MQ=255)
                                                cgcAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTc                                                    <  1:2469474/100‑1 (MQ=255)
                                                   aGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                 >  2:1344804/1‑100 (MQ=255)
                                                    gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg                                                >  2:446642/1‑100 (MQ=255)
                                                      cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc                                              >  2:2251178/1‑100 (MQ=255)
                                                       aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                             <  2:2077678/100‑1 (MQ=255)
                                                       aaGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                             <  2:1226621/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)
                                                                         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                        <  2:335749/100‑1 (MQ=255)
                                                                            cAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtc                        <  1:1914033/100‑1 (MQ=255)
                                                                              ctaaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAtctc                      >  2:2153342/11‑100 (MQ=37)
                                                                               taaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATctct                     >  1:1349059/10‑100 (MQ=38)
                                                                                aaaaaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                    >  2:1837895/9‑100 (MQ=38)
                                                                                aTTCGCGCACCACGGGCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACCCTTCTTTCAtct                       >  1:560909/1‑97 (MQ=255)
                                                                                 ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                   <  2:1544794/100‑1 (MQ=255)
                                                                                 ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                   <  2:333054/100‑1 (MQ=255)
                                                                                 ttCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                   <  2:822654/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)
                                                                                     cgcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACGCTTTCCACCCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATc               <  1:514080/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)
                                                                                         ccccGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCAtc           <  1:2014607/100‑2 (MQ=255)
                                                                                           ccGGTCTAGCCAGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa         >  2:298885/1‑100 (MQ=255)
                                                                                             ggTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCa         <  1:1116645/98‑1 (MQ=255)
                                                                                             ggTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc       <  2:400422/100‑1 (MQ=255)
                                                                                                cTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGa     >  1:2377711/1‑99 (MQ=255)
                                                                                                cTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGAc    <  1:2290651/100‑1 (MQ=255)
                                                                                                  aGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGt  >  2:300374/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
CCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGT  >  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: