Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,239 A→G 5.9% E67G (GAA→GGA)  PP_0890 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,2390AG 5.9% 201.9 / 2.8 68E67G (GAA→GGA) PP_0890hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base G (2/2);  total (32/36)
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.61e-01

AGGAACGCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCAC  >  NC_002947/1033092‑1033370
                                                                                                                                                   |                                                                                                                                   
gttcctgCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGaa                                                                                                                                    <  4:47388/143‑1 (MQ=255)
aGGAACGCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTGCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGAACTGCCAAGGGTAAGTGAAGCAGTGGCTGCCACTGCGCCTGCTggagga                                                                                                                                     >  1:304496/1‑148 (MQ=255)
 ggAACGCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGaaa                                                                                                                                   <  8:324134/149‑1 (MQ=255)
     cGCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAAcag                                                                                                                               >  5:157977/1‑149 (MQ=255)
       ccAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTTCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGTAGGAAAAcagc                                                                                                                              >  3:251766/1‑148 (MQ=255)
        cAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAAcagc                                                                                                                              >  5:358192/1‑147 (MQ=255)
                  ctgACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGa                                                                                                                  <  1:358406/147‑1 (MQ=255)
                   ggACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAg                                                                                                                 >  8:219372/1‑149 (MQ=255)
                         ccccGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCAc                                                                                                             <  6:145371/147‑1 (MQ=255)
                              gCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGaa                                                                                                      <  7:317618/149‑1 (MQ=255)
                              gCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTCCGACAAGGTGTACTGCCAAGTGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGaa                                                                                                      <  7:26469/149‑1 (MQ=255)
                                    tctcGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCt                                                                                                 >  4:6967/1‑148 (MQ=255)
                                    tctcGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTg                                                                                                >  1:181205/1‑149 (MQ=255)
                                    tctcGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTCCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTg                                                                                                >  6:339974/1‑149 (MQ=255)
                                      tcGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGACACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTggg                                                                                              >  6:97256/1‑149 (MQ=255)
                                        ggCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTggg                                                                                              <  6:259099/147‑1 (MQ=255)
                                        ggCCTACTCCGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTACCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTgggg                                                                                             <  8:119802/148‑1 (MQ=255)
                                           cTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGggcggc                                                                                         <  1:157642/149‑1 (MQ=255)
                                           cTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGggcggc                                                                                         <  5:339974/149‑1 (MQ=255)
                                             aCTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGgcggcg                                                                                        >  2:27272/1‑148 (MQ=255)
                                              cTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGgcggcg                                                                                        >  5:328742/1‑147 (MQ=255)
                                              cTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGgcggcg                                                                                        >  4:359828/1‑147 (MQ=255)
                                               tCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAg                                                                                     >  6:47725/1‑149 (MQ=255)
                                                 ggCCAAGTAATATTGTTCCTGCCACTTTGTCATGGGCTTCGACAATGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTc                                                                                   <  3:108043/149‑1 (MQ=255)
                                                  gCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCa                                                                                  <  7:354864/149‑1 (MQ=255)
                                                    cAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGGAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAgg                                                                                <  7:53057/149‑1 (MQ=255)
                                                    cAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGGAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAgg                                                                                <  7:58846/149‑1 (MQ=255)
                                                    cAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGGAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAgg                                                                                >  8:53057/1‑149 (MQ=255)
                                                    cAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGGAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAgg                                                                                >  8:58846/1‑149 (MQ=255)
                                                     aaGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAggc                                                                               <  6:157977/149‑1 (MQ=255)
                                                      aGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAggcg                                                                              >  3:34428/1‑149 (MQ=255)
                                                         aaTATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGAGGCGGCGCAGTCAGGCGGcc                                                                           <  3:6967/149‑1 (MQ=255)
                                                                ttCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAg                                                                    <  3:359828/149‑1 (MQ=255)
                                                                             gTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTTGCTGCCACTGAGCCTGCTGTATGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAgg                                                        <  6:421947/148‑1 (MQ=255)
                                                                                  gggCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAg                                                         >  4:309231/1‑142 (MQ=255)
                                                                                  gggCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAg                                                         <  3:309231/142‑1 (MQ=255)
                                                                                   ggCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGcc                                                  <  7:219372/148‑1 (MQ=255)
                                                                                   ggCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCg                                                 <  1:27272/149‑1 (MQ=255)
                                                                                    gCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGc                                                <  5:47725/149‑1 (MQ=255)
                                                                                      ttCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCtt                                              >  7:153171/1‑149 (MQ=255)
                                                                                       tCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTtg                                             >  2:398534/1‑149 (MQ=255)
                                                                                        cGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCtt                                              <  6:328742/147‑1 (MQ=255)
                                                                                        cGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTtg                                             <  6:358192/148‑1 (MQ=255)
                                                                                         gACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTtgct                                           >  3:34804/1‑149 (MQ=255)
                                                                                             aGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCa                                       <  1:102987/149‑1 (MQ=255)
                                                                                               gtgtACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCAt                                      >  3:323243/1‑148 (MQ=255)
                                                                                               gtgtACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCAt                                      <  2:181205/148‑1 (MQ=255)
                                                                                                  tACTGCCAAGGTTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGTGGGCTGCCGCTTGCTGCCATGAAt                                  <  4:251766/149‑1 (MQ=255)
                                                                                                  tACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGaa                                   <  8:153171/148‑1 (MQ=255)
                                                                                                  tACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGaa                                   <  4:170323/148‑1 (MQ=255)
                                                                                                     tGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATg                                 >  2:46880/1‑147 (MQ=255)
                                                                                                     tGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATg                                 <  1:46880/147‑1 (MQ=255)
                                                                                                      gCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAg                              >  5:204349/1‑149 (MQ=255)
                                                                                                       ccAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAgg                             >  6:422642/1‑149 (MQ=255)
                                                                                                       caaaGGGTACGTGAAGCAGTGGCTGCAACTGCGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGaa                                                                                                      >  6:381801/3‑76 (MQ=255)
                                                                                                            ggTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAAAGCAGGCggg                         >  7:4117/1‑148 (MQ=255)
                                                                                                                  tGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCAGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTCCCGCTTGCTGCCATCAATGCAGGCGGCTaacg                    >  7:182333/1‑145 (MQ=255)
                                                                                                                        aGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTggcggc            >  6:159352/1‑149 (MQ=255)
                                                                                                                         gTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGcc           >  3:237392/1‑149 (MQ=255)
                                                                                                                          tGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCt          <  2:118206/149‑1 (MQ=255)
                                                                                                                          tGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCt          <  4:34804/149‑1 (MQ=255)
                                                                                                                             cTGCCACTGAGCCTGCTGGAGGAATACAGCGAGCAGCAGCGAGTGAGCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt       <  8:182333/149‑1 (MQ=255)
                                                                                                                             cTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt       <  1:398534/149‑1 (MQ=255)
                                                                                                                                ccACTGATCCTGCTGGAGGAAAACAGCCATCATCAGCGAGTCAACGATGAATGCCTTGGGCGTCTCAGTCATGCTTCATGGAAAGTCGCGAGTGAGGGCTTCAGCTTGATGCCATGAATGCATTCGTGTAATTTTTGCTTCCTGTTTTc    <  2:304496/149‑1 (MQ=255)
                                                                                                                                 cACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGTAGGCCTGGGGCGGCGGAGTCAGGCGGCCTGGCATGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGaa                                   <  3:159201/117‑1 (MQ=255)
                                                                                                                                 cACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTgg        >  3:10386/1‑144 (MQ=255)
                                                                                                                                 cACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTgg        <  4:10386/144‑1 (MQ=255)
                                                                                                                                 cACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTCGGGCCGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGCCTGCCGCTTCCTGCCATGat                                   >  4:159201/1‑116 (MQ=255)
                                                                                                                                  aCTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt       <  7:393736/144‑1 (MQ=255)
                                                                                                                                  aCTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt       >  8:393736/1‑144 (MQ=255)
                                                                                                                                    tGAGCCTCCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCAc  >  5:210537/1‑147 (MQ=255)
                                                                                                                                               gaggaAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGcc                                                                           >  5:14191/1‑63 (MQ=255)
                                                                                                                                               gaggaAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGcc                                                                           <  6:14191/63‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                   
AGGAACGCCAGGCGCTGGCGGACTTCCCCGGCATACTCTCGGCCTACTCGGCCAAGGAATATTGTTCCTGCCGCTTTGTCATGGGCTTCGACAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCAC  >  NC_002947/1033092‑1033370

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: ‑