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.

GAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGT  >  NC_000913/4017117‑4017313
                                                                                                    |                                                                                                 
   tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcg                                                                                                <  2:3184148/100‑1 (MQ=255)
    gCTGAGACGATGAAACAAACCGAAAGCCATTCGGTGAAAATCGTGGTGGAAGCGACGCGACGTCTGCAGTAATTCTCTTCTCCTGTCTGCAGGCCGAcgc                                                                                               >  2:3248109/1‑100 (MQ=255)
    gCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcgc                                                                                               >  2:1019699/1‑100 (MQ=255)
     cTGCGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcgc                                                                                               >  2:4023884/1‑99 (MQ=255)
     cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcgcg                                                                                              >  2:3063292/1‑100 (MQ=255)
     cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcgcg                                                                                              >  2:2195193/1‑100 (MQ=255)
      tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGt                                                                                             >  2:1998467/1‑100 (MQ=255)
        agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGtc                                                                                            <  2:1933059/100‑2 (MQ=255)
         gaCGATGAAACAAACTGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                          <  2:4345116/100‑1 (MQ=255)
         gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                          >  1:3764699/1‑100 (MQ=255)
            gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGcc                                                                                       <  1:2116426/100‑1 (MQ=255)
              tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCtt                                                                                     <  1:1568156/100‑1 (MQ=255)
              tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCtt                                                                                     >  2:2484390/1‑100 (MQ=255)
                agacaaacCGATAGCCATGCGGTGAAAATCGTGGTGGAAGCTGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                   <  1:3432156/98‑3 (MQ=255)
                aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                   <  1:554343/100‑3 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                  <  2:2428019/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                  <  2:859651/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                  <  1:2478415/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                  <  2:3767456/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                  <  1:4274486/100‑4 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                  <  1:3863094/100‑1 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                  >  1:13285/1‑100 (MQ=255)
                 aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                  >  2:1396777/1‑100 (MQ=255)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagaca                                                                               <  2:2719191/100‑7 (MQ=38)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagaca                                                                               <  1:2257086/100‑7 (MQ=38)
                    aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTAGTCTGCTGTAATTCTCTTCTCCTGTCTGACGGCCGACGCGTTCGGCCTTTTGtct                                                                               >  1:2031559/1‑98 (MQ=255)
                     aacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacag                                                                              <  2:2255373/100‑8 (MQ=38)
                      acCGAAAGCCCTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGGCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacagg                                                                             <  1:3940853/100‑9 (MQ=37)
                       ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttt                                                                            >  2:2690340/1‑100 (MQ=255)
                        cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttttt                                                                           >  2:3628546/1‑100 (MQ=255)
                         gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                          <  1:850547/100‑1 (MQ=255)
                           aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCg                                                                        <  1:4165986/100‑1 (MQ=255)
                            aGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGt                                                                       <  1:1126262/100‑1 (MQ=255)
                              ccATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAg                                                                     >  1:2540293/1‑100 (MQ=255)
                              ccATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAg                                                                     >  1:821423/1‑100 (MQ=255)
                               cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGGAGGCCTACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                    >  1:2981608/1‑100 (MQ=255)
                               cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                    <  1:1619378/100‑1 (MQ=255)
                                   cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCt                                                                >  2:1331450/1‑100 (MQ=255)
                                   cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCt                                                                >  1:2784874/1‑100 (MQ=255)
                                    ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                               >  1:2605450/1‑100 (MQ=255)
                                     gTGAAAATAGTGGTGGCAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTCTGCGTAGCGCCTCg                                                              >  1:2624083/1‑100 (MQ=255)
                                       gAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCa                                                            <  1:4114472/100‑1 (MQ=255)
                                         aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                          >  2:178730/1‑100 (MQ=255)
                                            tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                       >  2:1335679/1‑100 (MQ=255)
                                              gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAAtn                                                     <  1:3147660/100‑2 (MQ=255)
                                               tggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                    <  2:2936821/100‑1 (MQ=255)
                                                ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                   >  1:506953/1‑100 (MQ=255)
                                                   ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCttt                                                <  2:1548117/100‑1 (MQ=255)
                                                    gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                               >  2:3996636/1‑100 (MQ=255)
                                                    gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                               <  2:3106680/100‑1 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  1:653422/1‑100 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  2:3358244/1‑100 (MQ=255)
                                                         ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                          >  1:2337548/1‑100 (MQ=255)
                                                          gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTg                                         <  1:3313632/100‑1 (MQ=255)
                                                              gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACg                                     >  2:1704939/1‑100 (MQ=255)
                                                                cgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGtt                                   <  1:3651997/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:2780912/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:1425225/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 <  2:3577115/100‑1 (MQ=255)
                                                                  tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                 >  1:1460292/1‑100 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              <  1:1823251/100‑1 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              >  2:4171846/1‑100 (MQ=255)
                                                                     gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                              <  2:4300285/100‑1 (MQ=255)
                                                                      ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                             >  2:3508834/1‑100 (MQ=255)
                                                                      ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                             <  2:2607612/100‑1 (MQ=255)
                                                                       tgTAATTCTCTTCTCCTGTCTGATGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAg                            <  2:3964274/100‑1 (MQ=255)
                                                                        gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                           <  1:358901/100‑1 (MQ=255)
                                                                        gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                           <  1:4389190/100‑1 (MQ=255)
                                                                         tAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                          >  2:487013/1‑100 (MQ=255)
                                                                         tAATTCTCTTCTCCTGTCTGAAGGCCAACCCGTTCGGCTTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                          >  1:765905/1‑100 (MQ=255)
                                                                          aaTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcac                         <  2:1584078/100‑1 (MQ=255)
                                                                          aaTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACGGGACGTTTGTACAGcac                         <  2:4138872/100‑1 (MQ=255)
                                                                           aTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTTTGTTTTGTTTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                        >  1:1379044/1‑100 (MQ=255)
                                                                               tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                    <  1:3972048/100‑1 (MQ=255)
                                                                               tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                    <  2:1772988/100‑1 (MQ=255)
                                                                                cttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCt                   <  2:4081556/100‑1 (MQ=255)
                                                                                cttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGCTTGTACAgatcggaaga                   >  2:1569555/1‑91 (MQ=255)
                                                                                    tcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt               <  2:2778183/100‑1 (MQ=255)
                                                                                      cTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtg             <  2:254382/100‑1 (MQ=255)
                                                                                       tGTCTTAAGGCCCACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTCGACGTTTGTACAACACAATTCTATTTtgt            >  1:3549465/1‑100 (MQ=255)
                                                                                        gTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgtg           >  2:3779315/1‑100 (MQ=255)
                                                                                         tCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGc          <  2:1066482/100‑1 (MQ=255)
                                                                                          caaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  2:1153549/4‑100 (MQ=255)
                                                                                          caaaaGGCCGAACGCGTCGGCCTTATGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGTACGTTTGTACAGCACAATTCTATTTTGTGc          >  1:1688817/4‑99 (MQ=255)
                                                                                          caaaaGGCCGAACGCGTCGGCATTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGAACGTTTGTACAGCACAATTCAATTTTTTGCg         >  1:1003701/4‑100 (MQ=255)
                                                                                          caaaaGGCCGAACGCGTCGGACTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  1:2694623/4‑100 (MQ=255)
                                                                                          cTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg         >  1:28587/1‑100 (MQ=255)
                                                                                          cTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGTTTGTAAAGCAAAATTCTATTTtgtg           >  2:2866077/1‑98 (MQ=255)
                                                                                           tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  2:1864242/1‑100 (MQ=255)
                                                                                           tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  2:2891224/1‑100 (MQ=255)
                                                                                           tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  1:277823/1‑100 (MQ=255)
                                                                                           naaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGACATGCCTTTCCAACTGGACGTTTGTACAGCCCAATTCTATTTtgtgttg        >  1:133/3‑97 (MQ=255)
                                                                                           aaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg        >  1:4057718/3‑100 (MQ=255)
                                                                                             acggCCGACGCGTCCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCCTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg       <  2:821424/97‑1 (MQ=255)
                                                                                             aaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt      >  2:3454849/1‑100 (MQ=255)
                                                                                             aaGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt      <  2:569282/100‑1 (MQ=255)
                                                                                              aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa     >  1:4099979/1‑100 (MQ=255)
                                                                                              aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa     <  2:3807822/100‑1 (MQ=255)
                                                                                                 ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgt  <  1:3073388/100‑1 (MQ=255)
                                                                                                 ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg   >  1:107161/1‑99 (MQ=255)
                                                                                                    |                                                                                                 
GAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGT  >  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: