Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,223 Δ1 bp 15.2% intergenic (+31/‑110) udp → / → rmuC uridine phosphorylase/DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,2220T.15.2% 628.9 / 111.3 217intergenic (+30/‑111)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (72/111);  new base . (16/17);  total (88/129)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 221 total aligned reads displayed.

TGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGT  >  NC_000913/4017123‑4017321
                                                                                                    |                                                                                                   
tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGt                                                                                                     <  2:3970575/100‑1 (MQ=255)
  agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                   <  2:1902269/100‑1 (MQ=255)
  agaCGATGAAACAAACCGAAAGCCATGCGGGGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                   <  1:1276840/100‑1 (MQ=255)
   gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                   >  1:4250806/1‑99 (MQ=255)
   gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                                  >  1:1436612/1‑100 (MQ=255)
    aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCgg                                                                                                 <  2:1583461/100‑1 (MQ=255)
    aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCgg                                                                                                 <  1:704753/100‑1 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGc                                                                                                >  2:2762314/1‑100 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGc                                                                                                >  2:3132832/1‑100 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGCCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGc                                                                                                >  1:2679209/1‑100 (MQ=255)
      gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGcc                                                                                               <  2:374483/100‑1 (MQ=255)
       aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCt                                                                                              >  1:4044812/1‑100 (MQ=255)
       aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCt                                                                                              <  1:2181044/100‑1 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                           <  2:3277903/100‑3 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                           <  2:4044841/100‑3 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCtttt                                                                                           >  2:90020/1‑100 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGCCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttc                                                                                            <  2:3361501/99‑2 (MQ=255)
           aacacacCGAAAGCCAGGCGGTGAAAATCGTGGTGGAGGCCGCGCCTCGCGTGCTGTAATTTTCGCCTCCTTCCTGAGGGCCGACGCGTTCGGCCtt                                                                                             <  1:290102/97‑1 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                          <  1:1120130/100‑4 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                          <  2:54523/100‑4 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                          >  1:660581/1‑100 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                          >  2:2946362/1‑100 (MQ=255)
            acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                           <  1:1002661/98‑3 (MQ=255)
            acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGt                                                                                         <  1:2592427/100‑1 (MQ=255)
             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGt                                                                                         >  2:1020055/1‑99 (MQ=255)
             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCAACGCGTTCGGCCTTTTGTa                                                                                        <  2:151550/100‑1 (MQ=255)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacagg                                                                                     <  2:1433323/100‑9 (MQ=38)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacagg                                                                                     <  1:2186986/100‑9 (MQ=38)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttt                                                                                     <  2:987399/100‑1 (MQ=255)
                  cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttttt                                                                                   <  2:801979/100‑1 (MQ=255)
                   gAAAGCCATGCGGTGGAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                                  <  2:4337536/100‑1 (MQ=255)
                   gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                                  >  1:3364677/1‑100 (MQ=255)
                    aaaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTAGGCCTTTTGTATTTTTGc                                                                                 <  2:1817559/100‑1 (MQ=255)
                     aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCg                                                                                <  2:2351399/100‑1 (MQ=255)
                       gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTa                                                                              >  2:296077/1‑100 (MQ=255)
                          aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                           <  1:1803983/100‑1 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  1:3749921/1‑100 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  2:1023387/1‑100 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  2:1213055/1‑100 (MQ=255)
                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTTATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                   <  2:3545150/100‑1 (MQ=255)
                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                   <  2:1476084/100‑1 (MQ=255)
                                   aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                                  >  2:2793878/1‑100 (MQ=255)
                                    aaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGa                                                                 <  2:613218/100‑1 (MQ=255)
                                     aTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                >  1:391753/1‑100 (MQ=255)
                                      tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                               <  1:1207080/100‑1 (MQ=255)
                                      tCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTATCTTCTCCTGTCTTAAGGCCTACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                               <  2:410881/100‑1 (MQ=255)
                                        gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTCTTTTTGCGTAGCGCCTCGCAGGAAATg                                                             >  1:1605809/1‑100 (MQ=255)
                                          ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCGCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                           <  1:214241/100‑1 (MQ=255)
                                            tggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCtt                                                         <  2:2147257/100‑1 (MQ=255)
                                              gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                       <  2:1303760/100‑1 (MQ=255)
                                              gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                       <  1:255189/100‑1 (MQ=255)
                                                aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCa                                                     <  1:3313740/100‑1 (MQ=255)
                                                  cggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                                   <  2:3327887/100‑1 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGCAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  1:1393613/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                                   <  1:46224/99‑1 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  2:2646632/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  1:106515/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  <  1:131093/100‑1 (MQ=255)
                                                      gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGa                                               <  2:3432080/100‑1 (MQ=255)
                                                       cgtcgtCTGCTGTAATTCTCTTCTCCTGTCGGAAGGCCGACGCGTTCGGCCTTTTGGATTTTTGCGTAGCGCCTCGCAGGAAAGGCATTTcctgcgaggc                                              <  2:1690550/100‑9 (MQ=255)
                                                        gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACg                                             <  1:933292/100‑1 (MQ=255)
                                                            tcTGCTGTAATTCTCTTCTCCTGTCTGAACGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                         <  1:2641784/100‑1 (MQ=255)
                                                             ctggtgTAATTCTCTTCTCCTGTCTGCAGGCCGGCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                         <  1:2045050/99‑1 (MQ=255)
                                                             ctgctgTAATTCTCTTCTCCTGTCTGAAGGCTGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                        <  1:2971813/100‑1 (MQ=255)
                                                             cagctgTAATTCTATTCTCCTGTCTGAAGGCAGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                        >  2:2971829/3‑100 (MQ=255)
                                                              tgctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTa                                       <  1:2241933/100‑1 (MQ=255)
                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                      <  2:487255/100‑1 (MQ=255)
                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                      >  1:696362/1‑100 (MQ=255)
                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                     <  2:512844/100‑1 (MQ=255)
                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                     <  1:3453417/100‑1 (MQ=255)
                                                                  gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                   <  1:1592031/100‑1 (MQ=255)
                                                                  gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                   >  1:3602742/1‑100 (MQ=255)
                                                                   tAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                                  <  2:917167/100‑1 (MQ=255)
                                                                        ctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAtt                             >  2:2464929/1‑100 (MQ=255)
                                                                            tctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAt                         >  2:2442782/1‑100 (MQ=255)
                                                                             ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtt                        >  2:4072381/1‑100 (MQ=255)
                                                                              aaaatacaaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                       >  2:3838345/10‑100 (MQ=37)
                                                                              aaaatacaaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                       >  2:3209750/10‑100 (MQ=37)
                                                                                ctctttGAGGGCCCACCCGTTCGGCCTTTTGTCATTTTTCGTAGCGCCTCTCAGGAAATGCCTTTTCAACTGAACGTTTGTCCAGCACAATTTTATCTtg                     <  1:4154084/95‑1 (MQ=255)
                                                                                  gTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgtg                   >  1:1463076/1‑100 (MQ=255)
                                                                                   tCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc                  <  2:226812/100‑1 (MQ=255)
                                                                                   acaaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc                  >  2:3141132/5‑100 (MQ=255)
                                                                                    cTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                 >  1:2435911/1‑100 (MQ=255)
                                                                                     tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTTCTACTGGTCGTTTGTACCGCACAATTCTATTTTGTGCgg                >  1:217717/1‑100 (MQ=255)
                                                                                     tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                >  2:1936913/1‑100 (MQ=255)
                                                                                     tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                 >  1:768899/1‑99 (MQ=255)
                                                                                      aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg               >  2:1092828/2‑100 (MQ=255)
                                                                                      aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg               >  2:3480733/2‑100 (MQ=255)
                                                                                      aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                >  1:3152789/2‑99 (MQ=255)
                                                                                       aaGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt              >  2:4224751/1‑100 (MQ=255)
                                                                                        aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa             <  2:2282108/100‑1 (MQ=255)
                                                                                         ggCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa            <  1:3397552/100‑1 (MQ=255)
                                                                                          gCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg           <  2:2842408/100‑1 (MQ=255)
                                                                                            cGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgtt         >  1:3246577/1‑100 (MQ=255)
                                                                                             gACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg        >  1:902429/1‑100 (MQ=255)
                                                                                                 gcgTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGttgttg     <  1:28403/100‑1 (MQ=255)
                                                                                                   gTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCg   <  2:1556829/100‑1 (MQ=255)
                                                                                                   gTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCg   >  1:1148869/1‑100 (MQ=255)
                                                                                                    ttCGGCCTTTTGTATTTTTTTGTATTTTTTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTTCGt  >  1:336962/1‑100 (MQ=255)
                                                                                                    ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGt  >  1:524506/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGT  >  NC_000913/4017123‑4017321

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