Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,217:1 +A 15.0% intergenic (+25/‑116) 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,2161.A15.0% 668.2 / 133.0 213intergenic (+24/‑117)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base . (97/84);  new base A (14/18);  total (111/102)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-02

Only 100 of 216 total aligned reads displayed.

GAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGT  >  NC_000913/4017117‑4017313
                                                                                                    |                                                                                                 
   tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acg                                                                                                <  2:3184148/100‑1 (MQ=255)
    gCTGAGACGATGAAACAAACCGAAAGCCATTCGGTGAAAATCGTGGTGGAAGCGACGCGACGTCTGCAGTAATTCTCTTCTCCTGTCTGCAGGCCG‑Acgc                                                                                               >  2:3248109/1‑100 (MQ=255)
    gCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acgc                                                                                               >  2:1019699/1‑100 (MQ=255)
     cTGCGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acgc                                                                                               >  2:4023884/1‑99 (MQ=255)
     cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acgcg                                                                                              >  2:3063292/1‑100 (MQ=255)
     cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑Acgcg                                                                                              >  2:2195193/1‑100 (MQ=255)
      tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGt                                                                                             >  2:1998467/1‑100 (MQ=255)
        agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGtc                                                                                            <  2:1933059/100‑2 (MQ=255)
         gaCGATGAAACAAACTGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                          <  2:4345116/100‑1 (MQ=255)
         gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                          >  1:3764699/1‑100 (MQ=255)
            gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGcc                                                                                       <  1:2116426/100‑1 (MQ=255)
              tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                     <  1:1568156/100‑1 (MQ=255)
              tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                     >  2:2484390/1‑100 (MQ=255)
                agacaaacCGATAGCCATGCGGTGAAAATCGTGGTGGAAGCTGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                   <  1:3432156/98‑3 (MQ=255)
                aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                   <  1:554343/100‑3 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                  <  2:2428019/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                  <  2:859651/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                  <  1:2478415/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                  <  2:3767456/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcag                                                                                  <  1:4274486/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                  <  1:3863094/100‑1 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                  >  1:13285/1‑100 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTg                                                                                  >  2:1396777/1‑100 (MQ=255)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagaca                                                                               <  2:2719191/100‑7 (MQ=38)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagaca                                                                               <  1:2257086/100‑7 (MQ=38)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTAGTCTGCTGTAATTCTCTTCTCCTGTCTGACGGCCG‑ACGCGTTCGGCCTTTTGtct                                                                               >  1:2031559/1‑98 (MQ=255)
                     aacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacag                                                                              <  2:2255373/100‑8 (MQ=38)
                      acCGAAAGCCCTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGGCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttcagacagg                                                                             <  1:3940853/100‑9 (MQ=37)
                       ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtttt                                                                            >  2:2690340/1‑100 (MQ=255)
                        cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAttttt                                                                           >  2:3628546/1‑100 (MQ=255)
                         gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTg                                                                          <  1:850547/100‑1 (MQ=255)
                           aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCg                                                                        <  1:4165986/100‑1 (MQ=255)
                            aGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGt                                                                       <  1:1126262/100‑1 (MQ=255)
                              ccATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAg                                                                     >  1:2540293/1‑100 (MQ=255)
                              ccATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAg                                                                     >  1:821423/1‑100 (MQ=255)
                               cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGGAGGCCT‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                    >  1:2981608/1‑100 (MQ=255)
                               cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                    <  1:1619378/100‑1 (MQ=255)
                                   cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCt                                                                >  2:1331450/1‑100 (MQ=255)
                                   cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCt                                                                >  1:2784874/1‑100 (MQ=255)
                                    ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                               >  1:2605450/1‑100 (MQ=255)
                                     gTGAAAATAGTGGTGGCAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTCTGCGTAGCGCCTCg                                                              >  1:2624083/1‑100 (MQ=255)
                                       gAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCa                                                            <  1:4114472/100‑1 (MQ=255)
                                         aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                          >  2:178730/1‑100 (MQ=255)
                                            tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                       >  2:1335679/1‑100 (MQ=255)
                                              gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAAtn                                                     <  1:3147660/100‑2 (MQ=255)
                                               tggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                    <  2:2936821/100‑1 (MQ=255)
                                                ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                   >  1:506953/1‑100 (MQ=255)
                                                   ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCttt                                                <  2:1548117/100‑1 (MQ=255)
                                                    gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                               >  2:3996636/1‑100 (MQ=255)
                                                    gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                               <  2:3106680/100‑1 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  1:653422/1‑100 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  2:3358244/1‑100 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  1:2337548/1‑100 (MQ=255)
                                                          gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTg                                         <  1:3313632/100‑1 (MQ=255)
                                                              gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACg                                     >  2:1704939/1‑100 (MQ=255)
                                                                cgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGtt                                   <  1:3651997/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:2780912/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:1425225/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:3577115/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 >  1:1460292/1‑100 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              <  1:1823251/100‑1 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              >  2:4171846/1‑100 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              <  2:4300285/100‑1 (MQ=255)
                                                                      ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                             >  2:3508834/1‑100 (MQ=255)
                                                                      ctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                             <  2:2607612/100‑1 (MQ=255)
                                                                       tgTAATTCTCTTCTCCTGTCTGATGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAg                            <  2:3964274/100‑1 (MQ=255)
                                                                        gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                           <  1:358901/100‑1 (MQ=255)
                                                                        gTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                           <  1:4389190/100‑1 (MQ=255)
                                                                         tAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                          >  2:487013/1‑100 (MQ=255)
                                                                         tAATTCTCTTCTCCTGTCTGAAGGCCA‑ACCCGTTCGGCTTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                          >  1:765905/1‑100 (MQ=255)
                                                                          aaTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcac                         <  2:1584078/100‑1 (MQ=255)
                                                                          aaTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACGGGACGTTTGTACAGcac                         <  2:4138872/100‑1 (MQ=255)
                                                                           aTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTTTGTTTTGTTTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                        >  1:1379044/1‑100 (MQ=255)
                                                                               tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                    <  1:3972048/100‑1 (MQ=255)
                                                                               tcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                    <  2:1772988/100‑1 (MQ=255)
                                                                                cttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCt                   <  2:4081556/100‑1 (MQ=255)
                                                                                cttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGCTTGTACAgatcggaaga                   >  2:1569555/1‑91 (MQ=255)
                                                                                    tcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt               <  2:2778183/100‑1 (MQ=255)
                                                                                      cTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtg             <  2:254382/100‑1 (MQ=255)
                                                                                       tGTCTTAAGGCCC‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTCGACGTTTGTACAACACAATTCTATTTtgt            >  1:3549465/1‑100 (MQ=255)
                                                                                        gTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgtg           >  2:3779315/1‑100 (MQ=255)
                                                                                         tCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc          <  2:1066482/100‑1 (MQ=255)
                                                                                          caaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  2:1153549/4‑100 (MQ=255)
                                                                                          caaaaGGCCGAACGCG‑TCGGCCTTATGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGTACGTTTGTACAGCACAATTCTATTTTGTGc          >  1:1688817/4‑99 (MQ=255)
                                                                                          caaaaGGCCGAACGCG‑TCGGCATTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGAACGTTTGTACAGCACAATTCAATTTTTTGCg         >  1:1003701/4‑100 (MQ=255)
                                                                                          caaaaGGCCGAACGCG‑TCGGACTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  1:2694623/4‑100 (MQ=255)
                                                                                          cTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  1:28587/1‑100 (MQ=255)
                                                                                          cTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGTTTGTAAAGCAAAATTCTATTTtgtg           >  2:2866077/1‑98 (MQ=255)
                                                                                           tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  2:1864242/1‑100 (MQ=255)
                                                                                           tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  2:2891224/1‑100 (MQ=255)
                                                                                           tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  1:277823/1‑100 (MQ=255)
                                                                                           naaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGACATGCCTTTCCAACTGGACGTTTGTACAGCCCAATTCTATTTtgtgttg        >  1:133/3‑97 (MQ=255)
                                                                                           aaaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  1:4057718/3‑100 (MQ=255)
                                                                                             acggCCG‑ACGCGTCCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCCTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg       <  2:821424/97‑1 (MQ=255)
                                                                                             aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt      >  2:3454849/1‑100 (MQ=255)
                                                                                             aaGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt      <  2:569282/100‑1 (MQ=255)
                                                                                              aGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa     >  1:4099979/1‑100 (MQ=255)
                                                                                              aGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa     <  2:3807822/100‑1 (MQ=255)
                                                                                                 ccG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgt  <  1:3073388/100‑1 (MQ=255)
                                                                                                 ccG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg   >  1:107161/1‑99 (MQ=255)
                                                                                                    |                                                                                                 
GAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGT  >  NC_000913/4017117‑4017313

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: ‑