Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,337,829 G→C 8.7% S330S (TCG→TCC) ‡ PP_2934 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,337,8290GC 8.7% 396.7 / 24.1 127S330S (TCG→TCC) ‡PP_2934hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (62/54);  new base C (6/5);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 133 total aligned reads displayed.

GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAA  >  NC_002947/3337680‑3337978
                                                                                                                                                      |                                                                                                                                                      
gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCg                                                                                                                                                        >  1:459767/1‑150 (MQ=255)
gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCg                                                                                                                                                        >  2:1090836/1‑150 (MQ=255)
      ttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgcg                                                                                                                                                  <  2:291682/150‑1 (MQ=255)
         gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCGGCGCTGGCGACCCCAGCGcctg                                                                                                                                              <  2:971690/150‑3 (MQ=255)
          tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCgg                                                                                                                                             >  2:1080938/1‑150 (MQ=255)
           aTTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGc                                                                                                                                            <  2:385714/150‑1 (MQ=255)
                  gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGt                                                                                                                                     >  1:2897353/1‑150 (MQ=255)
                    ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTc                                                                                                                                   >  2:2723758/1‑150 (MQ=255)
                    ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCGGCGCTGGCGACCCCAGCGCCTGGCTGaa                                                                                                                                        <  2:2192739/145‑2 (MQ=38)
                     cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                                                                                                  <  1:510090/150‑1 (MQ=255)
                     cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                                                                                                  >  2:2394618/1‑150 (MQ=255)
                     cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                                                                                                  >  2:2579541/1‑150 (MQ=255)
                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAg                                                                                                                                 >  1:1475129/1‑150 (MQ=255)
                          cGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGc                                                                                                                             <  2:337927/150‑1 (MQ=255)
                            aTGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGa                                                                                                                           <  2:2529574/150‑1 (MQ=255)
                             tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGAc                                                                                                                          <  1:632378/150‑1 (MQ=255)
                             tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCCAc                                                                                                                          >  2:226635/1‑150 (MQ=255)
                              gCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGcc                                                                                                                                                      >  1:1064035/1‑122 (MQ=255)
                               ccGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACgg                                                                                                                        <  1:1080938/150‑1 (MQ=255)
                                       tACGCGCGAATACCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACggg                                                                                                                       <  1:2095640/143‑1 (MQ=255)
                                       tACGCGCGAATACCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACggg                                                                                                                       >  2:2095640/1‑143 (MQ=255)
                                          gcgcgAATGCCGGCGGTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGcgcgccgcga                                                                                                             <  2:493716/150‑2 (MQ=255)
                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTc                                                                                                                                   <  1:255883/128‑1 (MQ=255)
                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                             >  2:253321/1‑150 (MQ=255)
                                            gcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACgg                                                                                                                        <  2:1273751/137‑1 (MQ=255)
                                            gcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCAAGCGCCGGCTGAGCGTTCAGCTGCGACgg                                                                                                                        >  1:1273751/1‑137 (MQ=255)
                                                 tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCgcgccgcgcccaggaa                                                                                                      <  2:976978/150‑7 (MQ=255)
                                                        gcTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTa                                                                                               <  1:2399442/150‑1 (MQ=255)
                                                          tgtgCAAGGGCGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCGGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                             <  1:2736443/150‑1 (MQ=255)
                                                          tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGc                                                                                                                                       >  1:1324179/1‑108 (MQ=255)
                                                          tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGc                                                                                                                                       <  2:1324179/108‑1 (MQ=255)
                                                          tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                             >  1:394400/1‑150 (MQ=255)
                                                          tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                             <  2:2787784/150‑1 (MQ=255)
                                                           gtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGc                                                                                            >  2:955778/1‑150 (MQ=255)
                                                            tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc                                                                                           >  2:462858/1‑150 (MQ=255)
                                                               aaGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGGGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAc                                                                                        >  1:650160/1‑150 (MQ=255)
                                                                 ggACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAccc                                                                                      >  1:272960/1‑150 (MQ=255)
                                                                       cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                                <  2:2044133/150‑1 (MQ=255)
                                                                         aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                              <  1:378357/150‑1 (MQ=255)
                                                                           ggCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                                    <  1:1838042/126‑1 (MQ=255)
                                                                           ggCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                                    >  2:1838042/1‑126 (MQ=255)
                                                                             cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                          >  2:1009250/1‑118 (MQ=255)
                                                                             cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                          <  1:1009250/118‑1 (MQ=255)
                                                                             cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGa                                                                          >  2:1061894/1‑150 (MQ=255)
                                                                              gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCGACAAAGTTGGGCACGCCCCGGGCAAGCCGAg                                                                         >  2:1426355/1‑150 (MQ=255)
                                                                              gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                                         >  1:888511/1‑150 (MQ=255)
                                                                                  aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                     >  1:1053781/1‑150 (MQ=255)
                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAg                                                                 >  1:2180285/1‑150 (MQ=255)
                                                                                         gTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCa                                                              <  1:2061029/150‑1 (MQ=255)
                                                                                               ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagcca                                                                                                       >  2:270331/1‑103 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCcacaca                                                          <  2:1487058/147‑1 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCAcaca                                                          >  1:1487058/1‑144 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                       >  1:787461/1‑150 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTCAAGCAGAAGGCGGCTTGGGTCGCCATCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                       <  1:687015/150‑1 (MQ=255)
                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                         >  1:2411732/1‑99 (MQ=255)
                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                         <  2:2411732/99‑1 (MQ=255)
                                                                                                   catTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGGCGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGAt                                                    >  2:2400186/1‑150 (MQ=255)
                                                                                                         cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                              >  1:1184309/1‑150 (MQ=255)
                                                                                                         cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                              <  2:787461/150‑1 (MQ=255)
                                                                                                         cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                              <  1:955778/150‑1 (MQ=255)
                                                                                                           gCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCAt                                            >  1:2468304/1‑150 (MQ=255)
                                                                                                            cAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCTCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATa                                           >  1:2882271/1‑150 (MQ=255)
                                                                                                             aGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAc                                          <  1:1188138/150‑1 (MQ=255)
                                                                                                               cGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGTTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCacagga                                        <  1:2410952/150‑6 (MQ=255)
                                                                                                               cGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAccc                                        >  2:2530426/1‑150 (MQ=255)
                                                                                                                tagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                       >  1:893182/2‑150 (MQ=255)
                                                                                                                gagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGTCGTGCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                       >  1:968941/1‑150 (MQ=255)
                                                                                                                gagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                       >  2:2182687/1‑150 (MQ=255)
                                                                                                                 agaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGACACGCAGAAGATCGAGGCATACCCCa                                      >  1:1946365/1‑150 (MQ=255)
                                                                                                                   aaaCGCCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAgg                                    <  1:515022/150‑1 (MQ=255)
                                                                                                                    aaCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGgt                                   <  1:253321/150‑1 (MQ=255)
                                                                                                                      cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCa                                                                  <  1:1746735/117‑1 (MQ=255)
                                                                                                                       ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCTCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTa                                <  2:2882271/150‑1 (MQ=255)
                                                                                                                         cAGCGTTCTTCAGCCAGGCGCTGGGGTCGTCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCGGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                              <  2:1013331/150‑1 (MQ=255)
                                                                                                                          aGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                     <  2:735127/110‑1 (MQ=255)
                                                                                                                             gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                                 >  1:224612/1‑127 (MQ=255)
                                                                                                                             gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                                 <  2:224612/127‑1 (MQ=255)
                                                                                                                              ttcttcAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTTAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCAc                         >  1:2845781/1‑150 (MQ=255)
                                                                                                                              ttcttcAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCAc                         >  1:1896164/1‑150 (MQ=255)
                                                                                                                               cgctcAGCCGGCGCTGGCGACCCCAGCGCCTGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAc                               >  2:1651681/4‑142 (MQ=255)
                                                                                                                               cgctcAGCCGGCGCTGGCGACCCCAGCGCCTGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAc                               <  1:1651681/139‑1 (MQ=255)
                                                                                                                                gctcAGCCGGCGCTGGCGACCCCAGCGCCTGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAgg                                                                             <  1:3041133/93‑1 (MQ=37)
                                                                                                                                gctcAGCCGGCGCTGGCGACCCCAGCGCCTGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAgg                                                                             >  2:3041133/3‑95 (MQ=25)
                                                                                                                                 ttcAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgc                      >  2:1223027/1‑150 (MQ=255)
                                                                                                                                  tcAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgcg                     >  2:2417644/1‑150 (MQ=255)
                                                                                                                                  tcAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgcg                     <  1:2530426/150‑1 (MQ=255)
                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCCGGGTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGAACGAGGCATACCCCAGCTGTACCGCCACGGCGccc                  <  2:650160/150‑1 (MQ=255)
                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGccc                  >  2:2469432/1‑150 (MQ=255)
                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGccc                  <  2:1896164/150‑1 (MQ=255)
                                                                                                                                      gcAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGACCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa                 >  2:1825377/2‑150 (MQ=255)
                                                                                                                                      ccAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                                         >  2:1675580/1‑94 (MQ=255)
                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACACGGGACAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCAt                >  1:1136390/1‑150 (MQ=255)
                                                                                                                                             gcTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGACCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCg           <  1:1825377/150‑1 (MQ=255)
                                                                                                                                              cTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCgg          <  1:1426355/150‑1 (MQ=255)
                                                                                                                                               tGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCggg         >  2:2037038/1‑150 (MQ=255)
                                                                                                                                                 gggTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGACTGCAGTGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACAGCCACGGCGCCCATCACCGGGcc       >  1:969297/1‑150 (MQ=255)
                                                                                                                                                   gTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCgg     >  1:1925673/1‑150 (MQ=255)
                                                                                                                                                   gTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCgg     <  1:2723758/150‑1 (MQ=255)
                                                                                                                                                      gCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCaa  <  1:2400186/150‑1 (MQ=255)
                                                                                                                                                      gCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCaa  <  2:888511/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCCGGCAA  >  NC_002947/3337680‑3337978

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