Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,399,559 C→A 5.1% A25A (GCC→GCA PP_3002 → shikimate dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,399,5590CA 5.1% 717.7 / 18.6 208A25A (GCC→GCAPP_3002shikimate dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (109/88);  new base A (5/6);  total (114/94)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

CGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTGGCCGACCG  >  NC_002947/3399410‑3399708
                                                                                                                                                     |                                                                                                                                                     
cgGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGcc                                                                                                                                                       <  1:3332659/150‑1 (MQ=255)
     ttggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGc                                                                                                                                                  >  1:55334/1‑150 (MQ=255)
     ttggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGc                                                                                                                                                  >  1:1458047/1‑150 (MQ=255)
     ttggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGAAGGCGTGTTGTTTACCCACGGAGATGCACGATGAGAGACCGCTACGCAGAGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTTTTAGCCCAGGc                                                                                                                                                  >  2:956910/1‑150 (MQ=255)
       ggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                                                                <  2:1204539/150‑1 (MQ=255)
       ggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                                                                <  2:3488701/150‑1 (MQ=255)
       ggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                                                                <  2:3462775/150‑1 (MQ=255)
        gcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAg                                                                                                                                               >  2:219931/1‑150 (MQ=255)
        gcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGGGCGAACAGGCCat                                                                                                                                               <  1:3578843/150‑2 (MQ=255)
          ttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCa                                                                                                                                             >  1:2468147/1‑150 (MQ=255)
           tCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCaa                                                                                                                                            >  1:1096549/1‑150 (MQ=255)
            cAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAAc                                                                                                                                           <  1:3466267/150‑1 (MQ=255)
            cAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAAc                                                                                                                                           <  2:1837997/150‑1 (MQ=255)
             aGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAAcc                                                                                                                                          >  1:666329/1‑150 (MQ=255)
              gCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGGGCGAACAGGCCatgaatca                                                                                                                                         <  1:2119615/150‑8 (MQ=255)
               cccTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCATGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagc                                                                                                                                       <  1:402764/150‑1 (MQ=255)
                ccTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagc                                                                                                                                       >  1:4098359/1‑150 (MQ=255)
                  tGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagcag                                                                                                                                     >  1:3025789/1‑150 (MQ=255)
                         aTCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAAt                                                                                                                              >  1:2555850/1‑150 (MQ=255)
                         aTCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAAt                                                                                                                              >  1:186871/1‑150 (MQ=255)
                              tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCACCCAGCAGCTGGAATACggg                                                                                                                         >  2:3399924/1‑150 (MQ=255)
                              tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACggg                                                                                                                         >  2:3276423/1‑150 (MQ=255)
                              tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACggg                                                                                                                         >  2:1537117/1‑150 (MQ=255)
                              tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACggg                                                                                                                         <  1:2766727/150‑1 (MQ=255)
                              tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATAGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCATGTTCGCACAGGCCAGCAACCAGCAGCTGGAATACggg                                                                                                                         >  2:2818267/1‑150 (MQ=255)
                               gCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACgggg                                                                                                                        <  2:1301721/150‑1 (MQ=255)
                                 cTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATATGGGGCt                                                                                                                      >  2:2710337/1‑150 (MQ=255)
                                        gCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAg                                                                                                               <  2:3880886/150‑1 (MQ=255)
                                          ggATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGc                                                                                                             >  1:1176169/1‑150 (MQ=255)
                                            aTGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTc                                                                                                           >  1:3700787/1‑150 (MQ=255)
                                             tGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCg                                                                                                          >  2:22292/1‑150 (MQ=255)
                                                cGTGTCGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTg                                                                                                       <  1:3887530/150‑1 (MQ=255)
                                                 gtgtTGTTTACCCACGGTGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTgg                                                                                                      <  1:3748296/150‑1 (MQ=255)
                                                 gtgtTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTgg                                                                                                      >  1:1933435/1‑150 (MQ=255)
                                                  tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGga                                                                                                     >  1:1722759/1‑150 (MQ=255)
                                                      gttTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGgacgac                                                                                                 >  1:1332988/1‑150 (MQ=255)
                                                      gttTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGgacgac                                                                                                 >  1:1850335/1‑150 (MQ=255)
                                                        ttACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACtt                                                                                               <  2:254562/150‑1 (MQ=255)
                                                         tACCCTCGGATATGCACGATGAGCGACCTCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGTCTCGCTGGACGACTTc                                                                                              <  2:2910440/150‑1 (MQ=255)
                                                           cccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGa                                                                                            <  2:1774710/150‑1 (MQ=255)
                                                            ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCATCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                                                                           <  1:823588/150‑1 (MQ=255)
                                                            ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                                                                           <  1:1635607/150‑1 (MQ=255)
                                                            ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                                                                           <  2:65202/150‑1 (MQ=255)
                                                              aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGc                                                                                         <  2:3245982/150‑1 (MQ=255)
                                                              aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGc                                                                                         >  1:1394376/1‑150 (MQ=255)
                                                              aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGc                                                                                         >  2:1019347/1‑150 (MQ=255)
                                                              aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGc                                                                                         >  1:4383219/1‑150 (MQ=255)
                                                                ggAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGccc                                                                                       <  2:2224604/150‑1 (MQ=255)
                                                                 gagaTGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCa                                                                                      <  2:3633679/150‑1 (MQ=255)
                                                                   gaTGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAgg                                                                                    <  1:2428317/150‑1 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGgacgac                                                                                                 >  2:2173231/1‑130 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGgacgac                                                                                                 <  1:2173231/130‑1 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCATGTACTGcg                                                                             >  2:4465545/1‑149 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCa                                                                             >  2:1495782/1‑150 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCa                                                                             >  2:1725326/1‑150 (MQ=255)
                                                                          gATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCa                                                                             >  2:3773116/1‑150 (MQ=255)
                                                                            tGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGt                                                                           <  1:2525310/150‑1 (MQ=255)
                                                                                  aCCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTa                                                                     <  1:602572/150‑1 (MQ=255)
                                                                                  aCCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTa                                                                     <  1:2686224/150‑1 (MQ=255)
                                                                                   ccGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAg                                                                    <  2:3567352/150‑1 (MQ=255)
                                                                                      cTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTa                                                                     >  2:609220/1‑146 (MQ=255)
                                                                                       tACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGaa                                                                >  2:851745/1‑150 (MQ=255)
                                                                                          gCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTAGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGc                                                             >  2:39322/1‑150 (MQ=255)
                                                                                            gttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCgg                                                           >  1:4451179/3‑150 (MQ=255)
                                                                                            aGTGATCGGCCGACCGATCAACCTCACCTAGTCCCCGCTGATTCTTGGCCTGTTCGCCCAGGCGAGCTACGAGCTGCTGGAATACGGGGGTTTCGAAGGGTCGCTGGACGACTTCGTAGCCCAGGTACTGCTGTTTGGTAGCGAAGgggg                                                           <  2:368304/150‑1 (MQ=255)
                                                                                            aGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCgg                                                           >  2:4124844/1‑150 (MQ=255)
                                                                                             ttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCggg                                                          >  1:4436525/2‑150 (MQ=255)
                                                                                             gTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCggg                                                          >  2:2097515/1‑150 (MQ=255)
                                                                                             gTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCggg                                                          <  1:2097515/150‑1 (MQ=255)
                                                                                                aTCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAg                                                       >  2:470540/1‑150 (MQ=255)
                                                                                                  cGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAggg                                                     >  2:1571861/1‑150 (MQ=255)
                                                                                                  cGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAggg                                                     <  2:1356168/150‑1 (MQ=255)
                                                                                                      cgaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCGGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATg                                                 <  1:4163024/150‑1 (MQ=255)
                                                                                                       gaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGa                                                >  2:509090/1‑150 (MQ=255)
                                                                                                           gaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTATCGAAGGCGGGAAGGGCATGAACAt                                            >  2:2302303/1‑150 (MQ=255)
                                                                                                           gaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACAt                                            >  2:3571153/1‑150 (MQ=255)
                                                                                                               aaCCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCAcc                                        <  2:3724669/150‑1 (MQ=255)
                                                                                                               aaCCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCAcc                                        <  1:3927231/150‑1 (MQ=255)
                                                                                                               aaCCACACCAACTCCCCGCTGATTCATGGCCTGTTCGCCCACGCCACCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAAGAACATCAcc                                        <  2:3199654/150‑1 (MQ=255)
                                                                                                                   acacCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGccc                                    >  1:3316524/1‑150 (MQ=255)
                                                                                                                      ccAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGt                                 >  1:267243/1‑150 (MQ=255)
                                                                                                                        aaGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTc                               >  1:3973784/1‑150 (MQ=255)
                                                                                                                        aaGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTc                               >  2:753900/1‑150 (MQ=255)
                                                                                                                              ccGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTg                         >  1:354156/1‑150 (MQ=255)
                                                                                                                                gCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCg                       >  1:3351343/1‑150 (MQ=255)
                                                                                                                                   gATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGc                    <  1:1895127/150‑1 (MQ=255)
                                                                                                                                       gctGGCCTGGGCGAACAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTc                >  2:2356179/3‑150 (MQ=255)
                                                                                                                                       cATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTc                <  2:1933435/150‑1 (MQ=255)
                                                                                                                                       cATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTc                <  1:863139/150‑1 (MQ=255)
                                                                                                                                         tGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGa              >  1:277665/1‑150 (MQ=255)
                                                                                                                                         tGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGa              >  2:1250912/1‑150 (MQ=255)
                                                                                                                                         tGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGa              >  1:2661694/1‑150 (MQ=255)
                                                                                                                                         tGGCCTGGGCGAACAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGa              >  2:2426224/1‑150 (MQ=255)
                                                                                                                                             cTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTg          <  2:3308588/150‑1 (MQ=255)
                                                                                                                                             cTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTg          <  1:4062054/150‑1 (MQ=255)
                                                                                                                                             cTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTg          >  1:155161/1‑150 (MQ=255)
                                                                                                                                              tGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGACGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTgg         >  2:1717922/1‑150 (MQ=255)
                                                                                                                                                   gCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTGGCCGAc    >  1:1427124/1‑150 (MQ=255)
                                                                                                                                                   gCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTGGCCGAc    >  2:3028247/1‑150 (MQ=255)
                                                                                                                                                     ccAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTGGCCGACCg  >  2:3665192/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCATGAACATCACCGCCCCGTTCAAGCTGCGTGCCTTCGAGCTGGCCGACCG  >  NC_002947/3399410‑3399708

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: