Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,327,941 T→A 6.5% C317S (TGC→AGC) ‡ PP_2922 → urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,327,9410TA 6.5% 814.3 / 29.5 242C317S (TGC→AGC) ‡PP_2922urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (119/107);  new base A (7/9);  total (126/116)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

GGGGGGCGGCTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGG  >  NC_002947/3327793‑3328090
                                                                                                                                                    |                                                                                                                                                     
        gcTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCACCCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCa                                                                                                                                              <  1:3056705/150‑1 (MQ=255)
          tACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAgg                                                                                                                                            >  1:879591/1‑150 (MQ=255)
          tACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAgg                                                                                                                                            >  1:2047081/1‑150 (MQ=255)
          tACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgc                                                                                                                                            <  1:4886856/150‑2 (MQ=255)
           aCGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGACGCTGGCCTGCCTGGCCAGGc                                                                                                                                           <  2:6815291/150‑1 (MQ=255)
             gCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGctc                                                                                                                                         <  2:911433/150‑1 (MQ=255)
             gCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGctc                                                                                                                                         >  2:475319/1‑150 (MQ=255)
              ccATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGctct                                                                                                                                        <  2:4051718/150‑1 (MQ=255)
              ccATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgcttct                                                                                                                                        <  1:1921283/150‑6 (MQ=255)
              ccATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgcttct                                                                                                                                        <  2:4808604/150‑6 (MQ=255)
               cATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAGcttctt                                                                                                                                       <  1:1482502/150‑1 (MQ=255)
                aTGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTc                                                                                                                                      <  1:2247760/150‑1 (MQ=255)
                    tCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCt                                                                                                                                  <  2:3386090/150‑1 (MQ=255)
                    tCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGGTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAGCCTGGCCAGGCAGGCCAGCttcttcttgt                                                                                                                                  <  2:4472641/150‑5 (MQ=255)
                      ggCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGa                                                                                                                                <  1:1875395/150‑1 (MQ=255)
                      ggCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGa                                                                                                                                <  1:3368893/150‑1 (MQ=255)
                       gCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAg                                                                                                                               >  1:2614164/1‑150 (MQ=255)
                         aCGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGca                                                                                                                             >  2:6204740/1‑150 (MQ=255)
                           ggTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGcatc                                                                                                                           >  1:5182015/1‑150 (MQ=255)
                             tGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGcatcat                                                                                                                         <  1:992759/150‑1 (MQ=255)
                               aTCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCACCCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGt                                                                                                                       <  2:2236287/150‑1 (MQ=255)
                               aTCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgct                                                                                                                                           >  1:3994508/1‑128 (MQ=255)
                                 cAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGttt                                                                                                                     >  1:1267526/1‑150 (MQ=255)
                                  aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                    <  1:6205198/150‑1 (MQ=255)
                                  aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                    <  1:279846/150‑1 (MQ=255)
                                  aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                    <  2:1920694/150‑1 (MQ=255)
                                      ccGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCt                                                                                                                <  1:2444980/150‑1 (MQ=255)
                                       cGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTg                                                                                                               >  1:92053/1‑150 (MQ=255)
                                        gACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCAGCATGTTTTCGCTGt                                                                                                              >  2:5870458/1‑150 (MQ=255)
                                            tCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGcc                                                                                                          >  2:3300015/1‑150 (MQ=255)
                                               aTCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCt                                                                                                       >  1:1869136/1‑150 (MQ=255)
                                                 cTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCACGTTTTCGCTGTAGCCGCTGc                                                                                                     <  2:2646672/150‑1 (MQ=255)
                                                  tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                    >  2:1099443/1‑150 (MQ=255)
                                                  tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                    >  2:1682924/1‑150 (MQ=255)
                                                  tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCAGGTTTTCGCTGTAGCCGCTGcc                                                                                                    >  2:3191075/1‑150 (MQ=255)
                                                   gATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCt                                                                                                   <  1:3657149/150‑1 (MQ=255)
                                                   gATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCt                                                                                                   >  2:4659509/1‑150 (MQ=255)
                                                    aTCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCtt                                                                                                  <  2:1310203/150‑1 (MQ=255)
                                                    aTCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCtt                                                                                                  >  1:537543/1‑150 (MQ=255)
                                                      cGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGc                                                                                                <  2:6173279/150‑1 (MQ=255)
                                                       gggCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCg                                                                                               <  1:3739835/150‑1 (MQ=255)
                                                         gCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGtc                                                                                             >  2:5326666/1‑150 (MQ=255)
                                                           aGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGtctc                                                                                           <  2:855569/150‑1 (MQ=255)
                                                              tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                        <  2:5575896/150‑1 (MQ=255)
                                                              tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                        <  2:723231/150‑1 (MQ=255)
                                                              tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                        <  2:1443651/150‑1 (MQ=255)
                                                              tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                        >  2:3580922/1‑150 (MQ=255)
                                                                cAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgg                                                                                      <  1:4670587/150‑1 (MQ=255)
                                                                cAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgg                                                                                      <  1:2493037/150‑1 (MQ=255)
                                                                  gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                    >  2:5567352/1‑150 (MQ=255)
                                                                  gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                    <  1:1149114/150‑1 (MQ=255)
                                                                  gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                    <  1:1004844/150‑1 (MQ=255)
                                                                  gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                    >  2:4166530/1‑150 (MQ=255)
                                                                     cAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTTCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCg                                                                                 >  2:5216536/1‑150 (MQ=255)
                                                                        ccAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCac                                                                              >  1:4967080/1‑149 (MQ=255)
                                                                        ccAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCaa                                                                              <  2:3938826/150‑1 (MQ=255)
                                                                        ccAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCaa                                                                              >  1:5258497/1‑150 (MQ=255)
                                                                         cagcagGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAAc                                                                             >  2:3077537/1‑150 (MQ=255)
                                                                             agGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGc                                                                         >  1:6373177/1‑150 (MQ=255)
                                                                                   ggTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCa                                                                   <  1:5907344/150‑1 (MQ=255)
                                                                                         tGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTaca                                                             >  2:5415400/1‑150 (MQ=255)
                                                                                         tGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTaca                                                             >  1:1681559/1‑150 (MQ=255)
                                                                                             ggTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTGCGTCCAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCACCAGGCTCCGCAGCTCCCCCGa                                                         >  2:2385612/1‑150 (MQ=255)
                                                                                              gTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAg                                                        <  1:6146280/150‑1 (MQ=255)
                                                                                               tGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGc                                                                         >  2:5790426/1‑132 (MQ=255)
                                                                                               tGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAgc                                                       <  1:5646774/150‑1 (MQ=255)
                                                                                                     tGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAg                                                 <  2:63863/150‑1 (MQ=255)
                                                                                                     tGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAg                                                 >  2:4361767/1‑150 (MQ=255)
                                                                                                          cAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTg                                            <  1:3503355/150‑1 (MQ=255)
                                                                                                           aGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGc                                           <  2:6313005/150‑1 (MQ=255)
                                                                                                           aGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGc                                           <  1:505746/150‑1 (MQ=255)
                                                                                                               aCTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGAt                                       <  2:5690079/150‑1 (MQ=255)
                                                                                                               aCTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGAt                                       >  2:2015922/1‑150 (MQ=255)
                                                                                                               aCTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGAt                                       <  1:5531467/150‑1 (MQ=255)
                                                                                                                cTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATg                                      <  2:6424136/150‑1 (MQ=255)
                                                                                                                cTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATg                                      <  2:6576782/150‑1 (MQ=255)
                                                                                                                      gCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTg                                >  1:5846565/1‑150 (MQ=255)
                                                                                                                      gCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTg                                >  1:6526988/1‑150 (MQ=255)
                                                                                                                       cACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGc                               >  2:2160683/1‑150 (MQ=255)
                                                                                                                        aCGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGcc                              >  1:318089/1‑150 (MQ=255)
                                                                                                                        aCGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGcc                              >  2:1851553/1‑150 (MQ=255)
                                                                                                                         cGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCg                             <  2:2133753/150‑1 (MQ=255)
                                                                                                                          ggCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCgg                            >  1:71142/1‑150 (MQ=255)
                                                                                                                                 tcagctgaagagcCTGGCCAGGCAGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGa                     >  2:3224308/7‑150 (MQ=255)
                                                                                                                                 gTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGa                     <  2:5182009/150‑1 (MQ=255)
                                                                                                                                   aCAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAgg                   >  1:5261665/1‑150 (MQ=255)
                                                                                                                                    cAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGc                  <  2:4987359/150‑1 (MQ=255)
                                                                                                                                     aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCa                                                  <  2:2374250/117‑1 (MQ=255)
                                                                                                                                     aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGcc                 >  1:451986/1‑150 (MQ=255)
                                                                                                                                     aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGcc                 >  1:2715275/1‑150 (MQ=255)
                                                                                                                                          gaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                        >  1:4286239/1‑74 (MQ=255)
                                                                                                                                           agccTGGCCAGGCAGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCa           >  2:6696610/4‑150 (MQ=255)
                                                                                                                                            agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa          <  1:1972436/150‑1 (MQ=255)
                                                                                                                                            agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa          <  1:475319/150‑1 (MQ=255)
                                                                                                                                            agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa          >  1:232493/1‑150 (MQ=255)
                                                                                                                                              cTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAgg        <  1:1094390/150‑1 (MQ=255)
                                                                                                                                                   ctgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTg   <  1:1600959/150‑1 (MQ=255)
                                                                                                                                                   ctgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTg   >  2:2846548/1‑150 (MQ=255)
                                                                                                                                                    tgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTgg  <  2:1783077/150‑1 (MQ=255)
                                                                                                                                                    tgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTgg  >  2:5399565/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGGGGCGGCTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGG  >  NC_002947/3327793‑3328090

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: