Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 15.8% 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,2220TG15.8% 203.8 / 35.8 89intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base G (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 = 4.92e-01

AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCG  >  NC_000913/2519123‑2519319
                                                                                                   |                                                                                                 
aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGACCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGat                                                                                                   >  2:1162820/1‑100 (MQ=255)
 ncAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                  >  1:1948236/2‑100 (MQ=255)
  cAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAGAAAACGGCAGGatat                                                                                                 >  1:2520434/1‑100 (MQ=255)
   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:1338102/1‑100 (MQ=255)
           gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                        >  2:2293284/1‑100 (MQ=255)
           gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                        >  1:1164850/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                                                                                   <  1:815875/100‑1 (MQ=39)
                gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                   <  2:814274/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:1639751/100‑1 (MQ=255)
                          tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                         <  1:363556/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                                                                    >  2:2496900/1‑100 (MQ=255)
                               tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                    <  1:2236482/100‑1 (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:1845900/100‑1 (MQ=255)
                                                           ctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                        <  1:577930/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               <  2:1704109/100‑1 (MQ=255)
                                                                                    taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc               >  1:1116518/1‑100 (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          <  2:1106680/100‑1 (MQ=255)
                                                                                         aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg          >  1:786837/1‑100 (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)
                                                                                                   |                                                                                                 
AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCG  >  NC_000913/2519123‑2519319

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: