Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,669,964 C→G 8.2% 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 8.2% 235.7 / 9.1 73H147Q (CAC→CAG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base G (2/4);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACT  >  NC_002947/5669821‑5670106
                                                                                                                                               |                                                                                                                                              
gCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGAGCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTg                                                                                                                                             <  8:215012/147‑1 (MQ=255)
     gTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGCGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTgg                                                                                                                                      <  4:64344/149‑1 (MQ=255)
      tACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGc                                                                                                                                     >  2:335915/1‑149 (MQ=255)
        cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                   <  2:49328/149‑5 (MQ=255)
        cAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatca                                                                                                                                   <  1:214942/149‑5 (MQ=255)
         aaGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGc                                                                                                                                  <  3:227427/149‑1 (MQ=255)
           ggCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatcagcg                                                                                                                                <  5:323022/149‑8 (MQ=255)
                ccGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTgg                                                                                                                                      <  8:130712/138‑1 (MQ=255)
                ccGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTgg                                                                                                                                      >  7:130712/1‑138 (MQ=255)
                   gCATCAAACGCATCGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAcc                                                                                                                                               <  5:233901/126‑1 (MQ=255)
                   gCATCAAACGCATCGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAcc                                                                                                                                               >  6:233901/1‑126 (MQ=255)
                    cATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGaaa                                                                                                                       <  5:28871/149‑1 (MQ=255)
                     aTCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAAt                                                                                                                      <  2:68795/149‑1 (MQ=255)
                     aTCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAAt                                                                                                                      <  8:345839/149‑1 (MQ=255)
                      tCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATg                                                                                                                     >  1:118269/1‑149 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatcag                                                                                                                                  <  6:89511/134‑6 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCAGGTGGAAGTgatcag                                                                                                                                  >  5:89511/1‑129 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAg                                                                                                                                         <  6:342433/127‑1 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAg                                                                                                                                         >  5:342433/1‑127 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                   >  3:335219/1‑149 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                   >  3:252323/1‑149 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                   >  6:332/1‑149 (MQ=255)
                              aTTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCg                                                                                                               >  7:172289/1‑147 (MQ=255)
                               ttGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAg                                                                                                            <  4:335219/149‑1 (MQ=255)
                                tGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAg                                                                                                            <  4:252323/148‑1 (MQ=255)
                                    gcgcTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGgcgcg                                                                                                       >  4:95381/1‑149 (MQ=255)
                                              gtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAg                                                                                             >  7:355532/1‑149 (MQ=255)
                                                 aCCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAc                                                                                                                  >  1:364361/1‑125 (MQ=255)
                                                 aCCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAc                                                                                                                  <  2:364361/125‑1 (MQ=255)
                                                    tGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGAt                                                                                       >  1:330384/1‑149 (MQ=255)
                                                    tGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGAt                                                                                       >  5:255457/1‑149 (MQ=255)
                                                        gTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTgg                                                                                   <  2:113644/149‑1 (MQ=255)
                                                         tCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACAGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGTATCTGGc                                                                                  >  3:275961/1‑149 (MQ=255)
                                                          cGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGcc                                                                                 <  5:64266/149‑1 (MQ=255)
                                                             gcatgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAAc                                                                              >  3:39805/1‑149 (MQ=255)
                                                              catgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACt                                                                             <  1:335915/149‑1 (MQ=255)
                                                              catgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACt                                                                             >  2:50088/1‑149 (MQ=255)
                                                               atgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCa                                                                                                     <  8:115539/124‑1 (MQ=255)
                                                               atgggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCa                                                                                                     >  7:115539/1‑124 (MQ=255)
                                                                            gtggtgAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTc                                                               <  2:125809/149‑1 (MQ=255)
                                                                               gtgAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCaa                                                             <  2:74618/148‑1 (MQ=255)
                                                                               gtgAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAg                                                            >  1:79343/1‑149 (MQ=255)
                                                                               gtgAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAg                                                            <  1:50088/149‑1 (MQ=255)
                                                                                  aaCTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGc                                                          <  6:255457/148‑1 (MQ=255)
                                                                                            tGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGa                                                 <  8:355532/147‑1 (MQ=255)
                                                                                               cAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAgcg                                            >  3:123925/1‑149 (MQ=255)
                                                                                                 gTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGcc                                          <  2:118269/149‑1 (MQ=255)
                                                                                                      cTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAc                                     >  7:296039/1‑149 (MQ=255)
                                                                                                      cTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAc                                     <  8:172289/149‑1 (MQ=255)
                                                                                                      cTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAc                                     >  1:229563/1‑149 (MQ=255)
                                                                                                         gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGTATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                    <  4:275961/147‑1 (MQ=255)
                                                                                                         gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCa                                         <  6:105394/142‑1 (MQ=255)
                                                                                                         gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCa                                         >  5:105394/1‑142 (MQ=255)
                                                                                                         gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                    <  2:79343/147‑1 (MQ=255)
                                                                                                         gTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCa                                  <  2:350789/149‑1 (MQ=255)
                                                                                                            gAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                              >  7:82756/1‑86 (MQ=255)
                                                                                                            gAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                              <  8:82756/86‑1 (MQ=255)
                                                                                                            gAGTACATCCGCCACCAGACCGCTGATAACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCACGCGCCCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATAACCCAGt                                >  8:138954/1‑148 (MQ=255)
                                                                                                               ttCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActt                            <  3:117800/149‑1 (MQ=255)
                                                                                                                   tCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACttcttc                        <  4:39805/149‑1 (MQ=255)
                                                                                                                             gaCCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGa                >  1:338816/1‑147 (MQ=255)
                                                                                                                             gaCCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACa              <  2:330384/149‑1 (MQ=255)
                                                                                                                                cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                              >  7:296108/1‑66 (MQ=255)
                                                                                                                                cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                              <  8:296108/66‑1 (MQ=255)
                                                                                                                                cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc            >  1:52235/1‑148 (MQ=255)
                                                                                                                                cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc            >  8:354445/1‑148 (MQ=255)
                                                                                                                                 gCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGcta          >  6:249176/1‑149 (MQ=255)
                                                                                                                                  cTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctac         <  3:95381/149‑1 (MQ=255)
                                                                                                                                   tGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctact        >  1:188980/1‑149 (MQ=255)
                                                                                                                                     aTCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc      >  5:378803/1‑149 (MQ=255)
                                                                                                                                     aTCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc      <  8:296039/149‑1 (MQ=255)
                                                                                                                                     aTCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctactt       <  2:229563/148‑1 (MQ=255)
                                                                                                                                       cACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGCCAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttcta    >  3:153143/1‑149 (MQ=255)
                                                                                                                                       cACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttct     >  7:175045/1‑148 (MQ=255)
                                                                                                                                         cTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttctact  >  1:151234/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACT  >  NC_002947/5669821‑5670106

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