Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 15.5% intergenic (+20/+32) 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,2250TA15.5% 385.3 / 59.3 156intergenic (+20/+32)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (67/64);  new base A (9/15);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 170 total aligned reads displayed.

AGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  NC_000913/2519126‑2519324
                                                                                                   |                                                                                                   
aGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatt                                                                                                     >  1:4192958/1‑100 (MQ=255)
 ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                    >  2:786877/1‑100 (MQ=255)
   aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                  <  2:978500/100‑1 (MQ=255)
        gTCGCCCAATTTTTGTCACGTGAAGTGATTCTCAATTTCCTCGCCCTGCTGGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATAGTAtctc                                                                                                >  1:2518706/1‑97 (MQ=255)
        gTCGAGCAATTTTTGTAACGTAAATTAATTCGCAATCTCACTGCTCGTCTCGATATTGTCAAATCACAATAAATAAAAAACGGCAGGATTTTAtctct                                                                                               <  1:2083922/98‑2 (MQ=255)
         tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                            <  2:3369383/100‑1 (MQ=255)
           gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                          <  2:404791/100‑1 (MQ=255)
            cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                         <  1:3224828/100‑1 (MQ=39)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                        <  2:2502543/100‑1 (MQ=255)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  2:3495985/100‑1 (MQ=255)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATGTCCGCGCCCGGCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                         <  1:277151/99‑1 (MQ=38)
              cAATTTTTGTCTCGTGAAGGAATTGGCAATGTCCGCCCACGTCTCGATATTGACGAATCAGCTTCAATAAAAAACGGCAGGAGATTATCTCGTGCCGttt                                                                                       >  1:658668/1‑100 (MQ=255)
              cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                       >  1:1025263/1‑100 (MQ=255)
               aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                      <  2:3660932/100‑2 (MQ=38)
               aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                      <  1:2408325/100‑2 (MQ=38)
                aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                     <  1:4402830/100‑3 (MQ=38)
                aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                     <  2:2698706/100‑3 (MQ=38)
                 tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttt                                                                                    <  2:2039020/100‑4 (MQ=38)
                    ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                                 <  2:3008346/100‑7 (MQ=37)
                     tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                                >  2:3227483/1‑100 (MQ=255)
                      gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                               <  2:712795/100‑1 (MQ=255)
                        cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                             >  1:4400640/1‑100 (MQ=255)
                         aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                            <  2:536326/100‑1 (MQ=255)
                             gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                        >  2:3780476/1‑100 (MQ=255)
                              aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                       <  1:4384639/100‑1 (MQ=255)
                                ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                     <  2:174868/100‑1 (MQ=255)
                                   aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                  <  2:420513/100‑1 (MQ=255)
                                   aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                  >  1:2061912/1‑100 (MQ=255)
                                    ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                 >  1:2733538/1‑100 (MQ=255)
                                    ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                 >  2:543964/1‑100 (MQ=255)
                                    ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                 >  1:1247129/1‑100 (MQ=255)
                                     tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                <  1:1473395/100‑1 (MQ=255)
                                      cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               <  1:971631/100‑1 (MQ=255)
                                      cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               >  2:443640/1‑100 (MQ=255)
                                       gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                              >  2:3254938/1‑100 (MQ=255)
                                       gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                              <  2:2297298/100‑1 (MQ=255)
                                        cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                             >  2:2578828/1‑100 (MQ=255)
                                          aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                           <  2:2834296/100‑1 (MQ=255)
                                          aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                           <  1:2723987/100‑1 (MQ=255)
                                           tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtt                                                          >  1:1211605/1‑100 (MQ=255)
                                            ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                         <  2:798143/100‑1 (MQ=255)
                                               cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                      <  2:2572483/100‑1 (MQ=255)
                                                 cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGGTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                    <  1:2347204/100‑1 (MQ=255)
                                                   cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                  >  1:2207842/1‑100 (MQ=255)
                                                    ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                 <  2:3890371/100‑1 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  2:3201083/100‑1 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  2:2733553/100‑1 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  1:14458/100‑1 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  2:435818/100‑1 (MQ=255)
                                                        ctcGCTATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                             >  1:701082/1‑100 (MQ=255)
                                                          cgctatTGCCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                           >  1:1134822/4‑100 (MQ=255)
                                                            atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                         >  1:924074/1‑100 (MQ=255)
                                                             tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        >  1:1629767/1‑100 (MQ=255)
                                                             tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        >  2:2476547/1‑100 (MQ=255)
                                                              atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                       >  1:1916405/1‑100 (MQ=255)
                                                               ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                      <  2:4302787/100‑1 (MQ=255)
                                                                 gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                    <  1:783415/100‑1 (MQ=255)
                                                                 gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                    >  2:1822616/1‑100 (MQ=255)
                                                                    gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                 >  2:1265238/1‑100 (MQ=255)
                                                                    gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                 <  2:4137163/100‑1 (MQ=255)
                                                                      aTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCTCGTTCAGCAATAAAATCCAgcgc                               >  1:357958/1‑100 (MQ=255)
                                                                       tCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                              >  2:3701990/1‑100 (MQ=255)
                                                                        cAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                             >  2:2071170/1‑100 (MQ=255)
                                                                         aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTTCTGATTTTCGCGATCAGCAATAAAATCCAGCGCtt                             >  1:4147927/1‑99 (MQ=255)
                                                                          gtataaataaAAAACGGCAGCATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTGTCGCGTTCAGCAAGAAAATCCAGCGCtt                             <  1:3381460/96‑1 (MQ=255)
                                                                          gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                           >  2:3612405/1‑100 (MQ=255)
                                                                          gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                           <  1:2980253/100‑1 (MQ=255)
                                                                          gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                           >  1:2206970/1‑100 (MQ=255)
                                                                           cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                          <  1:3156245/100‑1 (MQ=255)
                                                                            ataaataaAAACCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  2:678078/100‑1 (MQ=255)
                                                                            ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         >  2:3293800/1‑100 (MQ=255)
                                                                            ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  1:2899597/100‑1 (MQ=255)
                                                                             naaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         >  1:2480643/2‑99 (MQ=255)
                                                                              aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATTCAGCGCTttgttg                        >  1:3636091/1‑99 (MQ=255)
                                                                              aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  2:4058245/1‑100 (MQ=255)
                                                                              aaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  2:1288135/8‑100 (MQ=37)
                                                                                 taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAtgc                    >  1:980868/1‑98 (MQ=255)
                                                                                 taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    >  2:2173582/1‑100 (MQ=255)
                                                                                 taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAACCAGCGCTTTGTTGATAc                    >  1:4271226/1‑100 (MQ=255)
                                                                                 agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    >  2:1223209/5‑100 (MQ=37)
                                                                                  gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGGTTTTTGCGTTTGGAACTAAAATCCAGCGCTTTGTTGATa                     >  1:3793120/4‑98 (MQ=21)
                                                                                   aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                  >  1:3896678/1‑100 (MQ=255)
                                                                                    aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                 <  2:1688865/100‑1 (MQ=255)
                                                                                      aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                <  1:2883481/99‑1 (MQ=255)
                                                                                       ncGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGCTTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACTCTCGa              >  1:2071043/2‑100 (MQ=38)
                                                                                        cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             <  1:3006673/100‑1 (MQ=255)
                                                                                          gCAGGATATTATCTCCTGCCGTTTATCTTTTTTCACGCTGATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc           <  1:3085504/100‑1 (MQ=255)
                                                                                          gCAGGATAGTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCTTTCAGCAATAAAATCCAGTGCTTTGTTGATACTTTCTATAc           >  1:1454622/1‑100 (MQ=255)
                                                                                             ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTAGAGAc           >  1:2051479/1‑97 (MQ=255)
                                                                                              gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       >  1:3374104/1‑100 (MQ=255)
                                                                                              gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       <  1:4023012/100‑1 (MQ=255)
                                                                                              gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       <  1:1648718/100‑1 (MQ=255)
                                                                                                tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg     <  2:2382228/100‑1 (MQ=255)
                                                                                                tattatCTCCTGCCGTTTATCTTTTTACACGCCAATTACGGCCGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg     <  1:3085964/100‑1 (MQ=255)
                                                                                                  ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc   >  1:1736930/1‑100 (MQ=255)
                                                                                                  ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc   <  1:1193044/100‑1 (MQ=255)
                                                                                                  ctatCTCCTGCCGTTTATCTTTCTACACGCTAATTACTGCTTATTTTCGCTCTCAGCTATACAATCCACTGCTTTGTTTATACGCTCGATACTGCGGGTc   <  1:4318750/99‑1 (MQ=255)
                                                                                                   tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  <  1:2611621/100‑1 (MQ=255)
                                                                                                   tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  >  2:1496518/1‑100 (MQ=255)
                                                                                                   tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  >  1:721226/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  NC_000913/2519126‑2519324

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: