Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,669,964 C→G 6.8% H147Q (CAC→CAG metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,669,9640CG 6.8% 232.4 / 5.6 73H147Q (CAC→CAGmetF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base G (3/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

TGCTGGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669816‑5670111
                                                                                                                                                       |                                                                                                                                                    
tgctgGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAc                                                                                                                                                      >  5:117583/1‑149 (MQ=255)
     gCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTgg                                                                                                                                                 <  8:126069/148‑1 (MQ=255)
       cGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAg                                                                                                                                              >  1:105136/1‑149 (MQ=255)
            aCAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTCGCt                                                                                                                                         >  1:163608/1‑149 (MQ=255)
             cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCt                                                                                                                                         <  1:62430/148‑1 (MQ=255)
             cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTg                                                                                                                                        <  8:145726/149‑1 (MQ=255)
              acggCAGCCGGCATCAAACGCTTGGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGGCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatcag                                                                                                                                       <  8:299341/147‑6 (MQ=255)
                   aGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAGCTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTAcc                                                                                                                                  >  5:57909/1‑149 (MQ=255)
                   aGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTAcc                                                                                                                                  >  1:36832/1‑149 (MQ=255)
                    gCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAGCTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTAccc                                                                                                                                 <  6:57909/149‑1 (MQ=255)
                          aTCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAAt                                                                                                                           <  4:143735/149‑1 (MQ=255)
                            cAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGaaa                                                                                                                            >  7:177045/1‑146 (MQ=255)
                            cAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGaaa                                                                                                                            <  8:177045/146‑1 (MQ=255)
                             aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                        <  1:150941/149‑1 (MQ=255)
                             aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                        >  2:284722/1‑149 (MQ=255)
                             aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                        >  2:152670/1‑149 (MQ=255)
                               aCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAcc                                                                                                                      >  3:107913/1‑149 (MQ=255)
                                         gcgcTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGATGGGCGAGTTCATCCGCGAGGAGCCCGCTGATCAATTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGgc                                                                                                               >  1:112312/1‑149 (MQ=255)
                                          cgcTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATg                                                                                                                          >  2:232069/1‑134 (MQ=255)
                                          cgcTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATg                                                                                                                          <  1:232069/134‑1 (MQ=255)
                                               cGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACtt                                                                                                      <  2:105136/149‑1 (MQ=255)
                                                    tgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAgg                                                                                                 >  8:30831/1‑149 (MQ=255)
                                                      aCCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGc                                                                                                <  6:316852/148‑1 (MQ=255)
                                                      aCCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCg                                                                                               >  2:22563/1‑149 (MQ=255)
                                                        cTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGa                                                                                             >  4:189387/1‑149 (MQ=255)
                                                        cTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACACGCAGGCGCTCAACTTCGAGGCCGa                                                                                             >  6:270650/1‑149 (MQ=255)
                                                          gCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATc                                                                                           >  4:294571/1‑149 (MQ=255)
                                                            cGTCGGGCATGGGCATGGGCAGTGGTTAACTGCGCTATGCCAGTGAGC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTg                                                                                         <  2:163608/149‑1 (MQ=255)
                                                              tCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGc                                                                                       <  1:152670/149‑1 (MQ=255)
                                                                ggggctgggcatgggcaGTGGTTAACTGCGCTATGCCATTGACC‑‑‑TGGGCGTTTTAAACCGCGATGTGAACGGTGCTCCCATTCACTTTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGTATCTGGc                                                                                       <  1:284722/143‑1 (MQ=255)
                                                                  gcatgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAAc                                                                                   >  6:350479/1‑149 (MQ=255)
                                                                    atgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACtt                                                                                 >  2:231540/1‑149 (MQ=255)
                                                                     tgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACttt                                                                                <  5:350479/149‑1 (MQ=255)
                                                                      gggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtg                                                                               >  4:113876/1‑149 (MQ=255)
                                                                       ggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtgt                                                                              <  2:112312/149‑1 (MQ=255)
                                                                       ggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtg                                                                               <  3:189387/148‑1 (MQ=255)
                                                                        gcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCGGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtgtg                                                                             <  7:117831/149‑1 (MQ=255)
                                                                        gcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACttt                                                                                >  4:218268/1‑146 (MQ=255)
                                                                        gcatggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACttt                                                                                <  3:218268/146‑1 (MQ=255)
                                                                            ggCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGcac                                                                          >  7:82740/1‑148 (MQ=255)
                                                                             gCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACaa                                                                        <  7:30831/149‑1 (MQ=255)
                                                                                    gtgAACTGCGCTATGGCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGGCTACCCTGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCaa                                                                  <  5:270650/148‑1 (MQ=255)
                                                                                        aCTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCg                                                             <  1:231540/149‑1 (MQ=255)
                                                                                         cTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCgg                                                            <  2:36832/149‑1 (MQ=255)
                                                                                            cgcTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGc                                                         <  3:294571/149‑1 (MQ=255)
                                                                                                   ccAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAgc                                                  >  6:200928/1‑149 (MQ=255)
                                                                                                           c‑‑‑TGGTCGAGTTCATCCGCGAGGCGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTCACGGGAAATGCACCCGCAGGCACGCGACTTGGAGGCGGATCTGGGCAACTTTGTGCGCAAGGTCAGGGCCGGTGCGGCCAGCGCCATctc                                          >  7:235182/1‑147 (MQ=255)
                                                                                                           c‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAgcgc                                                <  4:339510/143‑1 (MQ=255)
                                                                                                           c‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAgcgc                                                >  3:339510/1‑143 (MQ=255)
                                                                                                           c‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAc                                          >  7:235178/1‑149 (MQ=255)
                                                                                                               tGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGt                                                           <  5:68450/131‑1 (MQ=255)
                                                                                                               tGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGt                                                           >  6:68450/1‑131 (MQ=255)
                                                                                                                ggTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAg                                                                 <  3:248131/124‑1 (MQ=255)
                                                                                                                ggTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGcc                                                        >  8:49644/1‑133 (MQ=255)
                                                                                                                ggTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGcc                                                        <  7:49644/133‑1 (MQ=255)
                                                                                                                ggTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCCCAACATCGAGGCGGATCAGGCCAACTTTGTGCACAAAGTCAAg                                                                 >  4:248131/1‑124 (MQ=255)
                                                                                                                 gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         <  8:82740/147‑1 (MQ=255)
                                                                                                                 gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         >  8:79505/1‑147 (MQ=255)
                                                                                                                 gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCa                                       <  4:107913/149‑1 (MQ=255)
                                                                                                                 gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGACTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCa                                       <  6:117583/149‑1 (MQ=255)
                                                                                                                          aTCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActtct                               <  5:200928/148‑1 (MQ=255)
                                                                                                                           tCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActtctt                              >  4:83580/1‑148 (MQ=255)
                                                                                                                              gcgAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAAc                          >  7:224490/1‑149 (MQ=255)
                                                                                                                                        cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  2:236163/1‑148 (MQ=255)
                                                                                                                                        cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  1:220731/1‑148 (MQ=255)
                                                                                                                                        cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGAGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  5:334159/1‑148 (MQ=255)
                                                                                                                                         gCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGACATAACCCAGTACTTCTTCAACGCCGACAGcta               >  7:169020/1‑149 (MQ=255)
                                                                                                                                           cagccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctact             >  7:200533/5‑149 (MQ=255)
                                                                                                                                            agccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACttcttc                             <  2:335324/129‑1 (MQ=255)
                                                                                                                                            agccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACttcttc                             >  1:335324/4‑132 (MQ=255)
                                                                                                                                              ccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc           >  5:323659/2‑148 (MQ=255)
                                                                                                                                               cACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttct          >  4:341094/1‑148 (MQ=255)
                                                                                                                                                  ttCCACCTGGAAGTGGCTTCCTACACGGAAATGCACCCGAAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCAAGTACTTCTTCAACGACGAAAGctacttctact       >  5:259851/1‑148 (MQ=255)
                                                                                                                                                      aCCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcgtc  <  3:113876/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TGCTGGCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669816‑5670111

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