Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 19.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,5500GA19.3% 391.5 / 76.6 175intergenic (+39/‑181)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (73/68);  new base A (13/21);  total (86/89)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 184 total aligned reads displayed.

TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227451‑4227649
                                                                                                   |                                                                                                   
   gAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc                                                                                                  >  2:941781/1‑100 (MQ=255)
     gCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                >  2:306647/1‑100 (MQ=255)
      cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAAACTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                               >  2:2360853/1‑100 (MQ=255)
       gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCATTACAACAAACAGGGCGCTCAATg                                                                                              >  2:3096081/1‑100 (MQ=255)
         tACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg                                                                                            >  2:693115/1‑100 (MQ=255)
         tACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg                                                                                            >  1:3614457/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          >  1:55470/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          >  2:2828055/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAgcg                                                                                          >  2:1657685/1‑98 (MQ=255)
              aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGccc                                                                                       <  1:1657694/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     <  2:3730319/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                     <  1:2270182/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                     <  1:367203/100‑1 (MQ=255)
                 ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                    <  1:1259828/100‑1 (MQ=255)
                  tGAGGGCTGGCTGGCACCGACTCTGTGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTGACAACGAACAGGGCACTCATTGAGCGCCCt                                                                                      <  2:1884294/97‑1 (MQ=255)
                  tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                   <  2:3083145/100‑1 (MQ=255)
                  tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                   >  1:1374738/1‑100 (MQ=255)
                     gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAGACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                <  1:418100/100‑5 (MQ=37)
                     gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                <  1:3273318/100‑5 (MQ=38)
                     gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                <  1:909635/100‑5 (MQ=38)
                      cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                               <  2:1534637/100‑1 (MQ=255)
                      cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                               <  1:568261/100‑6 (MQ=38)
                       gctggcttGCACCGAGTCTGGGGTATGACGCGTTCTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                              <  2:570660/100‑7 (MQ=25)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                             <  2:3336699/100‑1 (MQ=255)
                        ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTCCCCTTACAACAAACAGGGCGCTCCATGAGCGCCCTGTCTCttt                                                                              <  2:176732/99‑1 (MQ=255)
                         tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                            <  1:1798621/100‑1 (MQ=255)
                           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)
                              ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                       >  2:537093/1‑100 (MQ=255)
                               gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                      >  2:1080498/1‑100 (MQ=255)
                                cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                     <  1:872161/100‑1 (MQ=255)
                                cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                     >  1:780177/1‑100 (MQ=255)
                                 aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                    >  1:2321439/1‑100 (MQ=255)
                                      aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                               >  1:1408488/1‑100 (MQ=255)
                                       tCTGGGTTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                              <  2:3135589/100‑1 (MQ=255)
                                       tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                              <  2:2185330/100‑1 (MQ=255)
                                       tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                              >  1:3162776/1‑100 (MQ=255)
                                       tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                              <  2:1646771/100‑1 (MQ=255)
                                         tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                            >  1:743888/1‑100 (MQ=255)
                                          ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                           >  1:2092818/1‑100 (MQ=255)
                                            ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                         >  1:3187987/1‑100 (MQ=255)
                                             gTATGACGCGGACTGATTCACAAATCTGTCACTTGTCCTTACAACAAACAGGGCGCTCAATGAGTGCGCTGTCTCTTTATTAATATGAAACACTTATACt                                                        <  2:2739812/100‑1 (MQ=255)
                                               aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                      <  2:647107/100‑1 (MQ=255)
                                               aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                      >  2:1450745/1‑100 (MQ=255)
                                                 gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                    <  1:3449636/100‑1 (MQ=255)
                                                  aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                   >  1:2543560/1‑100 (MQ=255)
                                                   cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                  >  1:1549191/1‑100 (MQ=255)
                                                     ctgACTCATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATAGGAAACACTTATACTGGAAACAg                                                <  2:3684422/98‑1 (MQ=255)
                                                     cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                                >  2:257589/1‑100 (MQ=255)
                                                       gACTGATTCACAAATCTGTCACTTTTCCTTACTACAAAAAGTGCGCCCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                              <  2:1958084/100‑1 (MQ=255)
                                                         cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                            <  2:2018851/100‑1 (MQ=255)
                                                          tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  1:3100099/100‑1 (MQ=255)
                                                           gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                          >  2:126662/1‑100 (MQ=255)
                                                              tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                       >  2:3399609/1‑100 (MQ=255)
                                                               cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                      >  2:1202176/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)
                                                                 caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                    <  2:781551/100‑1 (MQ=255)
                                                                  aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                   <  2:2092807/100‑1 (MQ=255)
                                                                  aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                   <  2:3369064/100‑1 (MQ=255)
                                                                        gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                             <  1:1566955/100‑1 (MQ=255)
                                                                        gTCAATTTTCCTTACAACAAACAGGGTGCTAAATGAGTGCCCTGTCGCTTTAGTAATATGAAACACTTATACTGGAAGCAGGCTGGAATAATCTTa                                 <  2:3913403/96‑1 (MQ=255)
                                                                         tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTGATACTGGAAACAGGCTGGGATAATCTTAGCCgg                            <  2:1850564/100‑1 (MQ=255)
                                                                             ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaa                        <  1:1823757/100‑1 (MQ=255)
                                                                               ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                      <  2:834597/100‑1 (MQ=255)
                                                                                tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                     <  1:2924280/100‑1 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   >  2:1778259/1‑100 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   >  2:1322499/1‑100 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   <  2:1796429/100‑1 (MQ=255)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   >  2:3124415/1‑100 (MQ=255)
                                                                                    tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGTAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 <  1:2782529/100‑1 (MQ=255)
                                                                                    tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 >  1:20259/1‑100 (MQ=255)
                                                                                    tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 >  1:1624717/1‑100 (MQ=255)
                                                                                      caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               >  2:712506/1‑100 (MQ=255)
                                                                                      caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               >  1:2708614/1‑100 (MQ=255)
                                                                                      aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               >  1:621675/7‑100 (MQ=37)
                                                                                      aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATCAAACACTTATACTGGAAACAGGCTGGAATACTCTTCGCCGGAAAAACGAATGTc               >  2:2197327/7‑100 (MQ=25)
                                                                                       aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa              <  1:962416/100‑1 (MQ=255)
                                                                                       aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa              >  2:1824434/1‑100 (MQ=255)
                                                                                       aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa              >  2:1248895/1‑100 (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)
                                                                                           aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          >  1:3440980/1‑100 (MQ=255)
                                                                                            acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         <  1:1740956/100‑1 (MQ=255)
                                                                                            acaGGGCGCTCAATGAGTGCCCTGTCTCTTTAGTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         <  2:3689762/100‑1 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         >  2:3626504/1‑99 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         >  1:1299674/1‑99 (MQ=255)
                                                                                             caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  1:2173828/1‑100 (MQ=255)
                                                                                               gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac      >  2:1174551/1‑100 (MQ=255)
                                                                                               gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac      <  1:114541/100‑1 (MQ=255)
                                                                                                 gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg    <  2:1974269/100‑1 (MQ=255)
                                                                                                   gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa  >  1:2439463/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227451‑4227649

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: