Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,604,251 C→G 7.8% P12R (CCG→CGG) ‡ PP_4076 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,604,2510CG 7.8% 269.9 / 12.3 90P12R (CCG→CGG) ‡PP_4076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/44);  new base G (4/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4604104‑4604399
                                                                                                                                                   |                                                                                                                                                    
tcGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGGACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGt                                                                                                                                                    >  2:1594481/1‑150 (MQ=255)
           gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATcaag                                                                                                                                         <  1:2052095/150‑4 (MQ=255)
             gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       <  2:1446184/150‑1 (MQ=255)
             gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       <  1:1490964/150‑1 (MQ=255)
             gTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                       <  2:591779/150‑1 (MQ=255)
                 tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCg                                                                                                                                   >  2:350517/1‑150 (MQ=255)
                    gggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTTGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtct                                                                                                                                >  1:69961/1‑150 (MQ=255)
                    gggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtct                                                                                                                                >  1:2320953/1‑150 (MQ=255)
                     ggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtctc                                                                                                                               <  1:1335009/150‑1 (MQ=255)
                          gtgccAGTTCTTGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATCTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGCGGCGTCTCGCCAg                                                                                                                          <  2:1181932/147‑1 (MQ=25)
                           tACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGGCTCGCCAgg                                                                                                                         >  1:829780/1‑150 (MQ=255)
                              cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGACCTTCCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      <  1:626919/150‑1 (MQ=255)
                              cAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTTGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      <  2:69961/150‑1 (MQ=255)
                                     cgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGcc                                                                                                                                               <  2:2003772/118‑1 (MQ=255)
                                     cgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGcc                                                                                                                                               >  1:2003772/1‑118 (MQ=255)
                                     cgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                               <  2:245573/150‑1 (MQ=255)
                                      gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGACAGGAAATCTCTTCg                                                                                                              >  2:1752672/1‑150 (MQ=255)
                                        tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGca                                                                                                            >  2:1124411/1‑150 (MQ=255)
                                          gcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGcaca                                                                                                          >  1:208085/1‑150 (MQ=255)
                                          gcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGcaca                                                                                                          >  1:1315991/1‑150 (MQ=255)
                                           ccGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCCGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGGCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                         <  1:536880/150‑1 (MQ=255)
                                               gcgcCCGCACTTCCCGTACATGGGCGCGAACTTCCACATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTc                                                                                                     <  1:1075492/150‑1 (MQ=255)
                                                 gcCCGCACTTCGAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgc                                                                                                   <  2:2155444/150‑1 (MQ=255)
                                                  cccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                  <  2:1826831/150‑1 (MQ=255)
                                                  cccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                  <  1:767097/150‑1 (MQ=255)
                                                         ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAg                                                                                                                     >  1:2042066/1‑124 (MQ=255)
                                                         ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAg                                                                                                                     <  2:2042066/124‑1 (MQ=255)
                                                              gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                          <  2:1528449/114‑1 (MQ=255)
                                                              gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                          >  1:1528449/1‑114 (MQ=255)
                                                              gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        <  1:2285843/148‑1 (MQ=255)
                                                              gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                        >  2:2285843/1‑148 (MQ=255)
                                                                  aTGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGCCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAcc                                                                                  <  1:648104/150‑1 (MQ=255)
                                                                  aTGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCATCTCGGCAGGAAAGTGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAcc                                                                                  <  1:647736/150‑1 (MQ=255)
                                                                   tGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCGTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGGGCCACCt                                                                                 >  2:1239597/1‑150 (MQ=255)
                                                                        gcgAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaaa                                                                            <  2:1420410/150‑1 (MQ=255)
                                                                        gcgAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGCTTCACTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaaa                                                                            >  1:2254049/1‑150 (MQ=255)
                                                                            aCTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGtt                                                                        <  1:1124411/150‑1 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  2:362312/150‑1 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  1:1594481/150‑1 (MQ=255)
                                                                              ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  2:197455/1‑150 (MQ=255)
                                                                               tCGAAGCTTCCCTGGTAGGCCGCCGCGACACGCCATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGCCCGTCATCCGGGAGGCCGCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCt                                                                     <  2:2079868/150‑1 (MQ=255)
                                                                                 gAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                          >  1:1982149/1‑95 (MQ=255)
                                                                                 gAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                          <  2:1982149/95‑1 (MQ=255)
                                                                                  aaaTTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGACGTCGTACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGaa                                                                  <  1:1845269/150‑1 (MQ=255)
                                                                                      ttCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCt                                                              >  1:334981/1‑150 (MQ=255)
                                                                                          cTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      >  2:2418795/1‑90 (MQ=255)
                                                                                          cTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGaaa                                                                                                                      <  1:2418795/90‑1 (MQ=255)
                                                                                               aGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCAAGGAACGATACTAGCCGTGCCTCATTAGCAGCAGGTGCTACCTCGAAACGTTGCTGAACGCTCGCCAAGGt                                                     <  1:46176/150‑1 (MQ=255)
                                                                                                  ccGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACCAGTTTGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGc                                                  <  1:1239597/150‑1 (MQ=255)
                                                                                                   cGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                                <  1:1433019/119‑1 (MQ=255)
                                                                                                   cGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTc                                                                                >  2:1433019/1‑119 (MQ=255)
                                                                                                    ggCGCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTc                                                <  2:494264/150‑1 (MQ=255)
                                                                                                         gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCaga                                           <  2:2254049/150‑1 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                    >  1:577514/1‑58 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                    <  2:577514/58‑1 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                        <  2:1284543/150‑1 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                        >  2:554199/1‑150 (MQ=255)
                                                                                                            acGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagca                                        >  1:1405674/1‑150 (MQ=255)
                                                                                                            acGCAATGAGTGGTGTCTTCGATCGCTTGCTGGCGGACCCGTCCTCCATCCGGGAGGCGTCTCCCCTGGCAAGCGCGTCGCCCAAGTTCGGGGGCATCAAGTGCCACCTCGAAACGCTGCTGCACGCTCCCCCAGGTTCCCCCCAGagca                                        >  2:760314/1‑150 (MQ=25)
                                                                                                              gCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                      >  1:1576618/1‑102 (MQ=255)
                                                                                                              gCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGc                                                                                      <  2:1576618/102‑1 (MQ=255)
                                                                                                              gCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGagcagc                                      >  2:442804/1‑150 (MQ=255)
                                                                                                                 aTGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGcccc                                   <  2:208085/150‑1 (MQ=255)
                                                                                                                    aGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAg                                >  1:1781993/1‑150 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTcc                                               <  2:155687/131‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            <  1:350517/150‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            <  2:829780/150‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTgg                            >  2:2166271/1‑150 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGCCCATCAAGTGCCACCTAACATCCCAGCTGAACGCTCGCCAAGGTTGCTCCCGGAGCAGCCCCGAGCTgg                            <  1:416453/150‑1 (MQ=255)
                                                                                                                        gTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTAGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTcc                                               >  1:155687/1‑131 (MQ=255)
                                                                                                                         tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTggg                           >  2:538290/1‑150 (MQ=255)
                                                                                                                                 aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGcgcg                   >  2:834620/1‑150 (MQ=255)
                                                                                                                                  tCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGa                  <  2:2320953/150‑1 (MQ=255)
                                                                                                                                   tcccggATGGCGGACGGGTCCGCCATCCGGTAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGAc                 >  1:1025421/6‑150 (MQ=25)
                                                                                                                                   tcccggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGAc                 >  1:2280365/6‑150 (MQ=255)
                                                                                                                                    gCTTGATGGCGGACCCGTCCGCCATACGGGAGGCGTCTCGCCATTAAAGCGCTGCGCACACGTTCGCGGCCATCCCCTGTTAGGGCAACCTCGCGTTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACt                >  1:1125224/1‑150 (MQ=19)
                                                                                                                                     cTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACtt               >  1:650431/1‑150 (MQ=255)
                                                                                                                                     cTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACtt               >  2:1581270/1‑150 (MQ=255)
                                                                                                                                      ttGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTGGACAGGCAAGCGCTTAGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGACCCGAGCTGGGCCTGCGCGACTTc              >  1:892552/1‑150 (MQ=25)
                                                                                                                                      ttGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc              >  1:2294542/1‑150 (MQ=255)
                                                                                                                                        gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  2:1890591/92‑1 (MQ=255)
                                                                                                                                        gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  1:1890591/1‑92 (MQ=255)
                                                                                                                                        gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa            >  2:1874748/1‑150 (MQ=255)
                                                                                                                                         aTGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAAc           <  1:1581270/150‑1 (MQ=255)
                                                                                                                                          tGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                               >  1:211712/1‑49 (MQ=37)
                                                                                                                                          tGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                               <  2:211712/49‑1 (MQ=37)
                                                                                                                                           ggCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACgg         <  1:554199/150‑1 (MQ=255)
                                                                                                                                            gCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACggg        <  1:538290/150‑1 (MQ=255)
                                                                                                                                            gCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACggg        >  2:1939432/1‑150 (MQ=255)
                                                                                                                                                aCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATg    >  1:1120762/1‑150 (MQ=255)
                                                                                                                                                 cccGTCCGGCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGa   <  2:1315991/150‑1 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      <  2:693595/82‑1 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                      >  1:693595/1‑82 (MQ=255)
                                                                                                                                                  ccGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAc  >  1:617811/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGGACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCATGAC  >  NC_002947/4604104‑4604399

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