Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,546 G→T 8.4% T123T (ACC→ACA amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,5460GT 8.4% 613.7 / 33.6 190T123T (ACC→ACAamaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (83/91);  new base T (6/10);  total (89/101)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 196 total aligned reads displayed.

TTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTC  >  NC_002947/6005397‑6005695
                                                                                                                                                     |                                                                                                                                                     
ttCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATg                                                                                                                                                       <  1:4269829/150‑1 (MQ=255)
 tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                                                                                      >  1:1373024/1‑150 (MQ=255)
 tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                                                                                      <  2:1355822/150‑1 (MQ=255)
    aTACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCa                                                                                                                                                   >  1:2701613/1‑150 (MQ=255)
     tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAg                                                                                                                                                  >  1:986141/1‑150 (MQ=255)
        caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                                                                               >  1:512307/1‑150 (MQ=255)
        caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                                                                               >  1:3239115/1‑150 (MQ=255)
          ccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGt                                                                                                                                             <  1:4983923/150‑1 (MQ=255)
          ccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGGTCCAGGCCCACACGGCAACCGGGCAGTTGAAGGCGCTGAACACGCCGACCACACCCAGCGGGTGCCAGCTTTCACGCATGTGGTGGCCCGGGCGCTCGGCAGCGATGGTCAGGCCGt                                                                                                                                             >  2:4955194/1‑150 (MQ=255)
           cGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTa                                                                                                                                            >  1:16739/1‑150 (MQ=255)
                 gCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTg                                                                                                                                      <  1:4249025/150‑1 (MQ=255)
                   cGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgc                                                                                                                                    >  1:4576967/1‑150 (MQ=255)
                     gCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgcgc                                                                                                                                  <  1:4469846/150‑1 (MQ=255)
                     gCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgcgc                                                                                                                                  >  2:622624/1‑150 (MQ=255)
                       caccaGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGa                                                                                                                                >  2:4382082/1‑150 (MQ=255)
                          cagggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAg                                                                                                                             >  1:2293182/1‑150 (MQ=255)
                           agggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAgg                                                                                                                            >  2:2483978/1‑150 (MQ=255)
                            gggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGc                                                                                                                           <  2:2174363/150‑1 (MQ=255)
                               ccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCaa                                                                                                                        <  2:2431702/150‑1 (MQ=255)
                                     cGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGa                                                                                                                  <  2:2098174/150‑1 (MQ=255)
                                      ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                                                 <  1:1555767/150‑1 (MQ=255)
                                      ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                                                 <  1:162056/150‑1 (MQ=255)
                                      ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGCTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                                                 <  2:5388692/150‑1 (MQ=255)
                                          ttCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCg                                                                                                             <  1:2715130/150‑1 (MQ=255)
                                          ttCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCg                                                                                                             >  2:3548848/1‑150 (MQ=255)
                                            ccAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCa                                                                                                           <  1:4846146/150‑1 (MQ=255)
                                            ccAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACACATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCa                                                                                                           <  2:4156800/150‑1 (MQ=255)
                                             cAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAg                                                                                                          >  1:2996862/1‑150 (MQ=255)
                                              aGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGa                                                                                                         <  2:4721978/150‑1 (MQ=255)
                                               ggCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGAt                                                                                                        <  1:669242/150‑1 (MQ=255)
                                               ggCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGCACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGAt                                                                                                        <  1:2338595/150‑1 (MQ=255)
                                                  ccACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTc                                                                                                     >  2:3030600/1‑150 (MQ=255)
                                                           aaCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAgg                                                                                                                                                 <  1:607938/97‑1 (MQ=255)
                                                           aaCAGGGAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTc                                                                                            <  2:3215494/150‑1 (MQ=255)
                                                             ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGACGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCt                                                                                          <  1:2382497/150‑1 (MQ=255)
                                                             ccGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCt                                                                                          <  2:772182/150‑1 (MQ=255)
                                                             cagggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCt                                                                                          <  2:4849328/148‑1 (MQ=255)
                                                              cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCTTGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGTCCGTACAGCTGGCGCGACATGCGTACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                                                         <  1:3367160/150‑1 (MQ=255)
                                                              cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                                                         <  2:2760164/150‑1 (MQ=255)
                                                               gggCAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                                                                                        >  2:3243942/1‑150 (MQ=255)
                                                               gggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGc                                                                                        >  1:4220529/1‑150 (MQ=255)
                                                                ggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                                                 >  1:2754415/1‑124 (MQ=255)
                                                                ggAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                                                 <  2:2754415/124‑1 (MQ=255)
                                                                ggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTgta                                                                                       >  2:4628756/1‑148 (MQ=255)
                                                                     ttGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACAGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCa                                                                                  <  2:2051882/150‑1 (MQ=255)
                                                                       gAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAcc                                                                                <  2:5427751/150‑1 (MQ=255)
                                                                          ggCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATc                                                                                                 >  2:4605290/1‑130 (MQ=255)
                                                                              gagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCt                                                                         <  2:4004611/148‑1 (MQ=255)
                                                                              gagATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCt                                                                         <  2:2002137/148‑1 (MQ=255)
                                                                                gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc                                                                       <  1:5256829/150‑1 (MQ=255)
                                                                                gATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTc                                                                       <  1:4621637/150‑1 (MQ=255)
                                                                                  tCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCt                                                                     >  2:2313092/1‑150 (MQ=255)
                                                                                    aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                   <  1:528525/150‑1 (MQ=255)
                                                                                    aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                   >  1:4866343/1‑150 (MQ=255)
                                                                                    aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                   <  2:4869976/150‑1 (MQ=255)
                                                                                    aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                   >  1:4722022/1‑150 (MQ=255)
                                                                                      gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                                                 <  2:5267297/150‑1 (MQ=255)
                                                                                       ccGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTg                                                                <  1:3866624/150‑1 (MQ=255)
                                                                                         gACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGAt                                                              <  2:4299628/150‑1 (MQ=255)
                                                                                         gACCACACCCAGCGGGTGCCAGGTTTCACGCAAGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                                              <  2:125719/150‑2 (MQ=255)
                                                                                            cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTAGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                                                           >  2:5219945/1‑150 (MQ=255)
                                                                                             acacCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGaagtgc                                                          <  2:5315282/150‑6 (MQ=255)
                                                                                              cacCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGc                                                         <  1:2482385/150‑1 (MQ=255)
                                                                                                cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                       <  2:5073871/150‑1 (MQ=255)
                                                                                                cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                       <  2:701939/150‑1 (MQ=255)
                                                                                                  cAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGa                                                                                                   >  1:2459317/1‑104 (MQ=255)
                                                                                                       ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                <  2:2954615/150‑1 (MQ=255)
                                                                                                        atGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcgag                                               <  2:4060326/149‑2 (MQ=255)
                                                                                                         tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                                              <  2:964454/150‑1 (MQ=255)
                                                                                                           ccAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcgagcga                                            >  1:2396596/1‑150 (MQ=255)
                                                                                                           ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                            >  1:192613/1‑150 (MQ=255)
                                                                                                           ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                            >  1:3770085/1‑150 (MQ=255)
                                                                                                            cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                            >  1:3160155/1‑149 (MQ=255)
                                                                                                            cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                           >  2:186214/1‑150 (MQ=255)
                                                                                                            cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                           >  2:4060600/1‑150 (MQ=255)
                                                                                                                tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa                                       >  2:1766675/1‑150 (MQ=255)
                                                                                                                          gtggtgGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAg                             >  1:3926149/1‑150 (MQ=255)
                                                                                                                            ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTAGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                               <  2:4098194/130‑1 (MQ=255)
                                                                                                                            ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTAGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                               >  1:4098194/1‑130 (MQ=255)
                                                                                                                             gtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                          <  2:4220529/150‑1 (MQ=255)
                                                                                                                              tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGa                                                               >  1:3987897/1‑112 (MQ=255)
                                                                                                                              tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                         >  1:1043364/1‑150 (MQ=255)
                                                                                                                                gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGCGATCTTGCCGGCTtcgatc                                              <  2:5203614/127‑1 (MQ=255)
                                                                                                                                  ccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCt                     >  1:2941850/1‑150 (MQ=255)
                                                                                                                                   cGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                    <  2:16739/150‑1 (MQ=255)
                                                                                                                                    gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtg                   <  2:3851136/150‑1 (MQ=255)
                                                                                                                                    gcgcgtTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTtg                   <  1:1948738/148‑1 (MQ=255)
                                                                                                                                     ggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                  <  1:4915525/150‑1 (MQ=255)
                                                                                                                                     ggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                  <  2:1430125/150‑1 (MQ=255)
                                                                                                                                          tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCgg                        <  2:2632038/139‑1 (MQ=255)
                                                                                                                                          tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa             >  1:3446782/1‑150 (MQ=255)
                                                                                                                                          tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa             >  1:2553464/1‑150 (MQ=255)
                                                                                                                                          tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa             >  1:2240998/1‑150 (MQ=255)
                                                                                                                                            ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg           >  2:4064524/1‑150 (MQ=255)
                                                                                                                                            ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg           >  2:1293873/1‑150 (MQ=255)
                                                                                                                                             gAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgc          >  1:5359240/1‑150 (MQ=255)
                                                                                                                                              aaGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgca         >  1:4099313/1‑150 (MQ=255)
                                                                                                                                                  gATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCCCTTCACCCAGGCCGTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCGCGCCGAGGTCGGCCGGGTGCTCACGCAGCAc     >  2:5442148/1‑150 (MQ=255)
                                                                                                                                                  gATGGTCAGGCCGTAAAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCAc     >  2:2211669/1‑150 (MQ=255)
                                                                                                                                                     ggTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTc  <  1:3290011/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTC  >  NC_002947/6005397‑6005695

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