Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,616,687 A→C 6.0% T199P (ACC→CCC) ‡ pet → protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,616,6870AC 6.0% 498.7 / 15.8 149T199P (ACC→CCC) ‡petprotein Pet18
Reads supporting (aligned to +/- strand):  ref base A (75/65);  new base C (3/6);  total (78/71)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 152 total aligned reads displayed.

CGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGC  >  NC_002947/3616538‑3616836
                                                                                                                                                     |                                                                                                                                                     
cGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCa                                                                                                                                                       >  2:2437859/1‑150 (MQ=255)
  aaaCCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGACCTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCAcc                                                                                                                                                     <  2:3231103/150‑1 (MQ=255)
       gATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCACGTGACTTCTACCAGCGCACCGTGgc                                                                                                                                                >  2:3440841/1‑150 (MQ=255)
           ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCCGCGCACCGTGGCGCTg                                                                                                                                            <  1:415590/150‑1 (MQ=255)
           ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTg                                                                                                                                            <  2:1120897/150‑1 (MQ=255)
            cGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTgg                                                                                                                                           >  1:265551/1‑150 (MQ=255)
               aCAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGc                                                                                                                                        <  2:2056115/150‑1 (MQ=255)
               aCAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGc                                                                                                                                        <  1:3983172/150‑1 (MQ=255)
                 aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaa                                                                                                                                      <  1:1490975/150‑6 (MQ=255)
                 aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaa                                                                                                                                      <  2:3817631/150‑6 (MQ=255)
                   cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGc                                                                                                                                    >  2:4625304/1‑150 (MQ=255)
                   cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGc                                                                                                                                    >  1:2003425/1‑150 (MQ=255)
                    ttCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCt                                                                                                                                   <  2:4945485/150‑1 (MQ=255)
                       gTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTAGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttc                                                                                                                                >  2:4297791/1‑150 (MQ=255)
                        tCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttct                                                                                                                               >  2:4493987/1‑150 (MQ=255)
                          cATGCCGAGTGGATGACCCTGCATGACAACCCGGCAGTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGGGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttcttc                                                                                                                             <  1:2289572/150‑1 (MQ=255)
                            tGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttctttgc                                                                                                                           <  1:5008977/150‑4 (MQ=255)
                            tGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGa                                                                                                                           <  2:4112444/150‑1 (MQ=255)
                            tGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGa                                                                                                                           <  2:2105177/150‑1 (MQ=255)
                                gAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCt                                                                                                                       >  1:1384836/1‑150 (MQ=255)
                                  gTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtgt                                                                                                                     <  2:3794609/150‑1 (MQ=255)
                                   tGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtgtg                                                                                                                    >  1:3460941/1‑150 (MQ=255)
                                      aTTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTAc                                                                                                                 >  2:4506493/1‑150 (MQ=255)
                                          cccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCgg                                                                                                             >  2:2364225/1‑150 (MQ=255)
                                           ccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGa                                                                                                            >  1:3072643/1‑150 (MQ=255)
                                           ccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGa                                                                                                            >  1:1679595/1‑150 (MQ=255)
                                            cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                           <  2:2654731/150‑1 (MQ=255)
                                            cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                           <  2:2309078/150‑1 (MQ=255)
                                             tGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAgg                                                                                                          >  2:2530273/1‑150 (MQ=255)
                                             tGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAgg                                                                                                          <  2:4531646/150‑1 (MQ=255)
                                                   aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAAt                                                                                                    >  1:4065828/1‑150 (MQ=255)
                                                   aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAAt                                                                                                    >  1:4065855/1‑150 (MQ=255)
                                                   aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAAt                                                                                                    <  1:704502/150‑1 (MQ=255)
                                                     aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGa                                                                                                  <  1:2555537/150‑1 (MQ=255)
                                                       cccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGc                                                                                                >  1:1326385/1‑150 (MQ=255)
                                                        ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               >  2:4308577/1‑150 (MQ=255)
                                                        ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               <  2:1851006/150‑1 (MQ=255)
                                                        ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               >  2:1772791/1‑150 (MQ=255)
                                                        ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               >  1:1032253/1‑150 (MQ=255)
                                                          ggCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                             >  1:2343567/1‑150 (MQ=255)
                                                          ggCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                             <  1:1370670/150‑1 (MQ=255)
                                                            cATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGa                                                                                           <  1:1571033/150‑1 (MQ=255)
                                                             aTTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGat                                                                                          <  2:2189936/150‑1 (MQ=255)
                                                                  aGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCtt                                                                                     >  2:4625552/1‑150 (MQ=255)
                                                                  aGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCtt                                                                                     >  1:903077/1‑150 (MQ=255)
                                                                   ggCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                    <  1:2976583/150‑1 (MQ=255)
                                                                   ggCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                    <  1:2853517/150‑1 (MQ=255)
                                                                    gCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCg                                                                                   >  2:2981280/1‑150 (MQ=255)
                                                                    gCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCg                                                                                   >  2:3907924/1‑150 (MQ=255)
                                                                     cATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGa                                                                                  <  1:4824723/150‑1 (MQ=255)
                                                                     cATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGa                                                                                  <  1:1785677/150‑1 (MQ=255)
                                                                          gtcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                             >  2:4323351/1‑150 (MQ=255)
                                                                           tcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCt                                                                            >  2:4486968/1‑150 (MQ=255)
                                                                               tttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGGCTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAg                                                                        >  1:3788322/1‑150 (MQ=255)
                                                                               tttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAg                                                                        <  2:2081218/150‑1 (MQ=255)
                                                                                   cTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCg                                                                    >  2:3174694/1‑150 (MQ=255)
                                                                                   cTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCg                                                                    >  2:2940879/1‑150 (MQ=255)
                                                                                         gCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgcca                                                              >  1:5035565/1‑150 (MQ=255)
                                                                                            gAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccacgc                                                           <  1:4228636/150‑1 (MQ=255)
                                                                                                tGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCcaca                                                       >  1:711928/1‑150 (MQ=255)
                                                                                                tGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCcaca                                                       >  2:2811165/1‑150 (MQ=255)
                                                                                                 ggACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACaa                                                      <  2:1940938/150‑1 (MQ=255)
                                                                                                 ggACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACaa                                                      <  1:2629975/150‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  1:560998/150‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  2:3626240/150‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  2:4381776/150‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  1:4441517/150‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  1:4441545/150‑1 (MQ=255)
                                                                                                       cgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAAc                                                <  2:3139009/150‑1 (MQ=255)
                                                                                                        gTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAAcc                                               <  1:4667363/150‑1 (MQ=255)
                                                                                                          cGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGc                                             <  1:2298622/150‑1 (MQ=255)
                                                                                                          cGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCACTGCCACGCCACAATGGAACCGc                                             <  2:4947521/150‑1 (MQ=255)
                                                                                                            gCCCGCAGGAGGCCGGTTTGGCAAGTGACTTCTAACAGCGCACATTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTAGACGGAGGGCAAATGAGAATGGATATCTTAGATCGCCTAAAGGCCGCTGACACGCCACAATGGAACCGCTa                                           >  2:2874922/1‑150 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACg                                         >  1:2936600/1‑150 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACg                                         >  1:4924808/1‑150 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACg                                         >  2:1229487/1‑150 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACg                                         >  1:378691/1‑150 (MQ=255)
                                                                                                                gCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTa                                       >  1:3030895/1‑150 (MQ=255)
                                                                                                                gCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTa                                       >  2:4709527/1‑150 (MQ=255)
                                                                                                                   ggaggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGAc                                    <  2:4541180/150‑1 (MQ=255)
                                                                                                                    gaggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGAcc                                   >  1:3938760/1‑150 (MQ=255)
                                                                                                                     aggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCa                                  <  1:4506493/150‑1 (MQ=255)
                                                                                                                          gCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACt                             <  2:4846919/150‑1 (MQ=255)
                                                                                                                            tttGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTc                           >  2:2557860/1‑150 (MQ=255)
                                                                                                                            tttGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTc                           >  1:4547598/1‑150 (MQ=255)
                                                                                                                                  aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                     >  2:4490913/1‑150 (MQ=255)
                                                                                                                                       tcTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGAt                <  1:163136/149‑1 (MQ=255)
                                                                                                                                         ttCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATgg              <  2:264329/150‑1 (MQ=255)
                                                                                                                                           gctccAGCGCCACGGTGCGCTGGTGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAAGGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACACTGGAACCGCTACGTAGACCAGGACTTCGTGCGGCAGATGGgc            >  2:3986082/4‑150 (MQ=255)
                                                                                                                                            tACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcg           >  1:1015004/1‑150 (MQ=255)
                                                                                                                                             tccAGCGCCACGGTGCGCTGGAGCTGGGTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcgc          >  2:569633/2‑150 (MQ=255)
                                                                                                                                             tccAGCGCCACGGTGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcgc          >  2:569286/2‑150 (MQ=255)
                                                                                                                                              ccAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCa         <  2:2289807/150‑1 (MQ=255)
                                                                                                                                              ccAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCa         <  1:4625304/150‑1 (MQ=255)
                                                                                                                                              ccAGCGCACCGTGGCGCGGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCa         >  1:3066699/1‑150 (MQ=255)
                                                                                                                                               cAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAg        <  2:3899174/150‑1 (MQ=255)
                                                                                                                                                 gcgcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGc      >  1:3024500/1‑150 (MQ=255)
                                                                                                                                                   gcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCAc    >  1:534715/1‑150 (MQ=255)
                                                                                                                                                   gcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCAc    <  2:745325/150‑1 (MQ=255)
                                                                                                                                                     aCCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGc  <  2:2066922/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGC  >  NC_002947/3616538‑3616836

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