Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,669,965 C→G 7.3% L148V (CTG→GTG)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,669,9650CG 7.3% 220.0 / 4.3 73L148V (CTG→GTG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (36/31);  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.43e-01

GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTCGA  >  NC_002947/5669821‑5670113
                                                                                                                                                    |                                                                                                                                                    
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)
                                                                                                                                                    cTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcgtc    >  6:194046/1‑147 (MQ=255)
                                                                                                                                                    cTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTCGa  >  8:324867/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACC‑‑‑TGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTCGA  >  NC_002947/5669821‑5670113

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: