Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,327,945 T→A 6.5% L318Q (CTG→CAG)  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,9450TA 6.5% 812.8 / 29.5 242L318Q (CTG→CAG) 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 261 total aligned reads displayed.

GGGCGGCTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGGCCGT  >  NC_002947/3327796‑3328094
                                                                                                                                                     |                                                                                                                                                     
     gcTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCa                                                                                                                                                  <  2:5747420/150‑1 (MQ=255)
      cTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAg                                                                                                                                                 >  1:2670820/1‑150 (MQ=255)
       tACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAgg                                                                                                                                                >  1:879591/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)
           ccATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGctct                                                                                                                                            <  2:5859470/150‑1 (MQ=255)
           ccATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgcttct                                                                                                                                            <  2:4808604/150‑6 (MQ=255)
            cATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAGcttctt                                                                                                                                           <  1:1482502/150‑1 (MQ=255)
              tGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCa                                                                                                                                         >  1:5209253/1‑150 (MQ=255)
                 tCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCt                                                                                                                                      >  1:304667/1‑150 (MQ=255)
                   ggCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGa                                                                                                                                    <  1:3368893/150‑1 (MQ=255)
                    gCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAg                                                                                                                                   >  2:2257874/1‑150 (MQ=255)
                    gCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAg                                                                                                                                   >  1:2614164/1‑150 (MQ=255)
                      aCGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGca                                                                                                                                 >  2:6204740/1‑150 (MQ=255)
                        ggTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGcatc                                                                                                                               >  1:5182009/1‑150 (MQ=255)
                          tGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGcatcat                                                                                                                             <  1:992759/150‑1 (MQ=255)
                            aTCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgct                                                                                                                                               >  1:3994508/1‑128 (MQ=255)
                            aTCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCAGGCAGGCCAgct                                                                                                                                               <  2:3994508/130‑3 (MQ=255)
                              cAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGttt                                                                                                                         >  1:1267526/1‑150 (MQ=255)
                               aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                        <  1:279846/150‑1 (MQ=255)
                               aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                        <  1:6205260/150‑1 (MQ=255)
                               aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                        >  1:161625/1‑150 (MQ=255)
                               aGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGtttt                                                                                                                        <  2:1920694/150‑1 (MQ=255)
                                  gcCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGc                                                                                                                     <  1:6643507/150‑1 (MQ=255)
                                         tCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGcc                                                                                                              >  2:3300015/1‑150 (MQ=255)
                                           gATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGc                                                                                                            >  2:178098/1‑150 (MQ=255)
                                           gATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGc                                                                                                            >  2:2756922/1‑150 (MQ=255)
                                            aTCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCt                                                                                                           >  1:1869136/1‑150 (MQ=255)
                                               tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                        >  2:4761007/1‑150 (MQ=255)
                                               tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                        >  2:3365286/1‑150 (MQ=255)
                                               tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                        >  2:1682924/1‑150 (MQ=255)
                                               tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGcc                                                                                                        <  1:202279/150‑1 (MQ=255)
                                               tGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCAGGTTTTCGCTGTAGCCGCTGcc                                                                                                        >  2:3191075/1‑150 (MQ=255)
                                                gATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCt                                                                                                       >  2:4659509/1‑150 (MQ=255)
                                                gATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCt                                                                                                       <  1:3657149/150‑1 (MQ=255)
                                                 aTCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCTTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCta                                                                                                      >  1:944614/1‑149 (MQ=255)
                                                 aTCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCtt                                                                                                      <  2:1310203/150‑1 (MQ=255)
                                                   cGGGCAGAGGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCGGCCTTGc                                                                                                    <  2:6185815/150‑1 (MQ=255)
                                                    gggCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCg                                                                                                   <  1:3739835/150‑1 (MQ=255)
                                                        aGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGtctc                                                                                               <  2:855569/150‑1 (MQ=255)
                                                           tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            >  2:3580922/1‑150 (MQ=255)
                                                           tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            <  2:723231/150‑1 (MQ=255)
                                                           tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            <  2:5575896/150‑1 (MQ=255)
                                                           tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            <  2:6420182/150‑1 (MQ=255)
                                                           tGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            <  2:1443651/150‑1 (MQ=255)
                                                            gCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACg                                                                                           >  1:1881904/1‑150 (MQ=255)
                                                             cAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgg                                                                                          <  1:2493037/150‑1 (MQ=255)
                                                             cAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgg                                                                                          <  1:4670587/150‑1 (MQ=255)
                                                               gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                        >  2:5567437/1‑150 (MQ=255)
                                                               gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                        <  1:1004844/150‑1 (MQ=255)
                                                               gCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACgggg                                                                                        <  1:1149114/150‑1 (MQ=255)
                                                                  cAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTTCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCg                                                                                     >  2:5216536/1‑150 (MQ=255)
                                                                     ccAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAATCTGGCCTGCCTGGCCAGGCTATTCAGCTGAGCATCATGTTATCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCaa                                                                                  >  1:5258631/1‑150 (MQ=255)
                                                                     ccAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCac                                                                                  >  1:4967080/1‑149 (MQ=255)
                                                                      cagcagGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAAc                                                                                 >  2:3077537/1‑150 (MQ=255)
                                                                          agGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGc                                                                             >  1:6373177/1‑150 (MQ=255)
                                                                              gcgGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCg                                                                         >  1:5335991/1‑150 (MQ=255)
                                                                                ggTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCa                                                                       <  1:5907344/150‑1 (MQ=255)
                                                                                ggTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCa                                                                       <  2:6349888/150‑1 (MQ=255)
                                                                                     cTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTac                                                                  <  2:4551665/150‑1 (MQ=255)
                                                                                      tGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTaca                                                                 >  1:1681559/1‑150 (MQ=255)
                                                                                           gTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAg                                                            <  1:6146280/150‑1 (MQ=255)
                                                                                              aGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAgcgc                                                         <  2:4118949/150‑1 (MQ=255)
                                                                                                ccTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGACGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAgcgcgc                                                       <  2:1119097/150‑1 (MQ=255)
                                                                                                ccTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCACAACAGGCTCCGCAGCTACACCGAgcgcgc                                                       >  2:5966979/1‑150 (MQ=255)
                                                                                                    gAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGc                                                   <  2:2483798/150‑1 (MQ=255)
                                                                                                       cAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTg                                                <  1:3503355/150‑1 (MQ=255)
                                                                                                        aGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGc                                               <  1:505746/150‑1 (MQ=255)
                                                                                                            aCTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGAt                                           >  2:2015922/1‑150 (MQ=255)
                                                                                                             cTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATg                                          <  2:6576782/150‑1 (MQ=255)
                                                                                                                   gCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTg                                    >  1:5393328/1‑150 (MQ=255)
                                                                                                                   gCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTg                                    >  1:5846565/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)
                                                                                                                       ggCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCgg                                >  1:71142/1‑150 (MQ=255)
                                                                                                                         cGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGtt                              <  1:1935710/150‑1 (MQ=255)
                                                                                                                              tcagctgaagagcCTGGCCAGGCAGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGa                         >  2:3224308/7‑150 (MQ=255)
                                                                                                                              gTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGa                         <  2:5182015/150‑1 (MQ=255)
                                                                                                                               tACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAg                        >  2:6461673/1‑150 (MQ=255)
                                                                                                                                aCAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAgg                       <  1:2093015/150‑1 (MQ=255)
                                                                                                                                 cAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGc                      <  2:4987359/150‑1 (MQ=255)
                                                                                                                                 cAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGc                      >  2:213407/1‑150 (MQ=255)
                                                                                                                                  aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCa                                                      >  1:2374250/1‑117 (MQ=255)
                                                                                                                                  aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCa                                                      <  2:2374250/117‑1 (MQ=255)
                                                                                                                                  aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGcc                     >  2:5334116/1‑150 (MQ=255)
                                                                                                                                  aagaagaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGcc                     >  1:451986/1‑150 (MQ=255)
                                                                                                                                       gaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            >  1:4286239/1‑74 (MQ=255)
                                                                                                                                       gaagCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGAc                                                                                            <  2:4286239/74‑1 (MQ=255)
                                                                                                                                         agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa              <  1:475319/150‑1 (MQ=255)
                                                                                                                                         agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa              >  1:232493/1‑150 (MQ=255)
                                                                                                                                         agCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCaa              <  1:5641927/150‑1 (MQ=255)
                                                                                                                                           cTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAgg            <  1:1094390/150‑1 (MQ=255)
                                                                                                                                             ggCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGc          >  2:5507982/1‑150 (MQ=255)
                                                                                                                                                ctgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTg       >  2:2846548/1‑150 (MQ=255)
                                                                                                                                                 tgcctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTgg      >  2:5399643/1‑150 (MQ=255)
                                                                                                                                                   cctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGGcc    <  2:6599437/150‑1 (MQ=255)
                                                                                                                                                   cctgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGGcc    <  2:3882717/150‑1 (MQ=255)
                                                                                                                                                     tgGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGGCCGt  >  2:4455248/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGGCTACGCCATGATCGGCACGGTGATCAGCGCCGACCTCGATCTGATCGGGCAGATGCAGCCCAACCAGCAGGCGCGGTTCCTGCCGGTGAGCCTGGAGCAGGCACTGGAGGCACGGCGTTCGTACAAGAAGAAGCTGGCCTGCCTGGCCAGGCTCTTCAGCTGAGCATCATGTTTTCGCTGTAGCCGCTGCCTTGCGTCTCGACGGGGCCGCAACAGGCTCCGCAGCTACACCGAGCGCGCAGGCTTGCTGATGGTCTTGCCGGTTTGCGAGGCCACCTCAAGGGCTTGGCCGT  >  NC_002947/3327796‑3328094

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: ‑