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.

TGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGT  >  NC_000913/4017123‑4017321
                                                                                                    |                                                                                                   
tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGt                                                                                                     <  2:3970575/100‑1 (MQ=255)
  agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                                   <  2:1902269/100‑1 (MQ=255)
  agaCGATGAAACAAACCGAAAGCCATGCGGGGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                                   <  1:1276840/100‑1 (MQ=255)
   gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTc                                                                                                   >  1:4250806/1‑99 (MQ=255)
   gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                                  >  1:1436612/1‑100 (MQ=255)
    aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCgg                                                                                                 <  2:1583461/100‑1 (MQ=255)
    aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCgg                                                                                                 <  1:704753/100‑1 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                                >  2:2762314/1‑100 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                                >  2:3132832/1‑100 (MQ=255)
     cGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGCCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGc                                                                                                >  1:2679209/1‑100 (MQ=255)
      gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGcc                                                                                               <  2:374483/100‑1 (MQ=255)
       aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCt                                                                                              >  1:4044812/1‑100 (MQ=255)
       aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCt                                                                                              <  1:2181044/100‑1 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                           <  2:3277903/100‑3 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                           <  2:4044841/100‑3 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                           >  2:90020/1‑100 (MQ=255)
          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGCCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                            <  2:3361501/99‑2 (MQ=255)
           aacacacCGAAAGCCAGGCGGTGAAAATCGTGGTGGAGGCCGCGCCTCGCGTGCTGTAATTTTCGCCTCCTTCCTGAGGGCCG‑ACGCGTTCGGCCtt                                                                                             <  1:290102/97‑1 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                          <  1:1120130/100‑4 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                          <  2:54523/100‑4 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                          >  1:660581/1‑100 (MQ=255)
           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                          >  2:2946362/1‑100 (MQ=255)
            acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                           <  1:1002661/98‑3 (MQ=255)
            acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGt                                                                                         <  1:2592427/100‑1 (MQ=255)
             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGt                                                                                         >  2:1020055/1‑99 (MQ=255)
             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCA‑ACGCGTTCGGCCTTTTGTa                                                                                        <  2:151550/100‑1 (MQ=255)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacagg                                                                                     <  2:1433323/100‑9 (MQ=38)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacagg                                                                                     <  1:2186986/100‑9 (MQ=38)
                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttt                                                                                     <  2:987399/100‑1 (MQ=255)
                  cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                                   <  2:801979/100‑1 (MQ=255)
                   gAAAGCCATGCGGTGGAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                                  <  2:4337536/100‑1 (MQ=255)
                   gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                                  >  1:3364677/1‑100 (MQ=255)
                    aaaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTAGGCCTTTTGTATTTTTGc                                                                                 <  2:1817559/100‑1 (MQ=255)
                     aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCg                                                                                <  2:2351399/100‑1 (MQ=255)
                       gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTa                                                                              >  2:296077/1‑100 (MQ=255)
                          aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                           <  1:1803983/100‑1 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  1:3749921/1‑100 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  2:1023387/1‑100 (MQ=255)
                              ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                       >  2:1213055/1‑100 (MQ=255)
                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTTATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                   <  2:3545150/100‑1 (MQ=255)
                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                   <  2:1476084/100‑1 (MQ=255)
                                   aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                                  >  2:2793878/1‑100 (MQ=255)
                                    aaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGa                                                                 <  2:613218/100‑1 (MQ=255)
                                     aTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                >  1:391753/1‑100 (MQ=255)
                                      tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                               <  1:1207080/100‑1 (MQ=255)
                                      tCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTATCTTCTCCTGTCTTAAGGCCT‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                               <  2:410881/100‑1 (MQ=255)
                                        gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTCTTTTTGCGTAGCGCCTCGCAGGAAATg                                                             >  1:1605809/1‑100 (MQ=255)
                                          ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCGCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                           <  1:214241/100‑1 (MQ=255)
                                            tggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCtt                                                         <  2:2147257/100‑1 (MQ=255)
                                              gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                       <  2:1303760/100‑1 (MQ=255)
                                              gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                       <  1:255189/100‑1 (MQ=255)
                                                aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCa                                                     <  1:3313740/100‑1 (MQ=255)
                                                  cggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                                   <  2:3327887/100‑1 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGCAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  1:1393613/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                                   <  1:46224/99‑1 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  2:2646632/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  >  1:106515/1‑100 (MQ=255)
                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                                  <  1:131093/100‑1 (MQ=255)
                                                      gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGa                                               <  2:3432080/100‑1 (MQ=255)
                                                       cgtcgtCTGCTGTAATTCTCTTCTCCTGTCGGAAGGCCG‑ACGCGTTCGGCCTTTTGGATTTTTGCGTAGCGCCTCGCAGGAAAGGCATTTcctgcgaggc                                              <  2:1690550/100‑9 (MQ=255)
                                                        gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACg                                             <  1:933292/100‑1 (MQ=255)
                                                            tcTGCTGTAATTCTCTTCTCCTGTCTGAACGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                         <  1:2641784/100‑1 (MQ=255)
                                                             ctggtgTAATTCTCTTCTCCTGTCTGCAGGCCG‑GCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                         <  1:2045050/99‑1 (MQ=255)
                                                             ctgctgTAATTCTCTTCTCCTGTCTGAAGGCTG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                        <  1:2971813/100‑1 (MQ=255)
                                                             cagctgTAATTCTATTCTCCTGTCTGAAGGCAG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                        >  2:2971829/3‑100 (MQ=255)
                                                              tgctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTa                                       <  1:2241933/100‑1 (MQ=255)
                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                      <  2:487255/100‑1 (MQ=255)
                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                      >  1:696362/1‑100 (MQ=255)
                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                     <  2:512844/100‑1 (MQ=255)
                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                     <  1:3453417/100‑1 (MQ=255)
                                                                  gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                   <  1:1592031/100‑1 (MQ=255)
                                                                  gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                   >  1:3602742/1‑100 (MQ=255)
                                                                   tAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                                  <  2:917167/100‑1 (MQ=255)
                                                                        ctcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAtt                             >  2:2464929/1‑100 (MQ=255)
                                                                            tctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAt                         >  2:2442782/1‑100 (MQ=255)
                                                                             ctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtt                        >  2:4072381/1‑100 (MQ=255)
                                                                              aaaatacaaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                       >  2:3838345/10‑100 (MQ=37)
                                                                              aaaatacaaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                       >  2:3209750/10‑100 (MQ=37)
                                                                                ctctttGAGGGCCC‑ACCCGTTCGGCCTTTTGTCATTTTTCGTAGCGCCTCTCAGGAAATGCCTTTTCAACTGAACGTTTGTCCAGCACAATTTTATCTtg                     <  1:4154084/95‑1 (MQ=255)
                                                                                  gTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgtg                   >  1:1463076/1‑100 (MQ=255)
                                                                                   tCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc                  <  2:226812/100‑1 (MQ=255)
                                                                                   acaaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc                  >  2:3141132/5‑100 (MQ=255)
                                                                                    cTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                 >  1:2435911/1‑100 (MQ=255)
                                                                                     tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTTCTACTGGTCGTTTGTACCGCACAATTCTATTTTGTGCgg                >  1:217717/1‑100 (MQ=255)
                                                                                     tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                >  2:1936913/1‑100 (MQ=255)
                                                                                     tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                 >  1:768899/1‑99 (MQ=255)
                                                                                      aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg               >  2:1092828/2‑100 (MQ=255)
                                                                                      aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg               >  2:3480733/2‑100 (MQ=255)
                                                                                      aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                >  1:3152789/2‑99 (MQ=255)
                                                                                       aaGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt              >  2:4224751/1‑100 (MQ=255)
                                                                                        aGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa             <  2:2282108/100‑1 (MQ=255)
                                                                                         ggCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa            <  1:3397552/100‑1 (MQ=255)
                                                                                          gCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg           <  2:2842408/100‑1 (MQ=255)
                                                                                            cG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgtt         >  1:3246577/1‑100 (MQ=255)
                                                                                             g‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg        >  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)
                                                                                                    |                                                                                                   
TGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGT  >  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: ‑