Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,616,688 C→A 6.1% T199N (ACC→AAC) ‡ pet → protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,616,6880CA 6.1% 493.5 / 14.6 148T199N (ACC→AAC) ‡petprotein Pet18
Reads supporting (aligned to +/- strand):  ref base C (75/64);  new base A (3/6);  total (78/70)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 150 total aligned reads displayed.

AAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGC  >  NC_002947/3616540‑3616836
                                                                                                                                                    |                                                                                                                                                    
aaaCCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGACCTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCAcc                                                                                                                                                     <  2:3231103/150‑1 (MQ=255)
     gATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCACGTGACTTCTACCAGCGCACCGTGgc                                                                                                                                                >  2:3440841/1‑150 (MQ=255)
       tGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGgcgc                                                                                                                                              <  2:4690188/150‑1 (MQ=255)
         ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTg                                                                                                                                            <  2:1120924/150‑1 (MQ=255)
         ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTg                                                                                                                                            <  2:1120897/150‑1 (MQ=255)
         ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTg                                                                                                                                            <  2:4863647/150‑1 (MQ=255)
         ccGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTg                                                                                                                                            <  1:4049286/150‑1 (MQ=255)
          cGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTgg                                                                                                                                           >  1:265551/1‑150 (MQ=255)
            gACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggta                                                                                                                                         <  1:4601599/150‑3 (MQ=255)
            gACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAg                                                                                                                                         >  2:2251750/1‑150 (MQ=255)
             aCAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGc                                                                                                                                        <  1:3983172/150‑1 (MQ=255)
               aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaa                                                                                                                                      <  2:3817631/150‑6 (MQ=255)
               aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaa                                                                                                                                      <  1:1490975/150‑6 (MQ=255)
               aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTg                                                                                                                                      >  2:3976210/1‑150 (MQ=255)
                aCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTgg                                                                                                                                     <  2:3866710/150‑1 (MQ=255)
                 cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGc                                                                                                                                    >  1:78829/1‑150 (MQ=255)
                 cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGc                                                                                                                                    >  1:2003425/1‑150 (MQ=255)
                     gTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTAGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttc                                                                                                                                >  2:4297791/1‑150 (MQ=255)
                     gTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTAGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttc                                                                                                                                >  2:4297757/1‑150 (MQ=255)
                        cATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttcttc                                                                                                                             <  2:451522/150‑1 (MQ=255)
                        cATGCCGAGTGGATGACCCTGCATGACAACCCGGCAGTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGGGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTttcttc                                                                                                                             <  1:2289572/150‑1 (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:2063929/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:2309082/150‑1 (MQ=255)
                                          cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                           <  2:2654731/150‑1 (MQ=255)
                                          cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                           <  2:4599292/150‑1 (MQ=255)
                                           tGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAgg                                                                                                          >  2:2530273/1‑150 (MQ=255)
                                                 aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAAt                                                                                                    >  1:4065855/1‑150 (MQ=255)
                                                   aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGa                                                                                                  <  1:2555537/150‑1 (MQ=255)
                                                     cccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGc                                                                                                >  1:1326385/1‑150 (MQ=255)
                                                      ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               >  2:1772791/1‑150 (MQ=255)
                                                      ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCa                                                                                               >  1:1032253/1‑150 (MQ=255)
                                                       cGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCAt                                                                                              <  1:690208/150‑1 (MQ=255)
                                                        ggCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                             <  1:1370670/150‑1 (MQ=255)
                                                        ggCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                             >  1:2343567/1‑150 (MQ=255)
                                                          cATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGa                                                                                           <  1:1571033/150‑1 (MQ=255)
                                                            ttCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGata                                                                                         >  1:4302472/1‑150 (MQ=255)
                                                                aGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCtt                                                                                     >  2:4625552/1‑150 (MQ=255)
                                                                aGGCATTCGTCGCTATCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGAGACTTCTTCCAGCGCACCGTGGCGCTGGAGCTCGCTTTCTTCGATCCTGTGTACCCGGAGGGCAAATGAGCATGGATATCta                                                                                     >  2:1875372/1‑149 (MQ=255)
                                                                 ggCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                    <  1:2976583/150‑1 (MQ=255)
                                                                 ggCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                                    <  1:2853517/150‑1 (MQ=255)
                                                                  gCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCg                                                                                   >  2:3907924/1‑150 (MQ=255)
                                                                  gCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCg                                                                                   >  2:2981280/1‑150 (MQ=255)
                                                                   cATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGa                                                                                  <  1:4493987/150‑1 (MQ=255)
                                                                   cATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGa                                                                                  <  1:1785677/150‑1 (MQ=255)
                                                                             tttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGGCTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAg                                                                        >  1:3788322/1‑150 (MQ=255)
                                                                             tttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAg                                                                        <  2:2081232/150‑1 (MQ=255)
                                                                              ttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAgg                                                                       <  1:5110472/150‑1 (MQ=255)
                                                                                 cTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCg                                                                    >  2:3174694/1‑150 (MQ=255)
                                                                                 cTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCg                                                                    >  2:2940879/1‑150 (MQ=255)
                                                                                  tGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGc                                                                   >  1:4133786/1‑150 (MQ=255)
                                                                                   gCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAATTGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCt                                                                  >  2:5080569/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                                                      <  1:2629975/150‑1 (MQ=255)
                                                                                               ggACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACaa                                                      <  2:1940938/150‑1 (MQ=255)
                                                                                                gACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAAt                                                     <  2:4492700/150‑1 (MQ=255)
                                                                                                  ccGCGGCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACGTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGGGTACCCGGAGGGCAAATGAGCATGGAGATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCCCCATgg                                                   >  2:3787758/1‑150 (MQ=255)
                                                                                                   cgcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGa                                                  >  2:4381967/1‑150 (MQ=255)
                                                                                                    gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  1:560998/150‑1 (MQ=255)
                                                                                                    gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  1:4441545/150‑1 (MQ=255)
                                                                                                    gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 >  1:382760/1‑150 (MQ=255)
                                                                                                    gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  2:4381776/150‑1 (MQ=255)
                                                                                                    gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGaa                                                 <  2:3626240/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)
                                                                                                              gCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTa                                       >  1:3030895/1‑150 (MQ=255)
                                                                                                                 ggaggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGAc                                    <  2:4541180/150‑1 (MQ=255)
                                                                                                                    ggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCAt                                 <  1:4936807/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:4311921/1‑150 (MQ=255)
                                                                                                                                aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg                     >  2:4490913/1‑150 (MQ=255)
                                                                                                                                     tcTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGAt                <  1:163136/149‑1 (MQ=255)
                                                                                                                                       ttCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATgg              >  1:3659901/1‑150 (MQ=255)
                                                                                                                                       ttCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATgg              <  2:264329/150‑1 (MQ=255)
                                                                                                                                          tACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcg           >  1:1015004/1‑150 (MQ=255)
                                                                                                                                           tccAGCGCCACGGTGCGCTGGAGCTGGGTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgcgc          >  2:569633/2‑150 (MQ=255)
                                                                                                                                            ccAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCa         <  2:2289807/150‑1 (MQ=255)
                                                                                                                                            ccAGCGCACCGTGGCGCGGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCa         >  1:3066699/1‑150 (MQ=255)
                                                                                                                                             cAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAg        >  2:1452597/1‑150 (MQ=255)
                                                                                                                                               gcgcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGc      >  1:469170/1‑150 (MQ=255)
                                                                                                                                               gcgcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGc      >  1:3024500/1‑150 (MQ=255)
                                                                                                                                                 gcACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCAc    <  2:745325/150‑1 (MQ=255)
                                                                                                                                                   aCCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGc  <  2:2066922/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAGGCACGC  >  NC_002947/3616540‑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: ‑