Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,244 C→T 15.0% intergenic (+50/+179) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2440CT15.0% 391.7 / 42.9 133intergenic (+50/+179)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (56/57);  new base T (18/2);  total (74/59)
Fisher's exact test for biased strand distribution p-value = 5.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 138 total aligned reads displayed.

ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NC_000913/3423145‑3423343
                                                                                                   |                                                                                                   
aTCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTc                                                                                                     <  2:1052560/100‑1 (MQ=255)
   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGTATGGCCTTTCGCt                                                                                                  >  2:1508058/1‑100 (MQ=255)
   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCt                                                                                                  >  2:2863912/1‑100 (MQ=255)
   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCt                                                                                                  >  2:1641549/1‑100 (MQ=255)
   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCt                                                                                                  >  1:1922960/1‑100 (MQ=255)
    tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCtt                                                                                                 >  1:1223992/1‑100 (MQ=255)
        tCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatt                                                                                             <  1:1413225/100‑1 (MQ=255)
          agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttg                                                                                           <  1:2215055/100‑1 (MQ=255)
           gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttga                                                                                          >  1:2831730/1‑100 (MQ=255)
            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                         <  1:101029/100‑1 (MQ=255)
            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                         <  2:2927225/100‑1 (MQ=255)
            agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                         <  1:3071098/100‑1 (MQ=255)
             gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                         <  2:951003/99‑1 (MQ=255)
             gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATg                                                                                        <  1:295527/100‑1 (MQ=255)
               gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTc                                                                                      >  2:871479/1‑100 (MQ=255)
                 aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTg                                                                                    <  2:727022/100‑1 (MQ=255)
                 aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTg                                                                                    <  2:714041/100‑1 (MQ=255)
                 aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTg                                                                                    >  2:2806868/1‑100 (MQ=255)
                 aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTg                                                                                    <  2:1078587/100‑1 (MQ=255)
                  cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTgg                                                                                   >  1:624921/1‑100 (MQ=255)
                    aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCa                                                                                 >  2:2300369/1‑100 (MQ=255)
                    aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCa                                                                                 >  1:2289471/1‑100 (MQ=255)
                      ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGt                                                                               <  1:2070649/100‑1 (MQ=255)
                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGttt                                                                             <  2:1134907/100‑1 (MQ=255)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTa                                                                            <  2:2496176/100‑1 (MQ=255)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTa                                                                            >  2:2890929/1‑100 (MQ=255)
                          tttttATCGCAGGTAATTCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAt                                                                           >  1:1481115/1‑100 (MQ=255)
                          tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAt                                                                           <  1:2267652/100‑1 (MQ=255)
                          tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAt                                                                           >  1:1661453/1‑100 (MQ=255)
                           ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATg                                                                          <  2:681655/100‑1 (MQ=255)
                            tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATgg                                                                         >  1:307654/1‑100 (MQ=255)
                             ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggc                                                                        <  1:1973710/100‑1 (MQ=255)
                               aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                      >  2:385123/1‑100 (MQ=255)
                               aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                      >  1:2587001/1‑100 (MQ=255)
                                 cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                    <  1:709752/100‑1 (MQ=255)
                                  gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggcn                                                                   <  1:2767407/100‑2 (MQ=255)
                                  gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggcg                                                                   >  1:1726522/1‑100 (MQ=255)
                                       tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTg                                                              >  1:1198037/1‑100 (MQ=255)
                                          tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGccc                                                           >  1:2234918/1‑100 (MQ=255)
                                           tcATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCg                                                          <  1:2463336/99‑1 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCTTCCTGccc                                                           >  2:2608246/1‑98 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCga                                                         >  2:2463312/1‑99 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGc                                                         >  1:3034804/1‑100 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGc                                                         >  2:1229590/1‑100 (MQ=255)
                                                 aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                    <  2:403753/100‑1 (MQ=255)
                                                 aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                    >  1:165720/1‑100 (MQ=255)
                                                 aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                    <  2:1458686/100‑1 (MQ=255)
                                                   tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTc                                                  >  1:1946416/1‑100 (MQ=255)
                                                    gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                 >  1:1307323/1‑100 (MQ=255)
                                                     aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                >  1:2839186/1‑100 (MQ=255)
                                                       tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                              <  1:320620/100‑1 (MQ=255)
                                                        gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                             <  1:2635883/100‑1 (MQ=255)
                                                         ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGTGCGTCCTGCCCGCCACCCTCCGGGcc                                            >  2:1909418/1‑100 (MQ=255)
                                                           aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                            >  2:2004923/1‑98 (MQ=255)
                                                           aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                           >  1:1616569/1‑99 (MQ=255)
                                                           aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                           >  1:2872286/1‑99 (MQ=255)
                                                            gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCggt                                         >  2:1891646/1‑98 (MQ=255)
                                                            gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                         >  2:387862/1‑100 (MQ=255)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                      >  1:262881/1‑100 (MQ=255)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                      <  2:1756371/100‑1 (MQ=255)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATAGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                      <  1:28507/100‑1 (MQ=255)
                                                                  aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                   <  1:903433/100‑1 (MQ=255)
                                                                    aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                 <  2:221752/100‑1 (MQ=255)
                                                                     gCAAAAAGGCCATCCTGTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTGTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                <  2:1424574/100‑1 (MQ=255)
                                                                      gcagaaagCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAAc                               >  2:53856/8‑100 (MQ=1)
                                                                      gcagaaGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                               >  2:250632/5‑100 (MQ=1)
                                                                       cagaaGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                              >  1:1768874/4‑100 (MQ=1)
                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCCCAAc                               >  2:2757645/1‑99 (MQ=255)
                                                                        aaaaGGCCATCCTTTCGGATGGCCTGTCGCTTGAGTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCAGGGCCGTTGCTTCGCaan                               <  2:1770092/98‑2 (MQ=255)
                                                                           aagCCATCCTGACAGATGGCCTTCTGCTTAAGTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  1:2094171/3‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  1:1221824/1‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  1:1774123/1‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  1:436367/1‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  1:319620/1‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  2:2563907/1‑100 (MQ=1)
                                                                           aGGCCATCCTGACGGATGGCCTTCTGCTTAATTCGATGCCTGGTAGTTTATTGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGTTc                           >  2:2988545/1‑99 (MQ=1)
                                                                            ggCCATCCTGACGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                         >  1:763752/1‑100 (MQ=1)
                                                                             gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATCTCTGGCAGTTTATGGCTGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                        <  2:1543064/100‑1 (MQ=255)
                                                                              ccATCCTTTCGGATGTCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCTGGCGTCCTGCCCGCCAACCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                       >  2:744182/1‑100 (MQ=255)
                                                                              ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAAt                       >  2:2571957/1‑100 (MQ=255)
                                                                              ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                          >  2:160543/1‑97 (MQ=255)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                      >  1:620449/1‑100 (MQ=255)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                      >  1:514041/1‑100 (MQ=255)
                                                                                aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                     <  1:507799/100‑1 (MQ=255)
                                                                                   cgaaaggATGGCCTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                  >  1:3074142/6‑100 (MQ=255)
                                                                                   cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                  >  2:1989599/1‑100 (MQ=255)
                                                                                   cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                  >  2:2565946/1‑100 (MQ=255)
                                                                                     ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                <  1:1352538/100‑1 (MQ=255)
                                                                                     ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                <  2:537572/100‑1 (MQ=255)
                                                                                      tCGGATGGCCTTTCGCTTGAATTGATGTCTGGCAGTGTATGGCGGGCGTCCTGCCCGCCACCCTCCTTGCTGTTGCTTCGCAACGTTCAAATCCGc                   >  2:574038/1‑96 (MQ=14)
                                                                                       cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg              <  1:1904263/100‑1 (MQ=255)
                                                                                        gggtGGGCTGTCGCTTGATGTGATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACGTTCAAATCCGCTCCcgg             <  2:2831701/97‑1 (MQ=11)
                                                                                           tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg          <  1:2541981/100‑1 (MQ=255)
                                                                                            ggCCTTTCGCTTGATTTGATTTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa         <  1:106845/100‑1 (MQ=255)
                                                                                            ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg           <  2:2845054/98‑1 (MQ=255)
                                                                                            ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc            >  2:655169/1‑97 (MQ=255)
                                                                                              ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  1:1008598/100‑1 (MQ=255)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  >  2:1213201/1‑100 (MQ=40)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:2890958/100‑1 (MQ=40)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:2423480/100‑1 (MQ=40)
                                                                                                   |                                                                                                   
ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NC_000913/3423145‑3423343

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: