Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 16.2% intergenic (+17/+35) yfeD → / ← gltX DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2220TG16.2% 384.4 / 70.0 155intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (66/64);  new base G (14/11);  total (80/75)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

Only 100 of 171 total aligned reads displayed.

GCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGG  >  NC_000913/2519124‑2519321
                                                                                                  |                                                                                                   
gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                    >  2:3925699/1‑100 (MQ=255)
   ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                 <  2:3134237/100‑1 (MQ=255)
    gAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                >  2:1090642/1‑100 (MQ=255)
      aaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                              >  1:4018583/1‑100 (MQ=255)
          gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                          <  2:3493449/100‑1 (MQ=255)
           ttgcacAATATTTGTCACGTGAAGGAATTCCAAATTCCATCGTTCTTCTTAGTATTGAAGAATCAGCATAAATAAAATACGGCAGTATATTATCTCCTGc                                                                                         <  1:1899849/95‑1 (MQ=255)
           tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                         <  1:4003223/100‑1 (MQ=255)
            cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                        <  1:1080341/100‑1 (MQ=255)
            cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                        >  2:1609973/1‑100 (MQ=255)
             gcgcAATTTTTTTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                       <  2:853734/100‑1 (MQ=255)
             gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                         <  1:1160705/98‑1 (MQ=255)
             gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                       <  1:931053/100‑1 (MQ=255)
             gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                       >  2:77269/1‑100 (MQ=255)
             gcgaAATTTTTGTCACGTGAAGGCATTCGCAATTTCCTCGCCCGTCTCGATATTGACGACTCAGCATAAATAAAAAACGGCAGCATATTATCTCCTGc                                                                                         >  1:1558993/1‑98 (MQ=255)
              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                      >  2:4148782/1‑100 (MQ=255)
              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                      >  2:1677252/1‑100 (MQ=255)
              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                      <  2:2580798/100‑1 (MQ=39)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                     <  2:179549/100‑1 (MQ=39)
                 aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                   <  2:3634936/100‑1 (MQ=255)
                 aaTGTTTGTCACGGGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                   <  1:1799354/100‑1 (MQ=255)
                   tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttt                                                                                 <  1:902247/100‑4 (MQ=38)
                   tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                     <  2:4231829/96‑1 (MQ=255)
                    ttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCt                                                                                <  2:952686/100‑1 (MQ=255)
                      ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                              <  2:238565/100‑7 (MQ=37)
                       tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCTTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                             <  2:1844186/100‑1 (MQ=255)
                       tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttt                                                                              <  2:173391/99‑1 (MQ=255)
                       aggcACGAGAAGGAATTCACAATTTGATTGCCTGTCTCGATATTGACAAATAAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttt                                                                              <  1:355424/96‑1 (MQ=255)
                        gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                            >  2:3358681/1‑100 (MQ=255)
                        gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                            >  2:3868794/1‑100 (MQ=255)
                         tCACGTGAAGGAATTCGCAATTTCCTCTACCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                           <  1:1368873/100‑1 (MQ=255)
                         tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGCCAGGATATTATCTCCTGCAGTTTATCTTTTTa                                                                           <  1:3681557/100‑1 (MQ=255)
                          cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                           >  1:1659124/1‑99 (MQ=255)
                           acatGAAGGAATTCGCAATTTCCTCGCACGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                         <  1:2733936/97‑1 (MQ=255)
                            cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                        <  1:3367045/100‑1 (MQ=255)
                             gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                       >  2:376542/1‑100 (MQ=255)
                              tGAAGGACTTCGCAATTTCATCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                      <  1:3213216/100‑1 (MQ=255)
                                aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAACAAACGGCAGTATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                    >  1:3052045/1‑100 (MQ=255)
                                 aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                   <  2:2362516/100‑1 (MQ=255)
                                  ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                  <  2:2739476/100‑1 (MQ=255)
                                  ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                  >  2:3035281/1‑100 (MQ=255)
                                    aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACAGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                   <  1:1770177/97‑1 (MQ=255)
                                     aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAtt                                                                 <  2:56555/98‑1 (MQ=255)
                                       tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                             >  1:3191573/1‑100 (MQ=255)
                                        cGCAATTTCCTCGCTCGTCTCGATATTGACGAAGCAGCATGAATAAAAAACGGCAGGATATAATCTCTTGCCGTGTATCTTTTTACACGCTAATTActgc                                                            >  1:2403645/1‑100 (MQ=255)
                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                            >  2:2258607/1‑100 (MQ=255)
                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                            >  2:1661230/1‑100 (MQ=255)
                                         gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                           >  2:4266945/1‑100 (MQ=255)
                                         gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                           <  1:3687054/100‑1 (MQ=255)
                                          cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                          >  2:1347625/1‑100 (MQ=255)
                                          cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGGTTATCTTTTTACACGCTAATTActgctg                                                          <  1:2064425/100‑1 (MQ=255)
                                                ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                    >  2:1929905/1‑100 (MQ=255)
                                                 cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                   <  2:1312879/100‑1 (MQ=255)
                                                  tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                  <  2:2806991/100‑1 (MQ=255)
                                                   cGCCCGTCTCGCTATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                 >  1:2473595/1‑100 (MQ=255)
                                                   cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGTCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                 <  2:1757626/100‑1 (MQ=255)
                                                     cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                               >  2:3320408/1‑100 (MQ=255)
                                                       cGTCTCGATATTGACGAATCAGCATAAATAAAAAACTGCAAGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGCTCa                                             <  1:201131/100‑1 (MQ=255)
                                                        gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                            <  1:3901144/100‑1 (MQ=255)
                                                         tctgGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAGCTCCGGCCGTTTATCTTCTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                           <  1:4324061/100‑1 (MQ=255)
                                                         tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                           <  1:2944158/100‑1 (MQ=255)
                                                         tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                           <  1:3165496/100‑1 (MQ=255)
                                                          ctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                          >  1:2791137/1‑100 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                            <  1:4066802/97‑1 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                         >  1:127246/1‑100 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                         >  2:539153/1‑100 (MQ=255)
                                                            cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                        <  2:2813763/100‑1 (MQ=255)
                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                 <  2:1953572/100‑1 (MQ=255)
                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                 >  2:4026473/1‑100 (MQ=255)
                                                                     cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                               >  1:2835014/1‑100 (MQ=255)
                                                                     cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                               >  2:2805556/1‑100 (MQ=255)
                                                                      gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                              <  2:3182397/100‑1 (MQ=255)
                                                                            gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGTGTTCAGCAATAAAATCCAGCGCTttg                        <  1:940231/100‑1 (MQ=255)
                                                                             gtaaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                       >  2:3926900/11‑100 (MQ=25)
                                                                             cATAAATAAAAAACGGCAGGATATTAGCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                       <  1:1507919/100‑1 (MQ=255)
                                                                               taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                     >  2:1391073/1‑100 (MQ=255)
                                                                               taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                     <  2:2266687/100‑1 (MQ=255)
                                                                               aaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                     >  2:472517/9‑100 (MQ=37)
                                                                               aaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                     >  2:2465247/9‑100 (MQ=25)
                                                                                 aataaaAACCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTTTTGAt                   >  1:1067530/1‑100 (MQ=255)
                                                                                  aagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttgct                   >  1:1214746/6‑97 (MQ=37)
                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTTTTGAt                   >  1:4293023/1‑98 (MQ=255)
                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTAAGCAATAAAATCCAGCGCTTTGTTGATAc                 >  1:3904089/1‑100 (MQ=255)
                                                                                    aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGCTACg                >  1:2619171/1‑100 (MQ=255)
                                                                                    aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                >  1:2230415/1‑100 (MQ=255)
                                                                                      taaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt              >  2:3684883/2‑100 (MQ=38)
                                                                                       aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             <  1:394201/100‑1 (MQ=255)
                                                                                        aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             >  1:514490/1‑99 (MQ=38)
                                                                                        aaCGGCAGGAGATAATATCCTGCCGTTTATATTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  1:2848523/1‑100 (MQ=38)
                                                                                         aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTTTTTTTTATTTATTGCTGATTTTCGCGTTCAGGAATATAATGCAGTGGTTTGTTGATACGCTCg            >  1:3454349/1‑99 (MQ=16)
                                                                                          cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          >  1:2432311/1‑100 (MQ=255)
                                                                                          cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          >  2:4300815/1‑100 (MQ=255)
                                                                                          cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          <  1:64169/100‑1 (MQ=255)
                                                                                          cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          >  1:1613618/1‑100 (MQ=255)
                                                                                            gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa         >  1:3245688/1‑99 (MQ=255)
                                                                                            gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc        >  1:3363496/1‑100 (MQ=255)
                                                                                              aGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa           <  2:1527883/95‑1 (MQ=255)
                                                                                                  tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTTATACGCTCGATACTGCg    >  1:3847738/1‑98 (MQ=255)
                                                                                                  tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  >  2:4138147/1‑100 (MQ=255)
                                                                                                  tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  >  2:3354079/1‑100 (MQ=255)
                                                                                                  tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  <  2:210555/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGG  >  NC_000913/2519124‑2519321

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