Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 15.3% intergenic (+39/‑181) 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,5500GA15.3% 286.7 / 34.1 123intergenic (+39/‑181)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (45/58);  new base A (9/10);  total (54/69)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 133 total aligned reads displayed.

TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGT  >  NC_000913/4227451‑4227648
                                                                                                   |                                                                                                  
tATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg                                                                                                    >  2:308533/1‑100 (MQ=255)
 aTGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc                                                                                                   <  2:2939698/100‑1 (MQ=255)
  tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCt                                                                                                  >  1:2632799/1‑100 (MQ=255)
     acGTGACCGAAGTTGAGTGCTGGCTTGCACCGAATCTGGGGTATGAGGCGTATTGATTGACAAATCTGGCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                               >  2:2888380/2‑100 (MQ=255)
      cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                              <  1:2846740/100‑1 (MQ=255)
       gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                             <  1:459207/100‑1 (MQ=255)
          aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                          >  2:362071/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                         >  2:1968927/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                         >  1:993334/1‑100 (MQ=255)
             gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                       >  1:1079562/1‑100 (MQ=255)
             gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                       <  2:2798974/100‑1 (MQ=255)
               aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                     <  2:1411152/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACTAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                    <  2:584125/100‑1 (MQ=255)
                  tGAGCGCTGGCTGGCACCGAATCTGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                  <  2:2419579/100‑2 (MQ=39)
                   gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACTGGTCGCTCAATGAGTGCCCTGtct                                                                                 <  2:109960/100‑1 (MQ=255)
                   gAGCGCTGGATGGCACCGACTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCGCAATGAGTGCCCTGgct                                                                                 >  2:1762572/1‑100 (MQ=255)
                    aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                >  2:122169/1‑100 (MQ=255)
                      cgcTGGCTGGCACCGAATCTGGGGTACGACGCCGACTGATTCACCAATCGGTCACTTTTCCTTATAACAAACAGGGCGCTCAATGAGTGCCCTGtctcct                                                                              <  2:600142/100‑3 (MQ=255)
                       gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                             <  1:2093394/100‑1 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                             <  1:1978748/99‑1 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                             <  2:466285/99‑1 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                            >  2:2599431/1‑100 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                            <  1:308536/100‑1 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                            <  2:2070448/100‑1 (MQ=255)
                         tggctggcATCGAATCTGGGGTATGACGCGGACTGATTCATAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                           <  2:2443520/100‑1 (MQ=255)
                          ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                          >  2:95793/1‑100 (MQ=255)
                          ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtagatc                                                                          >  1:147218/1‑91 (MQ=38)
                           gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                         <  2:1554475/100‑11 (MQ=25)
                           gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                         <  1:1723931/100‑11 (MQ=37)
                           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)
                              ggcACCGAATCTGGAGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                      <  2:2769646/100‑1 (MQ=255)
                               gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                     <  2:2366648/100‑1 (MQ=255)
                                cACCGAATCTGGTGTATGACGCGGACTGATTCACAACTCTGTCACTTTTCCTTACAACAACCCGGGCGCCCAATGAGTGCCCTGTCTCTTCATTAATATg                                                                    <  2:561312/100‑1 (MQ=255)
                                cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCATGTCTCTTTATTAATATg                                                                    >  1:761955/1‑100 (MQ=255)
                                cACCAAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTACCCTGTCTCTTTATTAATATg                                                                    >  1:1373993/1‑100 (MQ=255)
                                 aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                   <  1:459374/100‑1 (MQ=255)
                                  ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                  >  2:1136380/1‑100 (MQ=255)
                                    gAATCTGGTGTATGACGCTGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                 <  1:2389693/99‑1 (MQ=255)
                                    gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:367581/100‑1 (MQ=255)
                                    gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:2568222/100‑1 (MQ=255)
                                      aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAGACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                              <  2:1682225/100‑1 (MQ=255)
                                      aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                              <  2:8718/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)
                                          ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGAGCACAAAGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                          >  1:2505635/1‑100 (MQ=255)
                                             gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                       >  2:2517990/1‑100 (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)
                                                       gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                             >  1:1253717/1‑100 (MQ=255)
                                                       gACTGATTCACAAATCTGTCAATTTTTCTGACAACAAACAGGTCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAATACTTATACTGGAAACAGGc                                             <  2:260543/100‑1 (MQ=255)
                                                           gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGCTCAATGAATGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                         >  1:2314791/1‑100 (MQ=255)
                                                            aTTCACAAATCTGTCACTTTTCCTTATATCATTCATTGTTTTCTATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                        <  1:2501731/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:2792936/100‑1 (MQ=255)
                                                                          cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                          >  2:443136/1‑100 (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:2613636/1‑100 (MQ=255)
                                                                               ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                     >  1:2582304/1‑100 (MQ=255)
                                                                               ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTCGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                     <  2:1987335/100‑1 (MQ=255)
                                                                                tCCTTACAACAAACATGGCGCTCAATGAGTGCCCTTTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                    <  1:1763666/100‑1 (MQ=255)
                                                                                tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                     >  1:340152/1‑99 (MQ=255)
                                                                                 ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                   >  2:1536240/1‑100 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                  <  1:1513102/100‑1 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                  >  1:1815233/1‑100 (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)
                                                                                        acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa            >  2:1903055/1‑100 (MQ=255)
                                                                                        acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa            >  1:2750000/1‑100 (MQ=255)
                                                                                          agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          >  1:2366493/3‑100 (MQ=38)
                                                                                           aacaGGGCGCTCAATGAGAGCCCTGTCTGTTTATTTATATGAATCATTTATACTTGAAACAGGCTGGAATAATTTTTGCCTGAAAAACGAATGTCAAAAg         <  2:1639661/100‑1 (MQ=255)
                                                                                             caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc       >  1:2559125/1‑100 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgac       >  2:2847605/1‑98 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc       >  1:2066267/1‑100 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc       >  1:2290184/1‑100 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc       >  1:2509332/1‑100 (MQ=255)
                                                                                               gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac     >  2:943946/1‑100 (MQ=255)
                                                                                               gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca      >  1:1205675/1‑99 (MQ=255)
                                                                                                 gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg   <  2:65985/100‑1 (MQ=255)
                                                                                                  cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt  >  1:2873620/1‑100 (MQ=255)
                                                                                                  cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt  >  1:890131/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGT  >  NC_000913/4227451‑4227648

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: