Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,039 (C)7→6 100% coding (585/720 nt) yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0330C.75.5% 208.7 / 96.1 94coding (591/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (8/14);  new base . (23/48);  total (31/63)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 101 total aligned reads displayed.

TATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCC  >  NC_000913/3100936‑3101128
                                                                                                 |                                                                                               
tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcc                                                                                                >  1:2554762/1‑99 (MQ=255)
  tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccc                                                                                             <  2:1344409/100‑1 (MQ=255)
   cTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccc                                                                                            >  2:1384706/1‑100 (MQ=255)
     tCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccag                                                                                          <  2:2105697/100‑3 (MQ=255)
      ccAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccagg                                                                                         >  1:1483216/1‑97 (MQ=255)
       cAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccaggc                                                                                        <  2:2498499/100‑5 (MQ=255)
       cAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccaggc                                                                                        >  2:2202600/1‑96 (MQ=255)
        agagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCt                                                                                      <  2:2078357/100‑1 (MQ=255)
        agagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCATGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCt                                                                                      >  1:1548299/1‑100 (MQ=255)
           gTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTcc                                                                                   >  1:2508370/1‑100 (MQ=255)
            tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTcc                                                                                   <  1:420294/100‑1 (MQ=255)
            tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTccc                                                                                  <  1:1484944/100‑1 (MQ=255)
            tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTccc                                                                                  <  2:1370012/100‑1 (MQ=255)
            tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTccc                                                                                  <  2:20193/100‑1 (MQ=255)
             cacaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCa                                                                                 <  2:222614/100‑1 (MQ=255)
             cacaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCa                                                                                 >  2:924919/1‑100 (MQ=255)
              acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAg                                                                                <  1:332500/100‑1 (MQ=255)
              acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAg                                                                                <  1:2166884/100‑1 (MQ=255)
              acaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAg                                                                                >  2:2442863/1‑100 (MQ=255)
               caACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTCTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGc                                                                               >  1:1682345/1‑100 (MQ=255)
                 aCGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGcac                                                                             <  2:2188987/100‑1 (MQ=255)
                  cgcgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGcac                                                                             <  2:366418/100‑1 (MQ=255)
                  cgcgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGcaca                                                                            <  2:748085/100‑1 (MQ=255)
                   gcgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACAt                                                                           >  2:1492527/1‑100 (MQ=255)
                      gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTc                                                                         <  2:1861683/100‑1 (MQ=255)
                      gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCt                                                                        <  2:2388707/100‑1 (MQ=255)
                          cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  2:480385/1‑100 (MQ=255)
                          cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  2:1946445/1‑100 (MQ=255)
                          cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  2:603782/1‑100 (MQ=255)
                          cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  1:1019169/1‑100 (MQ=255)
                            aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATg                                                                  <  1:2152818/100‑1 (MQ=255)
                              cATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                                <  1:142664/100‑1 (MQ=255)
                              cATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                                <  2:1579159/100‑1 (MQ=255)
                                 cTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTa                                                             <  2:1344335/100‑1 (MQ=255)
                                    tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccg                                                          <  2:1183084/100‑1 (MQ=255)
                                     gCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCCCTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccg                                                          <  2:1080453/99‑1 (MQ=255)
                                      ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                         <  2:705172/100‑1 (MQ=255)
                                      ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                        <  2:2575386/100‑1 (MQ=255)
                                      ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                        >  2:2075685/1‑100 (MQ=255)
                                       cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                        <  1:2253476/100‑1 (MQ=255)
                                       cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                       >  1:1209632/1‑100 (MQ=255)
                                         aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACcgccgc                                                      >  2:1983584/1‑100 (MQ=255)
                                           ctgctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtt                                                    <  2:370982/100‑1 (MQ=255)
                                             gctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                 <  2:286472/100‑1 (MQ=255)
                                             gctgGAAATCTTTTTCCTGCGTTTTCCACGCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCTCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                 <  1:1115996/100‑1 (MQ=255)
                                              ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                 >  1:779459/1‑100 (MQ=255)
                                              ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                 <  2:1579953/100‑1 (MQ=255)
                                              ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                <  1:1946435/100‑1 (MQ=255)
                                              ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                <  1:1835786/100‑1 (MQ=255)
                                              ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                <  1:188876/100‑1 (MQ=255)
                                               tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGATGATTGCAG‑CCCCCCGTGCTTTCCAGCACATTCTGTTATGGGTTACCGCCGCtttt                                                  <  1:2491520/97‑1 (MQ=255)
                                               tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                <  1:2468173/100‑1 (MQ=255)
                                               tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                               <  1:1554146/100‑1 (MQ=255)
                                                ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              >  2:2357100/1‑100 (MQ=255)
                                                ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              >  1:248548/1‑100 (MQ=255)
                                                  aaaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                            <  1:1558152/100‑1 (MQ=255)
                                                  aaaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                            <  1:966292/100‑1 (MQ=255)
                                                   aaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAcc                                           <  2:1981051/100‑1 (MQ=255)
                                                    aTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                          >  2:491581/1‑100 (MQ=255)
                                                     tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCg                                          >  2:1211139/1‑100 (MQ=255)
                                                     tCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGc                                         >  2:1737150/1‑100 (MQ=255)
                                                     tCTTTTGCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTCCAGCACCGc                                         <  1:2200848/100‑1 (MQ=255)
                                                           tCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgc                                   <  2:2401100/100‑1 (MQ=255)
                                                            ccTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                  <  2:292454/100‑1 (MQ=255)
                                                            ccTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcgcg                                  <  2:1715783/100‑1 (MQ=255)
                                                              tGCTTTTTCTACTCGGTTTGGATTGAGGATTGCAGCCCGCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCCCTTTcgcg                                  >  1:430752/1‑99 (MQ=255)
                                                              tGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                <  2:820071/100‑1 (MQ=255)
                                                              tGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                <  2:249430/100‑1 (MQ=255)
                                                              tGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                <  2:1161944/100‑1 (MQ=255)
                                                              tGCTTTTTCCACTCGGTTTGGATTGAGGATTGC‑GGCCCCCCGGGCTTCCCACCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCgcgc                                 <  1:2291436/99‑1 (MQ=255)
                                                               gCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                >  1:1644937/1‑100 (MQ=255)
                                                                cTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCa                              <  1:1769700/100‑1 (MQ=255)
                                                                 tttttCCACTCGGTTCGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                             <  2:2554781/100‑1 (MQ=255)
                                                                   tttCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                           >  1:316528/1‑100 (MQ=255)
                                                                    ttCCACTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTc                          <  2:248550/100‑1 (MQ=255)
                                                                       cACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                        <  1:1915491/100‑1 (MQ=255)
                                                                       cACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                        <  2:1463375/100‑1 (MQ=255)
                                                                        aCTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAtt                       <  1:1150288/100‑1 (MQ=255)
                                                                         cTCGGTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                     >  2:649629/1‑100 (MQ=255)
                                                                          tCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                     <  1:383168/100‑1 (MQ=255)
                                                                           cGGTTTGGATTGAGGGTTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                   <  2:2557208/100‑1 (MQ=255)
                                                                            ggTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGc                  <  2:1645814/100‑1 (MQ=255)
                                                                            ggTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGc                  <  1:2357085/100‑1 (MQ=255)
                                                                             gTTTGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                 >  1:1062223/1‑100 (MQ=255)
                                                                               ttGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                <  1:1169974/100‑1 (MQ=255)
                                                                               ttGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGt               <  1:232003/100‑1 (MQ=255)
                                                                                tGGATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTc              <  2:5682/100‑1 (MQ=255)
                                                                                 ggATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTc              >  2:1497308/1‑100 (MQ=255)
                                                                                 ggATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTcc             <  2:2500926/100‑1 (MQ=255)
                                                                                 ggATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTcc             <  2:316530/100‑1 (MQ=255)
                                                                                 ggATTGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTcc             <  1:2505213/100‑1 (MQ=255)
                                                                                  gATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTcc             <  2:1927589/100‑1 (MQ=255)
                                                                                   aTTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTCTCCt            >  1:2527339/1‑100 (MQ=255)
                                                                                     tGAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa         >  1:863850/1‑100 (MQ=255)
                                                                                      gAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa         >  2:1320688/1‑100 (MQ=255)
                                                                                      gAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa        >  1:621276/1‑100 (MQ=255)
                                                                                      gAGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa        >  2:856099/1‑100 (MQ=255)
                                                                                       aGGATTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGATTTCGCGCCCATTTCATTAATGCTGTCCTGTaaa       <  1:1272062/100‑1 (MQ=255)
                                                                                          aTTGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAAtt    <  1:1313797/100‑1 (MQ=255)
                                                                                            tGCAG‑CCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTcc  <  1:2075671/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCC  >  NC_000913/3100936‑3101128

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