Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,228 C→A 16.2% intergenic (+23/+29) 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,2280CA16.2% 236.7 / 33.6 91intergenic (+23/+29)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base A (6/9);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTA  >  NC_000913/2519129‑2519326
                                                                                                   |                                                                                                  
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:1241728/100‑1 (MQ=255)
                               aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                     <  1:1619795/100‑1 (MQ=255)
                                 ttCGCAATTTCTTCGTCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                   <  2:569799/100‑1 (MQ=255)
                                  tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAACTAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                  >  2:1267438/1‑100 (MQ=255)
                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                 >  2:322725/1‑100 (MQ=255)
                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                 >  1:1526363/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:1740432/100‑1 (MQ=255)
                                                  cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                  <  1:1294645/100‑1 (MQ=255)
                                                  cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                  <  2:308648/100‑1 (MQ=255)
                                                   gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                                 >  1:427504/1‑100 (MQ=255)
                                                    ttttgATATTGATGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGTGttta                                                  <  2:1213939/94‑3 (MQ=255)
                                                    tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                <  2:1052923/100‑1 (MQ=255)
                                                    tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                <  1:1615774/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                      >  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                 <  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)
                                                                                              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)
                                                                                                  tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa  <  2:1914702/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTA  >  NC_000913/2519129‑2519326

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: