Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,730 A→G pseudogene (1081/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,7300AG99.8% 1988.6 / NA 510pseudogene (1081/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (222/287);  minor base T (0/1);  total (222/288)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 516 total aligned reads displayed.

GGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCC  >  NC_008463/3566581‑3566876
                                                                                                                                                     |                                                                                                                                                    
ggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCgaa                                                                                                                                                    <  2:374278/151‑1 (MQ=255)
ggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCgaa                                                                                                                                                    >  1:4547348/1‑151 (MQ=255)
  aGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                                                                                                  >  2:2054902/1‑151 (MQ=255)
   gCCTGGTCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCCGGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACATCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGtt                                                                                                                                                 <  2:7398284/151‑1 (MQ=255)
   gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCAAGCCGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGACCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGACGtt                                                                                                                                                 <  2:10420990/151‑1 (MQ=255)
     cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                                                                                               >  2:10014814/1‑151 (MQ=255)
          cGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCg                                                                                                                                          <  2:7169159/151‑1 (MQ=255)
            tgctgcTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                                                                                        <  1:7597659/151‑1 (MQ=255)
              ctgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                                                                                      <  2:9220990/151‑1 (MQ=255)
               tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAc                                                                                                                                     <  2:6746185/151‑1 (MQ=255)
                gctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                                                                                                    <  2:10150956/151‑1 (MQ=255)
                gctgAGCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACC‑CCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                                                                                                    <  2:8355192/151‑1 (MQ=255)
                   gATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgag                                                                                                                                 >  1:7853951/1‑148 (MQ=255)
                       gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAccg                                                                                                                             >  2:3106397/1‑150 (MQ=255)
                        aCATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTg                                                                                                                            <  2:8480594/151‑1 (MQ=255)
                         cATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCTCATCGCGCCGGAGCGCGCGGACAACGCGGCGAGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                                                                                           >  2:11111900/1‑151 (MQ=255)
                         cATCCCGCAGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                                                                                           <  2:5828823/151‑1 (MQ=255)
                          aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                                                                                          >  2:10339714/1‑151 (MQ=255)
                          aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                                                                                          >  2:4264326/1‑151 (MQ=255)
                               gCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                                                                                                     <  2:9570096/151‑1 (MQ=255)
                                     ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGag                                                                                                               >  2:9670226/1‑150 (MQ=255)
                                     ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGag                                                                                                               >  2:7871911/1‑150 (MQ=255)
                                       tgtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGagaa                                                                                                             >  2:1092151/1‑151 (MQ=255)
                                         ttCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAGGTGGGAGAACt                                                                                                           >  2:1158404/1‑151 (MQ=255)
                                          tctcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACtt                                                                                                          <  1:3404065/151‑1 (MQ=255)
                                            tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                                                                                        >  2:9062187/1‑151 (MQ=255)
                                             ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCt                                                                                                       >  1:7469055/1‑151 (MQ=255)
                                             ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCt                                                                                                       <  1:9828541/151‑1 (MQ=255)
                                              cTTCATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTg                                                                                                      >  1:802329/1‑151 (MQ=255)
                                                 cATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                                                                                                   <  2:1392502/151‑1 (MQ=255)
                                                 cATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                                                                                                   >  1:6747288/1‑151 (MQ=255)
                                                 cATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                                                                                                   >  2:1871522/1‑151 (MQ=255)
                                                 cATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                                                                                                   <  2:5982868/151‑1 (MQ=255)
                                                 cATGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCT‑GGGCATGTGGGAGAACTTCCTGATcc                                                                                                  <  2:4623615/151‑1 (MQ=255)
                                                  cagCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCACCGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                                                                                                  <  1:10112407/149‑1 (MQ=255)
                                                  aTGCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                                                                                                  >  2:4386568/1‑151 (MQ=255)
                                                    gCAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTg                                                                                                <  2:6844317/151‑1 (MQ=255)
                                                     cAGGACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                                                                                               >  2:859825/1‑151 (MQ=255)
                                                        gACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                                                                                            >  1:3828035/1‑150 (MQ=255)
                                                        gACCTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                                                                                            >  1:3360027/1‑150 (MQ=255)
                                                           cTGCCCAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                                                                                         <  1:1699055/151‑1 (MQ=255)
                                                               ccAGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct                                                                                     <  1:9568653/151‑1 (MQ=255)
                                                                 aGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctct                                                                                   <  2:3975270/151‑1 (MQ=255)
                                                                 aGCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctct                                                                                   <  2:10763310/151‑1 (MQ=255)
                                                                  gCCCCGACTTCCACCA‑CCGCAAGCCGGGGGTGACCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTg                                                                                  <  2:8099241/151‑1 (MQ=255)
                                                                  gCCACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTg                                                                                  >  1:1594400/1‑151 (MQ=255)
                                                                    cACGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc                                                                                <  1:4680997/151‑1 (MQ=255)
                                                                     aCGACTTCCACCA‑CCGCTAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcg                                                                               <  1:2388577/151‑1 (MQ=255)
                                                                     aCGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcg                                                                               <  1:556942/151‑1 (MQ=255)
                                                                      cGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgc                                                                              <  1:4825770/151‑1 (MQ=255)
                                                                      cGACTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgc                                                                              <  1:9035628/151‑1 (MQ=255)
                                                                        aCTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                                                                                               >  1:8616343/1‑132 (MQ=255)
                                                                        aCTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg                                                                            >  1:3601958/1‑151 (MQ=255)
                                                                         cTTCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGc                                                                           <  1:6574002/151‑1 (MQ=255)
                                                                           tCCACCA‑CCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa                                                                         <  2:9834481/151‑1 (MQ=255)
                                                                             cacca‑cCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGt                                                                       <  2:7259088/151‑1 (MQ=255)
                                                                                     gCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAccc                                                                <  1:4432362/151‑1 (MQ=255)
                                                                                     gCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAccc                                                                >  2:2865350/1‑151 (MQ=255)
                                                                                     gCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGAccc                                                                <  2:7323484/151‑1 (MQ=255)
                                                                                      cAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCg                                                               >  1:8068701/1‑151 (MQ=255)
                                                                                         gCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTc                                                            <  1:2639186/151‑1 (MQ=255)
                                                                                           cGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑gg                                                         >  2:2219945/1‑151 (MQ=255)
                                                                                               ggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GG‑aatc                                                    <  2:3847961/151‑4 (MQ=255)
                                                                                                gTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GG‑aatcg                                                   <  1:6223091/151‑5 (MQ=255)
                                                                                                   aaCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAAtcgtcg                                                 <  1:199207/151‑1 (MQ=255)
                                                                                                     cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGAc                                               <  1:718038/151‑1 (MQ=255)
                                                                                                     cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGAc                                               <  1:6318764/151‑1 (MQ=255)
                                                                                                       tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTg                                             <  1:3924271/151‑1 (MQ=255)
                                                                                                       tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTNCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTg                                             >  1:5884486/1‑151 (MQ=255)
                                                                                                         tGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGt                                           >  1:8901899/1‑151 (MQ=255)
                                                                                                         tGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGt                                           >  2:5692967/1‑151 (MQ=255)
                                                                                                               aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACcgcg                                     <  1:10695518/151‑1 (MQ=255)
                                                                                                                  gcgcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCgaga                                  <  1:8298722/151‑1 (MQ=255)
                                                                                                                      cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGTTCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGGTC‑GGAATCGTCGACTGGTACCGCGAGAACCa                              >  2:3190791/1‑151 (MQ=255)
                                                                                                                       ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCAt                             <  2:6314689/151‑1 (MQ=255)
                                                                                                                        gatcgcgcgTCCATCGAGGAGCGCACCTCGATGTTCGGCGCGATGTCCGAGTCCTGGGGCTTGTGGGAGATCTTCCTTATCCTTCTCGACCTCCTCTTCAGCGGCCAGTCGGTCCCGTTCTGA‑ATCGTCGACTGGTACCGCGAGAACCATg                            <  2:3349138/148‑1 (MQ=255)
                                                                                                                         aGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGc                           <  1:10754873/151‑1 (MQ=255)
                                                                                                                           cgcgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGccc                         <  1:10265800/151‑1 (MQ=255)
                                                                                                                             cgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                       >  1:9384733/1‑151 (MQ=255)
                                                                                                                             cgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                       <  1:6038467/151‑1 (MQ=255)
                                                                                                                             cgcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCgg                       >  1:1231319/1‑151 (MQ=255)
                                                                                                                              gcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgc                      <  1:2172048/151‑1 (MQ=255)
                                                                                                                                ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                    <  1:274791/151‑1 (MQ=255)
                                                                                                                                 gCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcg                   <  1:10449772/151‑1 (MQ=255)
                                                                                                                                  ccAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgc                  >  1:5428489/1‑151 (MQ=255)
                                                                                                                                  ccAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgc                  >  2:9435073/1‑151 (MQ=255)
                                                                                                                                   caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                 <  1:10537677/151‑1 (MQ=255)
                                                                                                                                      cgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCg              <  1:10330465/149‑1 (MQ=255)
                                                                                                                                       gaGGCGCGCCCCTCGAAGTTCGGCCCGATGCCCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg             <  1:10960065/150‑1 (MQ=255)
                                                                                                                                       gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg             >  1:5176896/2‑151 (MQ=255)
                                                                                                                                       gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg             <  2:2653988/150‑1 (MQ=255)
                                                                                                                                         ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag           <  1:1445654/151‑1 (MQ=255)
                                                                                                                                          gcgcgcCCCACGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgaga          <  2:8745799/151‑1 (MQ=255)
                                                                                                                                            gcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGAcc        <  2:5176329/151‑1 (MQ=255)
                                                                                                                                              gcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGa      >  2:1628396/1‑151 (MQ=255)
                                                                                                                                               ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt     >  2:7546988/1‑151 (MQ=255)
                                                                                                                                               ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt     <  1:282991/151‑1 (MQ=255)
                                                                                                                                               ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt     <  1:3173744/151‑1 (MQ=255)
                                                                                                                                                 ccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgc   >  1:3113222/1‑151 (MQ=255)
                                                                                                                                                  cTCGAAGTTCGGCCCGATGACCGAGACCAGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgcc  <  2:7635871/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTC‑GG‑ATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCC  >  NC_008463/3566581‑3566876

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