Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 19.0% intergenic (+44/‑176) metH → / → yjbB homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,5550AT19.0% 415.5 / 81.9 177intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (74/69);  new base T (13/21);  total (87/90)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 179 total aligned reads displayed.

GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4227456‑4227653
                                                                                                   |                                                                                                  
gCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                    >  2:306647/1‑100 (MQ=255)
 cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAAACTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                   >  2:2360853/1‑100 (MQ=255)
  gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCATTACAACAAACAGGGCGCTCAATg                                                                                                  >  2:3096081/1‑100 (MQ=255)
    tACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg                                                                                                >  1:3614457/1‑100 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                              >  1:55470/1‑100 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                              >  2:2828055/1‑100 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         <  2:3730319/100‑1 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         >  1:882555/1‑100 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                         <  1:367203/100‑1 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                         <  1:2270182/100‑1 (MQ=255)
            ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                        <  1:1259828/100‑1 (MQ=255)
            ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                        <  1:812603/100‑1 (MQ=255)
            ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                        <  1:2310965/100‑1 (MQ=255)
             tGAGGGCTGGCTGGCACCGACTCTGTGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTGACAACGAACAGGGCACTCATTGAGCGCCCt                                                                                          <  2:1884294/97‑1 (MQ=255)
             tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                       >  2:2089288/1‑100 (MQ=255)
             tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                       <  2:3083145/100‑1 (MQ=255)
              gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                      <  1:1952494/100‑1 (MQ=255)
                gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAGACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                    <  1:418100/100‑5 (MQ=37)
                gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                    <  1:3273318/100‑5 (MQ=38)
                 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                                   <  1:568261/100‑6 (MQ=38)
                  gctggcttGCACCGAGTCTGGGGTATGACGCGTTCTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                  <  2:570660/100‑7 (MQ=25)
                  gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                  <  1:2803771/100‑1 (MQ=255)
                  gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                  <  2:1845911/100‑7 (MQ=38)
                     ggctggcACCTGATCTGGGGTATGACGCGGACTGATTCATAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaa                                                                               <  2:724709/100‑10 (MQ=25)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatta                                                                              <  2:1790489/100‑1 (MQ=255)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatta                                                                              <  1:3920970/100‑1 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           >  2:537093/1‑100 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           >  1:1091520/1‑100 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           >  2:2311576/1‑100 (MQ=255)
                          gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           <  2:822010/99‑1 (MQ=255)
                           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                         >  1:780177/1‑100 (MQ=255)
                           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                         >  1:354261/1‑100 (MQ=255)
                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                        >  2:3661842/1‑100 (MQ=255)
                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                        >  1:2321439/1‑100 (MQ=255)
                               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                     <  1:1709230/100‑1 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                   >  2:1804976/1‑100 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                   >  2:121903/1‑100 (MQ=255)
                                  tCTGGGTTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                  <  2:3135589/100‑1 (MQ=255)
                                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                  <  2:2185330/100‑1 (MQ=255)
                                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                  >  1:3162776/1‑100 (MQ=255)
                                    tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                >  1:743888/1‑100 (MQ=255)
                                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                               >  1:2092818/1‑100 (MQ=255)
                                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                             <  1:2763134/100‑1 (MQ=255)
                                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                             >  1:3187987/1‑100 (MQ=255)
                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTGTCCTTACAACAAACAGGGCGCTCAATGAGTGCGCTGTCTCTTTATTAATATGAAACACTTATACt                                                            <  2:2739812/100‑1 (MQ=255)
                                          aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                          <  2:647107/100‑1 (MQ=255)
                                          aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                          >  2:1450745/1‑100 (MQ=255)
                                           tGACGCGGACTGATTCACAAATCTGTCACTTTTTCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                         >  2:2568316/1‑100 (MQ=255)
                                            gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTCCCACAAACAGGGCGCTCCATGAGTGCCCTGTCTCTTTATTCATATGAAACACTTATACCGGa                                                         <  2:1164430/99‑1 (MQ=255)
                                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                       >  1:2543560/1‑100 (MQ=255)
                                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                       <  2:1857944/100‑1 (MQ=255)
                                              cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                      >  1:1549191/1‑100 (MQ=255)
                                              cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                      >  1:1632110/1‑100 (MQ=255)
                                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                     <  2:3851546/100‑1 (MQ=255)
                                                ctgACTCATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATAGGAAACACTTATACTGGAAACAg                                                    <  2:3684422/98‑1 (MQ=255)
                                                cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                                    >  2:257589/1‑100 (MQ=255)
                                                  gACTGATTCACAAATCTGTCACTTTTCCTTACTACAAAAAGTGCGCCCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                                  <  2:1958084/100‑1 (MQ=255)
                                                     tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                               <  1:3100099/100‑1 (MQ=255)
                                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                            >  1:1302465/1‑100 (MQ=255)
                                                          cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGATCAATAACTGCCCTGTCTCTTCATCAATATGAAACACACATACTGGAAACACGCTCGaaaa                                           >  2:986213/1‑97 (MQ=255)
                                                           acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                         >  1:1766826/1‑100 (MQ=255)
                                                           acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                         <  1:3058182/100‑1 (MQ=255)
                                                             aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                       <  2:3369064/100‑1 (MQ=255)
                                                               aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATGAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                     <  2:1692520/100‑1 (MQ=255)
                                                                 cTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                                   >  2:2007578/1‑100 (MQ=255)
                                                                   gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                 <  1:1566955/100‑1 (MQ=255)
                                                                   gTCAATTTTCCTTACAACAAACAGGGTGCTAAATGAGTGCCCTGTCGCTTTAGTAATATGAAACACTTATACTGGAAGCAGGCTGGAATAATCTTa                                     <  2:3913403/96‑1 (MQ=255)
                                                                       cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                             <  2:1537746/100‑1 (MQ=255)
                                                                            ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                        >  1:997550/1‑100 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                       >  2:3124415/1‑100 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                       >  2:1322499/1‑100 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                       >  2:1778259/1‑100 (MQ=255)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                     <  2:3405765/100‑1 (MQ=255)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                     >  1:20259/1‑100 (MQ=255)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                     >  1:1624717/1‑100 (MQ=255)
                                                                                 caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   >  1:2708614/1‑100 (MQ=255)
                                                                                 caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   >  2:712506/1‑100 (MQ=255)
                                                                                 caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   <  2:3887379/100‑1 (MQ=255)
                                                                                 aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   >  1:621675/7‑100 (MQ=37)
                                                                                   acacacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTACTATGAAACACTGATACTGGACACAGGCTGGAATAATCTGAGCCGGAAAAACGAATGTCaa                 >  2:3474881/1‑100 (MQ=255)
                                                                                   acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTTTTTTTTTTTTTCACTTATACTGGAAACAGGCTGGAATAATTTTAGCCGGAAAAACGAATGTCaa                 >  2:3215561/1‑100 (MQ=255)
                                                                                   acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                 <  1:941786/100‑1 (MQ=255)
                                                                                    gagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                >  2:3168954/4‑100 (MQ=38)
                                                                                    caaactGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                 <  1:1197628/99‑1 (MQ=255)
                                                                                     agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa               >  1:14021/3‑100 (MQ=38)
                                                                                      gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg              >  2:1250159/2‑100 (MQ=255)
                                                                                      aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg              >  1:3440980/1‑100 (MQ=255)
                                                                                       acaGGGCGCTCAATGAGTGCCCTGTCTCTTTAGTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg             <  2:3689762/100‑1 (MQ=255)
                                                                                        caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc            <  1:1804985/100‑1 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg             >  1:1299674/1‑99 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc            >  1:2173828/1‑100 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCGGGAATAATCTTAGCCGCAGAAACGAATGTCAAAAGGc            >  2:2148379/1‑100 (MQ=255)
                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac          >  2:1174551/1‑100 (MQ=255)
                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac          <  2:882551/100‑1 (MQ=255)
                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac          <  1:114541/100‑1 (MQ=255)
                                                                                          gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac          >  2:165542/1‑100 (MQ=255)
                                                                                              gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa      >  1:2439463/1‑100 (MQ=255)
                                                                                              gcTCAATGAGTGCCCTGTCTCTTTATTAAAAAGAAACACTTATACTGGAAACAGGCTGGAATAAACTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt       <  2:1522117/99‑1 (MQ=255)
                                                                                                 cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaacg   >  2:2392581/1‑98 (MQ=255)
                                                                                                  aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGa  <  1:3173779/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4227456‑4227653

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: