Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,560 T→C 15.0% intergenic (+49/‑171) 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,5600TC15.0% 300.0 / 44.1 135intergenic (+49/‑171)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base T (49/65);  new base C (10/10);  total (59/76)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 138 total aligned reads displayed.

ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAG  >  NC_000913/4227461‑4227659
                                                                                                   |                                                                                                   
 ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                    >  2:2038120/1‑100 (MQ=255)
 ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                    <  1:2758351/100‑1 (MQ=255)
   gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                                  >  1:1079562/1‑100 (MQ=255)
   gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                                  <  2:2798974/100‑1 (MQ=255)
    aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                 <  2:2415652/100‑1 (MQ=255)
     aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                <  2:1411152/100‑1 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACTAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                               <  2:584125/100‑1 (MQ=255)
       ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                              <  1:2915605/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)
         gAGCGCTGGATGGCACCGACTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCGCAATGAGTGCCCTGgct                                                                                            >  2:1762572/1‑100 (MQ=255)
          aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                           >  2:122169/1‑100 (MQ=255)
            cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                                         <  2:1945996/100‑6 (MQ=38)
            cgcTGGCTGGCACCGAATCTGGGGTACGACGCCGACTGATTCACCAATCGGTCACTTTTCCTTATAACAAACAGGGCGCTCAATGAGTGCCCTGtctcct                                                                                         <  2:600142/100‑3 (MQ=255)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                        <  2:466285/99‑1 (MQ=255)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                        <  1:1978748/99‑1 (MQ=255)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                       <  1:308536/100‑1 (MQ=255)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                       >  2:2599431/1‑100 (MQ=255)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                       >  1:2317883/1‑100 (MQ=255)
                ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                     >  2:95793/1‑100 (MQ=255)
                 gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                                    <  1:1523994/100‑11 (MQ=25)
                 gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                                    <  1:1723931/100‑11 (MQ=37)
                  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)
                    ggcACCGAATCTGGAGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                 <  2:2769646/100‑1 (MQ=255)
                     gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                                <  2:2366648/100‑1 (MQ=255)
                      cACCGAATCTGGTGTATGACGCGGACTGATTCACAACTCTGTCACTTTTCCTTACAACAACCCGGGCGCCCAATGAGTGCCCTGTCTCTTCATTAATATg                                                                               <  2:561312/100‑1 (MQ=255)
                      cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                               >  1:1721133/1‑100 (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:708913/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)
                               tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                      >  2:1763646/1‑100 (MQ=255)
                                ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGAGCACAAAGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                                     >  1:2505635/1‑100 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  <  1:2439178/100‑1 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  >  2:2517990/1‑100 (MQ=255)
                                    tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                                 <  1:2594856/100‑1 (MQ=255)
                                      tGACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                               <  1:648945/100‑1 (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)
                                                 gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGCTCAATGAATGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                                    >  1:2314791/1‑100 (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)
                                                       caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                              <  2:826168/100‑1 (MQ=255)
                                                         aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                            >  2:1678970/1‑100 (MQ=255)
                                                              gTCACTTTTCCTTACAACAAACAGGGCGTTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                       <  1:347697/100‑1 (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: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                              <  1:2115474/100‑1 (MQ=255)
                                                                       ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                              >  2:1536240/1‑100 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             <  1:1513102/100‑1 (MQ=255)
                                                                         ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                            <  2:464605/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                       >  1:2750000/1‑100 (MQ=255)
                                                                              acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                       >  2:1903055/1‑100 (MQ=255)
                                                                                agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa                     >  1:2366493/3‑100 (MQ=38)
                                                                                   caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:2559125/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:2290184/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:2066267/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:2509332/1‑100 (MQ=255)
                                                                                     gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                >  2:943946/1‑100 (MQ=255)
                                                                                      ggCGCTCAAGGAGTGCCCTGTCTCTTTATTTATATGAAACACTTATACTGTAAACAGGCTGGAATAATCTTATCCGGAAAATCGAATGTCAAAAGGcac                <  2:1080205/99‑1 (MQ=255)
                                                                                       gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg              <  2:65985/100‑1 (MQ=255)
                                                                                       gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg              >  1:465177/1‑100 (MQ=255)
                                                                                       gcacTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg              >  2:1492389/4‑100 (MQ=255)
                                                                                        cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt             >  2:2830458/1‑100 (MQ=255)
                                                                                        cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt             >  1:890131/1‑100 (MQ=255)
                                                                                            cAATGAGTGCCCTGTCTCTTTATTAATATGAAACCCTTATACTGGAAACAGGCTGGAATAATCTTAGCAGGAAAAACAAATGTCAAAAGGCACAGTAAAg         >  1:2558019/1‑100 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  2:520707/100‑1 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  1:2830490/100‑1 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      >  2:2111707/1‑100 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  1:2558895/100‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGTCTGGAATAAAAAAAAAAAAAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc    <  2:2929815/100‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc    <  2:340149/100‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc    <  2:841910/100‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGc     <  1:2098263/99‑1 (MQ=255)
                                                                                                   tGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg  <  2:2386361/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAG  >  NC_000913/4227461‑4227659

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: