Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,826,373 T→G 5.1% L49R (CTC→CGC)  PP_5102 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,826,3730TG 5.1% 1265.7 / 35.3 365L49R (CTC→CGC) PP_5102hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (185/161);  new base G (6/13);  total (191/174)
Fisher's exact test for biased strand distribution p-value = 9.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 371 total aligned reads displayed.

ACCCATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTAT  >  NC_002947/5826224‑5826522
                                                                                                                                                     |                                                                                                                                                     
  ccATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcg                                                                                                                                                     >  2:3523106/1‑150 (MQ=255)
  ccATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcg                                                                                                                                                     >  1:2653430/1‑150 (MQ=255)
   cATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgc                                                                                                                                                    >  2:3072881/1‑150 (MQ=255)
    aTGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgcg                                                                                                                                                   >  1:3823471/1‑150 (MQ=255)
    aTGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACAGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgcg                                                                                                                                                   <  2:4005628/150‑1 (MQ=255)
        aTCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGgcg                                                                                                                                               <  2:2973316/150‑1 (MQ=255)
            ggTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGccca                                                                                                                                           <  2:2330158/150‑2 (MQ=255)
               aCTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTgcg                                                                                                                                        <  1:8666617/150‑1 (MQ=255)
                  tGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGt                                                                                                                                     <  2:4956246/150‑1 (MQ=255)
                     tCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGcc                                                                                                                                  <  1:7930275/150‑1 (MQ=255)
                       tGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGt                                                                                                                                >  2:4812244/1‑150 (MQ=255)
                        gCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGtt                                                                                                                               <  2:4081852/150‑1 (MQ=255)
                          tCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCa                                                                                                                             >  1:8259407/1‑150 (MQ=255)
                          tCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCa                                                                                                                             >  1:1300086/1‑150 (MQ=255)
                           cgcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCAt                                                                                                                            <  2:5656986/150‑1 (MQ=255)
                            gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                           <  2:5418463/150‑1 (MQ=255)
                            gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                           >  1:3057846/1‑150 (MQ=255)
                             cgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATgg                                                                                                                          <  2:8354064/150‑1 (MQ=255)
                              gACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGgc                                                                                                                         >  2:4557738/1‑150 (MQ=255)
                              gACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGgc                                                                                                                         >  1:6973407/1‑150 (MQ=255)
                              gACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGcccagg                                                                                                                                         <  2:2077481/134‑4 (MQ=255)
                                     tCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATg                                                                                                                  <  1:4001002/150‑1 (MQ=255)
                                        cGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTc                                                                                                               <  1:5692012/150‑1 (MQ=255)
                                         gCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCa                                                                                                              <  2:8721905/150‑1 (MQ=255)
                                          cTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCAt                                                                                                             <  1:7713970/150‑1 (MQ=255)
                                             gAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGAt                                                                                                          <  1:4840832/150‑1 (MQ=255)
                                                ccGCAGCACCCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCgg                                                                                                       >  2:6580786/1‑150 (MQ=255)
                                                 cGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGt                                                                                                      <  1:2025571/150‑1 (MQ=255)
                                                     gcaACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCg                                                                                                  >  1:6366932/1‑150 (MQ=255)
                                                      caACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCgg                                                                                                 <  2:7373848/150‑1 (MQ=255)
                                                             tCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTgcg                                                                                          >  1:7571405/1‑150 (MQ=255)
                                                              cGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTgcgc                                                                                         <  2:2450922/150‑1 (MQ=255)
                                                               ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGTCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                        <  1:9151772/150‑1 (MQ=255)
                                                               ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                        >  1:255590/1‑150 (MQ=255)
                                                               ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                        >  2:7738111/1‑150 (MQ=255)
                                                               ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                        >  1:1481227/1‑150 (MQ=255)
                                                               ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCG‑CAGGGCGCCCGCGAGGCCGcccag                                                                                                                                          >  1:5175878/1‑97 (MQ=37)
                                                               ggAGATTTCACTGGCCGGGGGAGGCGGCCTGCTGGTGCTGTTCGGCACCGATGTCACAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGAGGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                        <  2:5940649/150‑1 (MQ=255)
                                                                      tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTTCGCTGATCTAt                                                                                 <  1:285282/150‑1 (MQ=255)
                                                                      tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                                 <  1:1742991/150‑1 (MQ=255)
                                                                      tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                                 <  1:4534604/150‑1 (MQ=255)
                                                                      tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                                 <  1:4733353/150‑1 (MQ=255)
                                                                        aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                               >  1:167132/1‑150 (MQ=255)
                                                                        aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                               >  2:2111266/1‑150 (MQ=255)
                                                                        aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGATCGGCCCCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                               >  2:5196190/1‑150 (MQ=255)
                                                                              cTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACAcc                                                                         >  2:1721668/1‑150 (MQ=255)
                                                                                tGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgt                                                                       <  1:3718366/150‑1 (MQ=255)
                                                                                    aGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgggt                                                                   <  2:5125765/150‑1 (MQ=255)
                                                                                    aGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgggt                                                                   <  2:2766452/150‑1 (MQ=255)
                                                                                      gCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgggtgg                                                                 <  2:1253746/150‑1 (MQ=255)
                                                                                        tgcctgcTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGgtggtg                                                               >  2:6764264/1‑150 (MQ=255)
                                                                                          cctgcTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGaa                                                             >  1:509824/1‑150 (MQ=255)
                                                                                           ctgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAg                                                            >  1:8310821/1‑150 (MQ=255)
                                                                                            tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                           <  1:6894555/150‑1 (MQ=255)
                                                                                            tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                           <  1:4404989/150‑1 (MQ=255)
                                                                                            tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                           <  2:6122973/150‑1 (MQ=255)
                                                                                             gctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGTGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCAGGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGGAggg                                                          >  1:7084343/1‑150 (MQ=255)
                                                                                             gctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAggg                                                          <  2:7900216/150‑1 (MQ=255)
                                                                                               tgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGc                                                        <  2:703209/150‑1 (MQ=255)
                                                                                                ggTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCt                                                       <  2:1386682/150‑1 (MQ=255)
                                                                                                 gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                      >  2:2666107/1‑150 (MQ=255)
                                                                                                 gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                      >  1:8511532/1‑150 (MQ=255)
                                                                                                 gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                      >  1:8210704/1‑150 (MQ=255)
                                                                                                  tGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCTCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCa                                                     <  2:4908926/150‑1 (MQ=255)
                                                                                                  tGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCa                                                     <  1:7172438/150‑1 (MQ=255)
                                                                                                       ttCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCTCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAg                                                >  1:5425592/1‑150 (MQ=255)
                                                                                                         cGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGtt                                              <  1:1119391/150‑1 (MQ=255)
                                                                                                          ggCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTc                                             <  1:5295109/150‑1 (MQ=255)
                                                                                                            cACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTcaa                                           <  1:8206782/150‑1 (MQ=255)
                                                                                                             aCCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTcaac                                          <  1:9000798/150‑1 (MQ=255)
                                                                                                                gATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCaacaac                                       <  1:8097891/150‑1 (MQ=255)
                                                                                                                 aTGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCaacaaca                                      >  1:5594486/1‑149 (MQ=255)
                                                                                                                      aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTCGCGGTGCTGCGCTGATCTATGCCACACCGGGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                 >  2:5495673/1‑150 (MQ=255)
                                                                                                                      aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                 >  1:601241/1‑150 (MQ=255)
                                                                                                                      aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                 >  1:6106434/1‑150 (MQ=255)
                                                                                                                      aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                 >  1:7531952/1‑150 (MQ=255)
                                                                                                                      aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                 >  2:2157896/1‑150 (MQ=255)
                                                                                                                             gggCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGgcgc                          <  1:4508596/150‑1 (MQ=255)
                                                                                                                              ggCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGcc                         >  1:3309775/1‑150 (MQ=255)
                                                                                                                               gCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCt                        >  2:2396544/1‑150 (MQ=255)
                                                                                                                                 cGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCtgt                      <  2:7587380/150‑1 (MQ=255)
                                                                                                                                  gAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCtgtg                     >  2:6810693/1‑150 (MQ=255)
                                                                                                                                   agcagcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTgtgt                    <  2:90161/150‑1 (MQ=255)
                                                                                                                                   agcagcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTgtgt                    <  1:367805/150‑1 (MQ=255)
                                                                                                                                     cagcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTa                  <  1:4621482/150‑1 (MQ=255)
                                                                                                                                         gcCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtgg              >  1:6237413/1‑150 (MQ=255)
                                                                                                                                           cTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtg            <  2:3866666/150‑1 (MQ=255)
                                                                                                                                           cTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtg            <  2:6441610/150‑1 (MQ=255)
                                                                                                                                            tGGGCGGCCTCACGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtgg           >  1:78562/1‑150 (MQ=255)
                                                                                                                                             gggCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTggtggt          <  2:1192606/150‑1 (MQ=255)
                                                                                                                                               gcggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGc        <  2:2578099/150‑1 (MQ=255)
                                                                                                                                                 ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGc        >  2:7469489/1‑148 (MQ=255)
                                                                                                                                                 ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg      <  2:1881888/150‑1 (MQ=255)
                                                                                                                                                 ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg      <  2:8232167/150‑1 (MQ=255)
                                                                                                                                                  gcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTgg     <  2:465592/150‑1 (MQ=255)
                                                                                                                                                    cTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTa   <  1:3101820/150‑1 (MQ=255)
                                                                                                                                                     tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat  >  1:2190776/1‑150 (MQ=255)
                                                                                                                                                     tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat  >  2:4405530/1‑150 (MQ=255)
                                                                                                                                                     tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat  >  2:7653632/1‑150 (MQ=255)
                                                                                                                                                     tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat  >  2:7986799/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTAT  >  NC_002947/5826224‑5826522

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