Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,217:1 +A 15.6% 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.6% 674.6 / 141.1 212intergenic (+24/‑117)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base . (71/108);  new base A (16/17);  total (87/125)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 216 total aligned reads displayed.

AATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTG  >  NC_000913/4017118‑4017315
                                                                                                   |                                                                                                   
    cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAcgcg                                                                                                >  1:4267070/1‑100 (MQ=255)
     tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGt                                                                                               <  2:3970575/100‑1 (MQ=255)
     tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGt                                                                                               >  1:2842392/1‑100 (MQ=255)
       agaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                             <  2:1902269/100‑1 (MQ=255)
       agaCGATGAAACAAACCGAAAGCCATGCGGGGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                             <  1:1276840/100‑1 (MQ=255)
        gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                            <  1:2145209/100‑1 (MQ=255)
        gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                            <  2:504127/100‑1 (MQ=255)
         aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCgg                                                                                           <  2:1583461/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                                                                                         <  1:2059557/100‑1 (MQ=255)
            aTGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGGCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCt                                                                                        <  1:2181044/100‑1 (MQ=255)
              gAAACAAACCGAAAGCCCTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTTAAGGTCTTCGCTTTTGGCCttt                                                                                      <  1:2562008/100‑1 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                     <  2:4044841/100‑3 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                     <  2:3277903/100‑3 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCtttt                                                                                     >  2:90020/1‑100 (MQ=255)
               aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGCCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttc                                                                                      <  2:3361501/99‑2 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                    <  1:1120130/100‑4 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                    <  2:54523/100‑4 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                    <  2:467687/100‑1 (MQ=255)
                aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTg                                                                                    >  2:2946362/1‑100 (MQ=255)
                 acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                     <  1:1002661/98‑3 (MQ=255)
                 acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGt                                                                                   >  2:1497494/1‑100 (MQ=255)
                  caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTa                                                                                  >  2:4087133/1‑100 (MQ=255)
                   aaagCGAAAGGCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                     <  1:4053060/96‑3 (MQ=255)
                    aacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtt                                                                                <  2:1989390/100‑1 (MQ=255)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacagg                                                                               <  1:2186986/100‑9 (MQ=38)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagacagg                                                                               <  2:1433323/100‑9 (MQ=38)
                     acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttt                                                                               <  2:987399/100‑1 (MQ=255)
                      ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttt                                                                              <  1:526254/100‑1 (MQ=255)
                       cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttttt                                                                             <  2:801979/100‑1 (MQ=255)
                        gacaggcATGCGGTGAAAATCGTGGTGGAAGCTGCGCGTCGTCTGCGGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                    <  1:2765090/86‑4 (MQ=37)
                        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)
                          aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGAGTCGTCTGCTCTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttttt                                                                            <  1:679246/98‑2 (MQ=255)
                            gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTa                                                                        >  2:296077/1‑100 (MQ=255)
                               aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                     <  1:1803983/100‑1 (MQ=255)
                                   ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                 >  2:1213055/1‑100 (MQ=255)
                                    gTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCg                                                                <  2:3943569/100‑1 (MQ=255)
                                     tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                               <  1:2496065/100‑1 (MQ=255)
                                     tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                               >  2:1082944/1‑100 (MQ=255)
                                       aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTTATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                             <  2:3545150/100‑1 (MQ=255)
                                        aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                            >  2:2246596/1‑100 (MQ=255)
                                        aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                             <  2:1010254/99‑1 (MQ=255)
                                             gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTCTTTTTGCGTAGCGCCTCGCAGGAAATg                                                       >  1:1605809/1‑100 (MQ=255)
                                               ggtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                     <  1:3132815/100‑1 (MQ=255)
                                                  ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCttt                                                  >  2:1376648/1‑100 (MQ=255)
                                                   gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                 <  2:1303760/100‑1 (MQ=255)
                                                     aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCa                                               <  1:3313740/100‑1 (MQ=255)
                                                       cggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                             <  2:3327887/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAAc                                             <  1:46224/99‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            <  2:3832228/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            >  2:2646632/1‑100 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            <  1:131093/100‑1 (MQ=255)
                                                        ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACt                                            >  1:106515/1‑100 (MQ=255)
                                                         gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGGAGGAAATGCCTTTCCAACTg                                           <  2:4140972/100‑1 (MQ=255)
                                                         gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTg                                           <  2:2454145/100‑1 (MQ=255)
                                                           gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGa                                         <  1:3247104/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)
                                                                  ctgctgTAATTCTCTTCTCCTGTCTGCAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  >  1:3606387/1‑100 (MQ=255)
                                                                  ctgctgTAATTCTCTTCTCCTGTCTGAAGGCTGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  <  1:2971813/100‑1 (MQ=255)
                                                                  cagctgTAATTCTATTCTCCTGTCTGAAGGCAGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                  >  2:2971829/3‑100 (MQ=255)
                                                                    gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                <  2:3199198/100‑1 (MQ=255)
                                                                    gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                >  1:699425/1‑100 (MQ=255)
                                                                     ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                               <  1:3453417/100‑1 (MQ=255)
                                                                      tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                               >  2:3429626/1‑99 (MQ=255)
                                                                      tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAg                              <  1:1458241/100‑1 (MQ=255)
                                                                       gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                             <  1:1592031/100‑1 (MQ=255)
                                                                       gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                             <  2:3472214/100‑1 (MQ=255)
                                                                        tAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                            <  2:917167/100‑1 (MQ=255)
                                                                          aTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                          >  2:28403/1‑100 (MQ=255)
                                                                             ctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAtt                       >  2:2464929/1‑100 (MQ=255)
                                                                               cttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCt                     >  1:2400837/1‑100 (MQ=255)
                                                                                   tcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                 <  1:148555/100‑1 (MQ=255)
                                                                                   aaaatacaaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAttt                 >  2:3209750/10‑100 (MQ=37)
                                                                                    ccTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                >  1:1157482/1‑100 (MQ=255)
                                                                                     cTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                <  2:3738318/99‑1 (MQ=255)
                                                                                       gTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTGGTACAGCACAAATCTATTTtgtg             >  1:422899/1‑100 (MQ=255)
                                                                                         caaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg           >  1:3501540/4‑100 (MQ=255)
                                                                                         caaaaGGCCGAACGCGTCGGCCTATTGTATTTTTGAGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGAGTGTACAACACAATTCTATTTtgtgag           >  1:3211728/4‑98 (MQ=255)
                                                                                          tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTTCTACTGGTCGTTTGTACCGCACAATTCTATTTTGTGCgg          >  1:217717/1‑100 (MQ=255)
                                                                                          tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg          >  2:2924966/1‑100 (MQ=255)
                                                                                          tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg          >  2:1936913/1‑100 (MQ=255)
                                                                                          tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg           >  1:768899/1‑99 (MQ=255)
                                                                                           aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTTCggg         >  1:4018294/2‑100 (MQ=255)
                                                                                             aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa       <  2:2282108/100‑1 (MQ=255)
                                                                                              ggCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa      >  2:659148/1‑100 (MQ=255)
                                                                                              ggCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa      <  1:3397003/100‑1 (MQ=255)
                                                                                              ggCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa       >  1:1428646/1‑99 (MQ=255)
                                                                                               gCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg     <  1:1213050/100‑1 (MQ=255)
                                                                                               gCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg     <  2:1020549/100‑1 (MQ=255)
                                                                                                ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgt    >  1:3859531/1‑100 (MQ=255)
                                                                                                  gACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg  >  1:902429/1‑100 (MQ=255)
                                                                                                  gACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAACACAATTCTATTTTGTGCGGGTAAgttg  <  2:4238284/100‑1 (MQ=255)
                                                                                                  gACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAAGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg  >  2:154992/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTG  >  NC_000913/4017118‑4017315

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: