Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 15.5% 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,5550AT15.5% 325.1 / 48.3 129intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (48/60);  new base T (10/10);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 130 total aligned reads displayed.

GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAG  >  NC_000913/4227456‑4227652
                                                                                                   |                                                                                                 
acGTGACCGAAGTTGAGTGCTGGCTTGCACCGAATCTGGGGTATGAGGCGTATTGATTGACAAATCTGGCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                   >  2:2888380/2‑100 (MQ=255)
 cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                  <  1:2846740/100‑1 (MQ=255)
  gTTACCGAAGTTGAGCGCTGTCTGGTACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCGTACAACAAACAGGGCGCTCAATg                                                                                                 <  1:1462102/100‑1 (MQ=255)
  gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                 <  1:459207/100‑1 (MQ=255)
     aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                              >  2:362071/1‑100 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                             <  1:2758351/100‑1 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                             >  2:1968927/1‑100 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                             >  2:2038120/1‑100 (MQ=255)
        gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                           >  1:1079562/1‑100 (MQ=255)
        gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                           <  2:2798974/100‑1 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACTAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                        <  2:584125/100‑1 (MQ=255)
            ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                       <  2:1259288/100‑1 (MQ=255)
             tGAGCGCTGGCTGGCACCGAATCTGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                      <  2:2419579/100‑2 (MQ=39)
             tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                      <  1:1906076/100‑2 (MQ=39)
              gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACTGGTCGCTCAATGAGTGCCCTGtct                                                                                     <  2:109960/100‑1 (MQ=255)
               aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                    >  2:122169/1‑100 (MQ=255)
                gatctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                   <  2:147216/97‑5 (MQ=38)
                 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                                  <  2:1945996/100‑6 (MQ=38)
                  gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                 <  1:2093394/100‑1 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                 <  1:1978748/99‑1 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                 <  2:466285/99‑1 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                <  2:2070448/100‑1 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                >  1:2317883/1‑100 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                >  2:2599431/1‑100 (MQ=255)
                   ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                <  1:308536/100‑1 (MQ=255)
                    tggctggcATCGAATCTGGGGTATGACGCGGACTGATTCATAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                               <  2:2443520/100‑1 (MQ=255)
                     ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                              >  2:95793/1‑100 (MQ=255)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                             <  2:1554475/100‑11 (MQ=25)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                             <  1:1523994/100‑11 (MQ=25)
                       ctggcTCCGAATCTTGGGTATGACGCAAACTGATTCACAAATCTATCACTTTTCCTTACAACTAACATGTCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                            <  2:2644901/100‑1 (MQ=255)
                       ctgccatggAATCTTTGGTATGACGCGCACTGTTTCACCAATCTGTCACTTTTCCTTCCAACACTCAGTGAGTTCAATGATTGCCCTGGCTCTTGAGTaa                                                                            <  2:1790758/92‑1 (MQ=255)
                        tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                           <  1:2876068/100‑1 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                          <  2:2769145/100‑1 (MQ=255)
                          gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                         >  2:599255/1‑100 (MQ=255)
                          gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                         <  2:2366648/100‑1 (MQ=255)
                           cACCGAATCTGGTGTATGACGCGGACTGATTCACAACTCTGTCACTTTTCCTTACAACAACCCGGGCGCCCAATGAGTGCCCTGTCTCTTCATTAATATg                                                                        <  2:561312/100‑1 (MQ=255)
                           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                        >  1:1721133/1‑100 (MQ=255)
                           cACCAAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTACCCTGTCTCTTTATTAATATg                                                                        >  1:1373993/1‑100 (MQ=255)
                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                       <  1:459374/100‑1 (MQ=255)
                               gAATCTGGTGTATGACGCTGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                     <  1:2389693/99‑1 (MQ=255)
                               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                    <  1:367581/100‑1 (MQ=255)
                               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                    <  1:708913/100‑1 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAGACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                  <  2:1682225/100‑1 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                  <  2:8718/100‑1 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                  <  2:1579139/100‑1 (MQ=255)
                                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                 <  2:1301313/100‑1 (MQ=255)
                                   cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                                <  1:1459090/100‑1 (MQ=255)
                                   cTGGGGTATGACGCGGACTGATTCACAAATCTGGCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTACCCTGTCTCTTTATTAATATGAAACACtt                                                                >  2:1816737/1‑100 (MQ=255)
                                    tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                               >  2:1763646/1‑100 (MQ=255)
                                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGAGCACAAAGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                              >  1:2505635/1‑100 (MQ=255)
                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                           >  2:2517990/1‑100 (MQ=255)
                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                           <  1:2439178/100‑1 (MQ=255)
                                         tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                          <  1:2594856/100‑1 (MQ=255)
                                           tGACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                        <  1:648945/100‑1 (MQ=255)
                                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTggg                                                        >  1:774834/1‑99 (MQ=255)
                                              ctcttaCTGATTCACAAATTTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                       <  1:1037840/93‑1 (MQ=255)
                                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                    <  2:1079555/100‑1 (MQ=255)
                                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAATCAGGc                                                 <  2:1972503/100‑1 (MQ=255)
                                                  gACTGATTCACAAATCTGTCAATTTTTCTGACAACAAACAGGTCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAATACTTATACTGGAAACAGGc                                                 <  2:260543/100‑1 (MQ=255)
                                                      gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGCTCAATGAATGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                             >  1:2314791/1‑100 (MQ=255)
                                                       aTTCACAAATCTGTCACTTTTCCTTATATCATTCATTGTTTTCTATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                            <  1:2501731/100‑1 (MQ=255)
                                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGTGCTCAATGAGTGCCCTGTCCCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                            <  2:1648133/100‑1 (MQ=255)
                                                        ttCGCAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                           >  2:228607/1‑100 (MQ=255)
                                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGCGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAACCACTTATACTGGAAACAGGCTGGaat                                           >  2:1379427/1‑100 (MQ=255)
                                                          cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                         >  1:466989/1‑100 (MQ=255)
                                                           acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                        >  1:827775/1‑100 (MQ=255)
                                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGGAACAGGCTGGAATAATc                                       <  2:774830/100‑1 (MQ=255)
                                                              aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                     >  2:1678970/1‑100 (MQ=255)
                                                                tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGTCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                   <  1:2360359/100‑1 (MQ=255)
                                                                   gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                <  1:337533/100‑1 (MQ=255)
                                                                     cACTTTTCCTTACAACAAACAGGGCGCTCACTGAGTGCCCTGTCTCTTTATAAATATGAAACACTAATAATGGAAACAGGCTGGAATAATCTTAGACGGa                              >  1:631778/1‑100 (MQ=255)
                                                                     cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                              >  2:443136/1‑100 (MQ=255)
                                                                     cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                              <  2:2792936/100‑1 (MQ=255)
                                                                     cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCCTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                              >  2:161572/1‑100 (MQ=255)
                                                                       cTTTTCCTTACAACAAACAGTGCGTTCAATGATTTCTCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                            <  1:1039746/100‑1 (MQ=255)
                                                                       cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                            >  2:2357128/1‑100 (MQ=255)
                                                                          ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                         >  1:2196488/1‑100 (MQ=255)
                                                                          ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                         >  1:2582304/1‑100 (MQ=255)
                                                                          ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                         >  1:2613636/1‑100 (MQ=255)
                                                                           tCCTTACAACAAACATGGCGCTCAATGAGTGCCCTTTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                        <  1:1763666/100‑1 (MQ=255)
                                                                           tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                         >  1:340152/1‑99 (MQ=255)
                                                                            ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                       <  1:2115474/100‑1 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                      >  1:1815233/1‑100 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                      <  1:1513102/100‑1 (MQ=255)
                                                                              ttAGAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                     <  2:2376727/100‑1 (MQ=255)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGGAACACTTATACTGGAAACAGGCTGTAATAATCTTAGCCGGAAAAACGAATg                    <  2:150862/100‑1 (MQ=255)
                                                                                taaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATTAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                   >  2:2235846/8‑100 (MQ=37)
                                                                                  aagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGACACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                 >  2:1643112/6‑100 (MQ=38)
                                                                                        caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc           >  1:2559125/1‑100 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc           >  1:2552315/1‑100 (MQ=255)
                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac         >  2:943946/1‑100 (MQ=255)
                                                                                          gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca          >  1:1205675/1‑99 (MQ=255)
                                                                                           ggCGCTCAAGGAGTGCCCTGTCTCTTTATTTATATGAAACACTTATACTGTAAACAGGCTGGAATAATCTTATCCGGAAAATCGAATGTCAAAAGGcac         <  2:1080205/99‑1 (MQ=255)
                                                                                            gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg       <  2:65985/100‑1 (MQ=255)
                                                                                            gcacTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg       >  2:1492389/4‑100 (MQ=255)
                                                                                             cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt      >  1:890131/1‑100 (MQ=255)
                                                                                             cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt      >  1:2873620/1‑100 (MQ=255)
                                                                                             cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt      >  2:2830458/1‑100 (MQ=255)
                                                                                               cTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGAAAAAACGAATGTCAAAAGGCACAGTaa    >  2:727061/1‑100 (MQ=255)
                                                                                                 cAATGAGTGCCCTGTCTCTTTATTAATATGAAACCCTTATACTGGAAACAGGCTGGAATAATCTTAGCAGGAAAAACAAATGTCAAAAGGCACAGTAAAg  >  1:2558019/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAG  >  NC_000913/4227456‑4227652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: