Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 15.7% 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.7% 216.5 / 31.1 89intergenic (+20/+32)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base A (5/9);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

AGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  NC_000913/2519126‑2519324
                                                                                                   |                                                                                                   
aGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatt                                                                                                     >  1:1821386/1‑100 (MQ=255)
 ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                    <  1:1827628/100‑1 (MQ=255)
   aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                  >  2:818040/1‑100 (MQ=255)
      aaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATAGTATCTcc                                                                                               >  1:1052931/1‑100 (MQ=255)
        gTCGCGCAATTTTTGTCACGTGAAGGAATTCGGAATTTCCTCGCCCGTCTCGATATTGACGAATCAGAATCAATAAAACACGGCAGGATATTAtctca                                                                                               >  2:1587361/1‑97 (MQ=255)
        gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                             >  1:1164850/1‑100 (MQ=255)
        gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                             >  1:1338102/1‑100 (MQ=255)
        gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                             >  2:2293284/1‑100 (MQ=255)
        gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAACATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                             >  1:228993/1‑100 (MQ=255)
          cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCTCGTCTTGATATCGACGAATCAGCATAAATAAACATCGGCAGGATATTATCTCCTGcc                                                                                           >  2:1982545/1‑100 (MQ=255)
           gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTGCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                          <  2:987631/100‑1 (MQ=255)
           gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAAAAAAAAAAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                          <  2:1883305/100‑1 (MQ=255)
            cgcAATTTTTGTCACTTTAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                         <  1:1491731/100‑1 (MQ=38)
            cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                         <  2:392939/100‑1 (MQ=39)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTTGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  1:1796121/100‑1 (MQ=38)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                        <  1:403428/100‑1 (MQ=255)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  2:814274/100‑1 (MQ=39)
             gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  1:815875/100‑1 (MQ=39)
               aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                      >  1:1775078/1‑100 (MQ=255)
                aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                     <  1:1446399/100‑3 (MQ=38)
                 tttgtGTCACGTGAAGGAATTCGCAATTTCCTCACCAGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTAt                                                                                     <  2:1088817/99‑1 (MQ=255)
                    ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                                 <  1:450783/100‑7 (MQ=37)
                     tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGACAGGATATTATCTCCTGCCGTTTATCtttt                                                                                <  1:942050/100‑1 (MQ=255)
                     tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGAGATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttattt                                                                                <  2:49179/100‑8 (MQ=25)
                       tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                              <  1:363556/100‑1 (MQ=255)
                       tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                              <  1:1639751/100‑1 (MQ=255)
                        cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATAATATCTCCTGCCGTTTATCTTTTTac                                                                             <  2:2249138/100‑1 (MQ=255)
                           gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGTTATTATCTCCTGCCGTTTATCTTTTTACACg                                                                          >  2:1383667/1‑100 (MQ=255)
                            tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAACAAACGCCAGGATATTATCTCCGCCAGATTATCtttacacc                                                                            >  2:1740420/1‑92 (MQ=255)
                            tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         <  1:2236482/100‑1 (MQ=255)
                            tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         >  2:2496900/1‑100 (MQ=255)
                              aaGGATTTCAAAATTTCCTCGCCCGTTTCAATATTAATTAATCAGTATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                       <  2:1539334/100‑1 (MQ=255)
                               aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                      >  2:1769384/1‑100 (MQ=255)
                                ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                     >  2:1430849/1‑100 (MQ=255)
                                  aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                   <  1:1619795/100‑1 (MQ=255)
                                  aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                   <  1:1241728/100‑1 (MQ=255)
                                    ttCGCAATTTCTTCGTCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                 <  2:569799/100‑1 (MQ=255)
                                     tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAACTAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                >  2:1267438/1‑100 (MQ=255)
                                      cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               >  1:1526363/1‑100 (MQ=255)
                                      cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               >  2:322725/1‑100 (MQ=255)
                                       gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAAGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               <  2:757309/99‑1 (MQ=255)
                                       gCAACTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGTCCGACTCTTATCTCCTGCTGGTTATCTTTTTACACGCTAATTAc                                                                  >  2:2362889/1‑96 (MQ=255)
                                             tCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATATCATGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                        >  2:1803035/1‑100 (MQ=255)
                                               cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                      <  2:86110/100‑1 (MQ=255)
                                                tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                     <  1:11993/100‑1 (MQ=255)
                                                 cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                    >  2:1923263/1‑100 (MQ=255)
                                                  gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                   >  2:1201579/1‑100 (MQ=255)
                                                    ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                 <  2:2099377/100‑1 (MQ=255)
                                                     cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                <  1:1294645/100‑1 (MQ=255)
                                                     cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                <  2:308648/100‑1 (MQ=255)
                                                     cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                <  1:1740432/100‑1 (MQ=255)
                                                      gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                               >  1:427504/1‑100 (MQ=255)
                                                       ttttgATATTGATGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGTGttta                                                <  2:1213939/94‑3 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  1:1615774/100‑1 (MQ=255)
                                                       tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  2:1052923/100‑1 (MQ=255)
                                                        ctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                             <  1:577930/100‑1 (MQ=255)
                                                        ctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                             <  1:1845900/100‑1 (MQ=255)
                                                          cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTGAGCAAt                                           >  1:1533586/1‑100 (MQ=255)
                                                           gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATa                                          <  1:748184/100‑1 (MQ=255)
                                                              atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                       <  2:506715/100‑1 (MQ=255)
                                                                tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTTCTTGTTTTTTACTGCTGATTTTTGCGTTCATCAATAAAAt                                      >  2:853284/1‑99 (MQ=255)
                                                                 gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                    <  2:1562967/100‑1 (MQ=255)
                                                                    gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                 <  1:303980/100‑1 (MQ=255)
                                                                      aTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcgc                               >  1:362819/1‑100 (MQ=255)
                                                                       tCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                              >  2:214137/1‑100 (MQ=255)
                                                                        cAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                             <  1:1878495/100‑1 (MQ=255)
                                                                          gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                            >  2:2355288/1‑99 (MQ=255)
                                                                            ataaataaAAAACGGCAGGATATTATCTCCTTCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  2:2113796/100‑1 (MQ=255)
                                                                            ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  2:1086829/100‑1 (MQ=255)
                                                                            ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  1:2282064/100‑1 (MQ=255)
                                                                              aaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  1:2203970/8‑100 (MQ=37)
                                                                               aataagAAACGGCAGGATATTAGCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTCGCGTTCAGCAATAAAATCCAGCGCTttgttgn                       <  2:1948222/98‑2 (MQ=255)
                                                                               aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                      <  2:2319676/100‑1 (MQ=255)
                                                                                ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                     >  1:1673429/1‑100 (MQ=255)
                                                                                 taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    >  1:1116518/1‑100 (MQ=255)
                                                                                 taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    <  2:1704109/100‑1 (MQ=255)
                                                                                  gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                   >  1:213871/4‑100 (MQ=38)
                                                                                  aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCAGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    >  2:1634347/1‑99 (MQ=255)
                                                                                   aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                  >  2:2164203/1‑100 (MQ=255)
                                                                                    aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                 >  1:317176/1‑100 (MQ=255)
                                                                                     aaaCGGCAGGATATTATCTCCTTCCGTTTATCTTTTTACACGCTAATTACTGCTTTTTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                <  2:302688/100‑1 (MQ=255)
                                                                                      aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  1:786837/1‑100 (MQ=255)
                                                                                      aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               <  2:1106680/100‑1 (MQ=255)
                                                                                      aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  1:1976494/1‑100 (MQ=39)
                                                                                      aaCGGCAGGAGATAATATCCTGCCGTGTATCTTTTTACACGCTAATTACTGCTGACTCTCTCGTTCAGACACAAAATCCAGCTCTttgtcg                        >  2:1424384/1‑89 (MQ=19)
                                                                                      aaCGGCAGCAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               <  1:191279/100‑1 (MQ=255)
                                                                                       aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  2:1752675/1‑99 (MQ=39)
                                                                                        cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             >  1:1258859/1‑100 (MQ=255)
                                                                                         ggCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa            >  2:90052/1‑100 (MQ=255)
                                                                                              gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCTAGCGCTTTGTTGATACGCTCGATACTg         >  1:22222/1‑98 (MQ=255)
                                                                                              gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       <  2:1850172/100‑1 (MQ=255)
                                                                                                 attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt    >  2:1106223/1‑100 (MQ=255)
                                                                                                  ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc   >  1:917081/1‑100 (MQ=255)
                                                                                                  ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc   <  1:1881769/100‑1 (MQ=255)
                                                                                                   tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  <  1:1466626/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: