Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,560 T→C 18.9% 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,5600TC18.9% 381.0 / 80.6 179intergenic (+49/‑171)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base T (75/70);  new base C (12/22);  total (87/92)
Fisher's exact test for biased strand distribution p-value = 9.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 179 total aligned reads displayed.

CCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCA  >  NC_000913/4227462‑4227658
                                                                                                  |                                                                                                  
ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                   >  2:2828055/1‑100 (MQ=255)
ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                   >  1:55470/1‑100 (MQ=255)
ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                   <  2:3691242/100‑1 (MQ=255)
   aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                <  2:3686721/100‑1 (MQ=255)
    aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                               <  1:3310433/100‑1 (MQ=255)
     gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                              >  1:882555/1‑100 (MQ=255)
     gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                              <  1:367203/100‑1 (MQ=255)
      ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                             <  1:1259828/100‑1 (MQ=255)
      ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                             <  1:2310965/100‑1 (MQ=255)
      ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                             <  1:2208373/100‑1 (MQ=255)
      ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                             <  1:812603/100‑1 (MQ=255)
       tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                            <  2:3083145/100‑1 (MQ=255)
       tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                            >  2:2089288/1‑100 (MQ=255)
        gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                           <  1:1952494/100‑1 (MQ=255)
        gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                           <  1:3931917/100‑1 (MQ=255)
          gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAGACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                         <  1:418100/100‑5 (MQ=37)
            gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                       <  1:2803771/100‑1 (MQ=255)
            gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                       <  2:1845911/100‑7 (MQ=38)
             ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                      <  2:3336699/100‑1 (MQ=255)
             ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                      <  2:3132467/100‑1 (MQ=255)
               ggctggcACCTGATCTGGGGTATGACGCGGACTGATTCATAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaa                                                                                    <  2:724709/100‑10 (MQ=25)
                gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTGCAACAAATAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                                   <  2:2500209/100‑11 (MQ=25)
                gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatta                                                                                   <  1:3920970/100‑1 (MQ=255)
                 ctggcACCGAATTTGGGGTATGACGCGGACTGGTTCACGACTCTGTCACTTTTCCGTACAACAAGCAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                       <  2:3457932/95‑7 (MQ=21)
                   ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                >  1:1091520/1‑100 (MQ=255)
                    gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                               >  2:1080498/1‑100 (MQ=255)
                    gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                <  2:822010/99‑1 (MQ=255)
                     cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                              >  1:780107/1‑100 (MQ=255)
                     cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                              >  1:354261/1‑100 (MQ=255)
                      aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                             >  2:3661842/1‑100 (MQ=255)
                       ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                            >  1:3250185/1‑100 (MQ=255)
                         gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                          <  1:1709230/100‑1 (MQ=255)
                           aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                        >  2:121903/1‑100 (MQ=255)
                            tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                       >  1:3162776/1‑100 (MQ=255)
                                 ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                  <  1:2763134/100‑1 (MQ=255)
                                  gTATGACGCGGACTGATTCACAAATCTGTCACTTGTCCTTACAACAAACAGGGCGCTCAATGAGTGCGCTGTCTCTTTATTAATATGAAACACTTATACt                                                                 <  2:2739812/100‑1 (MQ=255)
                                   tATGACGCGGACTTATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                                <  1:693119/100‑1 (MQ=255)
                                       aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                            <  2:1857944/100‑1 (MQ=255)
                                       aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                            >  1:2543560/1‑100 (MQ=255)
                                        cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                           >  1:1632110/1‑100 (MQ=255)
                                        cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                           >  1:1549191/1‑100 (MQ=255)
                                          cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                                         >  2:257589/1‑100 (MQ=255)
                                              cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                                     <  1:209792/100‑1 (MQ=255)
                                              cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                                     <  2:2018851/100‑1 (MQ=255)
                                               tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                    <  1:3100099/100‑1 (MQ=255)
                                                gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATCTGAAACACTTATACTGGAAACAGGCt                                                      >  2:1893376/1‑97 (MQ=255)
                                                gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                                   >  2:126662/1‑100 (MQ=255)
                                                  ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                 >  2:550760/1‑100 (MQ=255)
                                                  ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                 >  1:1302465/1‑100 (MQ=255)
                                                  ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                 >  2:579018/1‑100 (MQ=255)
                                                   tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                                >  2:3399609/1‑100 (MQ=255)
                                                    cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                                >  2:3756099/1‑99 (MQ=255)
                                                    cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGATCAATAACTGCCCTGTCTCTTCATCAATATGAAACACACATACTGGAAACACGCTCGaaaa                                                >  2:986213/1‑97 (MQ=255)
                                                     acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                              <  1:3058182/100‑1 (MQ=255)
                                                     acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                              >  1:1766826/1‑100 (MQ=255)
                                                      caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                             <  2:3859531/100‑1 (MQ=255)
                                                      caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                             >  2:2335845/1‑100 (MQ=255)
                                                         aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATGAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                          <  2:1692520/100‑1 (MQ=255)
                                                            tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                       <  1:3180773/100‑1 (MQ=255)
                                                             gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                      <  1:1566955/100‑1 (MQ=255)
                                                                 cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                                  <  2:1537746/100‑1 (MQ=255)
                                                                       cgtACAACAAACAAGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                            >  1:2654766/3‑100 (MQ=255)
                                                                       cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                            >  2:1778259/1‑100 (MQ=255)
                                                                       cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                            <  2:1796429/100‑1 (MQ=255)
                                                                         tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                          >  1:1624717/1‑100 (MQ=255)
                                                                           caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                        <  2:3887379/100‑1 (MQ=255)
                                                                           caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                        >  2:712506/1‑100 (MQ=255)
                                                                           aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATCAAACACTTATACTGGAAACAGGCTGGAATACTCTTCGCCGGAAAAACGAATGTc                        >  2:2197327/7‑100 (MQ=25)
                                                                            aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                       >  2:1824434/1‑100 (MQ=255)
                                                                            aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                       <  1:962416/100‑1 (MQ=255)
                                                                             acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                      <  1:941786/100‑1 (MQ=255)
                                                                              gagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAGAACGAATGTCaaa                     >  1:2292287/4‑100 (MQ=38)
                                                                              gagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                     >  2:3168954/4‑100 (MQ=38)
                                                                              caaactGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                      <  1:1197628/99‑1 (MQ=255)
                                                                              caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                     >  2:1372798/1‑100 (MQ=255)
                                                                               agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa                    >  1:14021/3‑100 (MQ=38)
                                                                                gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                   >  2:1250159/2‑100 (MQ=255)
                                                                                gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                   >  2:3157787/2‑100 (MQ=39)
                                                                                aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                   >  2:46744/1‑100 (MQ=255)
                                                                                aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                   >  1:3440980/1‑100 (MQ=255)
                                                                                  caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                 <  1:1804985/100‑1 (MQ=255)
                                                                                  caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg                  >  1:1299674/1‑99 (MQ=255)
                                                                                  caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg                  >  2:3626504/1‑99 (MQ=255)
                                                                                  caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCGGGAATAATCTTAGCCGCAGAAACGAATGTCAAAAGGc                 >  2:2148379/1‑100 (MQ=255)
                                                                                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               >  2:1174551/1‑100 (MQ=255)
                                                                                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               <  2:882551/100‑1 (MQ=255)
                                                                                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               >  2:308001/1‑100 (MQ=255)
                                                                                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               >  1:604279/1‑100 (MQ=255)
                                                                                    gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               <  1:114541/100‑1 (MQ=255)
                                                                                    gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac               >  2:165542/1‑100 (MQ=255)
                                                                                      gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAACGGCACAg             <  1:2074922/100‑1 (MQ=255)
                                                                                      gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg             <  2:1974269/100‑1 (MQ=255)
                                                                                        gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa           >  1:2439463/1‑100 (MQ=255)
                                                                                        gcTCAATGAGTGCCCTGTCTCTTTATTAAAAAGAAACACTTATACTGGAAACAGGCTGGAATAAACTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt            <  2:1522117/99‑1 (MQ=255)
                                                                                           cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaacg        >  2:2392581/1‑98 (MQ=255)
                                                                                              tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTTGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg     <  1:2335857/100‑1 (MQ=255)
                                                                                              tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg     >  2:3115226/1‑100 (MQ=255)
                                                                                              tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg     >  1:3086620/1‑100 (MQ=255)
                                                                                                aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc   <  1:1838451/100‑1 (MQ=255)
                                                                                                 gTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCa  >  1:737662/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCA  >  NC_000913/4227462‑4227658

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: