Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,394,501 T→A 7.0% T411T (ACA→ACT PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,394,5010TA 7.0% 952.7 / 13.5 336T411T (ACA→ACTPP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base T (127/162);  new base A (27/4);  total (168/168)
Fisher's exact test for biased strand distribution p-value = 4.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 348 total aligned reads displayed.

TAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAG  >  NC_002947/4394352‑4394650
                                                                                                                                                     |                                                                                                                                                     
tAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACt                                                                                                                                                       <  2:353687/150‑1 (MQ=255)
   gCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACtgtg                                                                                                                                                    <  2:1427893/150‑1 (MQ=255)
        gATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgc                                                                                                                                               <  2:7372640/150‑1 (MQ=255)
         aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCCCGGCGCTCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCCCGGTGGCcgcg                                                                                                                                              >  2:5083479/1‑150 (MQ=255)
          tCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCCCCCGCGTGGCCACTGTGGCCGCGa                                                                                                                                             >  2:4259949/1‑150 (MQ=255)
          tCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCCCGAGCCCCCCGGCGATCAGGCCAAACGGCGGGGCCCCCCAGGCCCCACCCCGCACCCCCGGGGCCACAGTGGCCGCGa                                                                                                                                             >  1:6597133/1‑150 (MQ=255)
           cAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGaa                                                                                                                                            <  2:1823738/150‑1 (MQ=255)
                cccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccg                                                                                                                                       >  2:2393118/1‑150 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCGGCCGCGCTCCCCCCGGCGCTCAGGCCCAAAGGCGTGGGCGCGCCGGCCCCACCCCGCACCCCCGGGGGCCCAGTGGCCGCCGACAcccc                                                                                                                                      >  1:1944532/1‑148 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccgc                                                                                                                                      <  1:4459303/150‑1 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccgc                                                                                                                                      >  2:6604357/1‑150 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAccgc                                                                                                                                      <  1:7077293/150‑1 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCCGCGATCAGGCCAATCGGCGTGGCCGCGCAGGCCCCACCGAGCCCCCGCGTGGCCACGGTGGCCGCGCACAcccc                                                                                                                                      >  2:3400830/1‑148 (MQ=255)
                 ccACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGAACCCCCCGGCGCTAAGGCCAAAACGCGGGGGCGCCAAGGCCACACCCCGCACCCCCGGGGCCCCAGTGGCCGCGAAAAcccc                                                                                                                                      >  1:2196786/1‑148 (MQ=255)
                  cacaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCCCCCCGAAACACCCCCCGAACAGCACAAACCCCCGTGCCCCCAAGGCCACACCCACCACCCCCCTGGCCACAATGTCCGCCAACAccccc                                                                                                                                     >  1:3000118/1‑150 (MQ=255)
                     aGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCcgca                                                                                                                                  >  2:2378254/1‑150 (MQ=255)
                           cTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAgc                                                                                                                            <  2:2798351/150‑1 (MQ=255)
                           cTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCCATCCCCCCGGCGATCAGGCCAATCGGCGTGGCCCCGAAGGCCCCACCCCGCACCCCCGTGGCCCCGGTGGCCGCGAACACCCCcccaccccgc                                                                                                                            >  2:2185774/1‑145 (MQ=255)
                            tGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCCCTGTGGCCGCGAACACCGCCGCACGCAgcg                                                                                                                           >  2:7205137/1‑150 (MQ=255)
                               agaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCgg                                                                                                                        <  2:3189590/150‑1 (MQ=255)
                                   cGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                    <  1:6172659/150‑1 (MQ=255)
                                    ggCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                                   <  2:2638516/150‑1 (MQ=255)
                                    ggCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                                   <  1:6849461/150‑1 (MQ=255)
                                    ggCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGAGCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                                   <  1:2765199/150‑1 (MQ=255)
                                             ttCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCg                                                                                                          <  2:5358391/150‑1 (MQ=255)
                                             ttCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGGGAAGGCCGCCCCg                                                                                                          >  2:5506099/1‑150 (MQ=255)
                                               ctcCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCCCGGCGCTCAGGCCAAACGGCGTGGCCGCGAAGGCCCCCCCCAGCACCCCCGTGGCCACAGTGGCCGCGCACACCCCCGCACGCAGCGCGGTGAAGGCccccccccc                                                                                                        >  1:5449953/1‑147 (MQ=255)
                                                 cccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                                                                      <  2:4172205/150‑1 (MQ=255)
                                                  ccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAt                                                                                                     <  1:5556850/150‑1 (MQ=255)
                                                     gTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGAc                                                                                                  <  2:1175396/150‑1 (MQ=255)
                                                        ggCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGAt                                                                                               <  1:2132045/150‑1 (MQ=255)
                                                           gacgaTCAGACCGGCCGCCGCCGCGATCCCCGCGGCGCTCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACGGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACCATGcc                                                                                            >  2:6535993/1‑150 (MQ=255)
                                                           gacgaTCAGACCGGCCGCCGCCGCGATCCCCCCGGCGATCAGGCCAATCGGCGTGGCCGCGAACGCCACACCCAGCACCCCCGTGGCCACTGTGGCCGCGAACACCGCCCCACGCACCCCGGTGAAGGCCCCCCCCGCAAAGACGAAGcc                                                                                            >  2:6273459/1‑150 (MQ=255)
                                                               aTCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCgcgc                                                                                        <  1:227442/150‑1 (MQ=255)
                                                                  aGACCGGCCGCCGCCGCGATCCCCCCCGCGCTCAGGCCAAACGGCGTGGCCCCGAACGCCCCACCCAGCACCCCCGTGGCCACAGGGGCCCCGAACACCGCCCCACCCACCCCGGGGAAGGCCGCCCCCGCCAAGAACATGCCCCGCAcc                                                                                     >  2:5700078/1‑150 (MQ=255)
                                                                      cGGCCGCCGCCGCGATACCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCCAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCCCACGCAGCGCGGTGAAGGCCGCCCCGGCAAAGACGATGCCCCGCACCAAca                                                                                 >  2:6797103/1‑149 (MQ=255)
                                                                        gccgccgccGCGATCCCCGCGGCGATCAGGCCAAAAGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                                                               >  2:1255925/1‑150 (MQ=255)
                                                                         ccgccgccgCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCg                                                                              <  1:4603766/150‑1 (MQ=255)
                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  1:3843956/150‑1 (MQ=255)
                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  1:501643/150‑1 (MQ=255)
                                                                          cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                             <  2:4918844/150‑1 (MQ=255)
                                                                           gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCggg                                                                            <  1:3387949/150‑1 (MQ=255)
                                                                           gccgccgcGATCCCCGCGGCGAGCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCggg                                                                            >  2:6111513/1‑150 (MQ=255)
                                                                             cgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                          <  2:1694447/150‑1 (MQ=255)
                                                                              gccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTg                                                                         >  1:3157605/1‑150 (MQ=255)
                                                                              gccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTg                                                                         >  1:7128371/1‑150 (MQ=255)
                                                                                cgcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTc                                                                       >  1:1758537/1‑150 (MQ=255)
                                                                                 gcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACCCCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCCACTGCGGGTTGTCg                                                                      <  2:3644591/150‑1 (MQ=255)
                                                                                  cgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGc                                                                     <  1:7276598/150‑1 (MQ=255)
                                                                                   gATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACTCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAGGACGATGACGCGCACCAACTGCGGGTTGTCgat                                                                    <  2:3149260/150‑3 (MQ=255)
                                                                                    aTCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                                   >  1:5482402/1‑150 (MQ=255)
                                                                                     tCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                                  <  1:1526588/150‑1 (MQ=255)
                                                                                     tCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                                  <  2:2495379/150‑1 (MQ=255)
                                                                                        ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                               >  2:6753982/1‑150 (MQ=255)
                                                                                        ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                               >  1:4855349/1‑150 (MQ=255)
                                                                                          gcggcgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatcatc                                                             >  1:900748/1‑150 (MQ=255)
                                                                                          gcggcgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatcatc                                                             >  1:2401983/1‑150 (MQ=255)
                                                                                                aTCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGcc                                                       >  1:4962898/1‑150 (MQ=255)
                                                                                                    ggCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                                   >  1:763032/1‑150 (MQ=255)
                                                                                                      ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATgg                                                 <  1:5145758/150‑1 (MQ=255)
                                                                                                      ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATgg                                                 <  1:6748434/150‑1 (MQ=255)
                                                                                                       cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                                >  2:2013486/1‑150 (MQ=255)
                                                                                                       cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                                <  2:5186776/150‑1 (MQ=255)
                                                                                                       cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                                >  1:1709465/1‑150 (MQ=255)
                                                                                                        aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCt                                               >  2:5770857/1‑150 (MQ=255)
                                                                                                             gCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATcc                                          <  1:4266168/150‑1 (MQ=255)
                                                                                                               gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                        >  1:1025708/1‑150 (MQ=255)
                                                                                                                  gCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGc                                     >  1:4755952/1‑150 (MQ=255)
                                                                                                                  gCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGc                                     >  2:5922850/1‑150 (MQ=255)
                                                                                                                     gcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGc                                  >  2:1811258/1‑150 (MQ=255)
                                                                                                                     gcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGc                                  >  2:636795/1‑150 (MQ=255)
                                                                                                                      cgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCa                                 <  1:56914/150‑1 (MQ=255)
                                                                                                                      cgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCa                                 <  2:7019995/150‑1 (MQ=255)
                                                                                                                       gAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAg                                <  1:2505738/150‑1 (MQ=255)
                                                                                                                       gAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAg                                >  2:2543121/1‑150 (MQ=255)
                                                                                                                       agtggCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAg                                >  2:2235398/4‑150 (MQ=255)
                                                                                                                          ggCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTg                             <  1:2378254/150‑1 (MQ=255)
                                                                                                                            ccACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCt                           >  1:1062211/1‑150 (MQ=255)
                                                                                                                            ccACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCt                           >  2:4605834/1‑150 (MQ=255)
                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGTGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  1:6283135/150‑1 (MQ=255)
                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  1:2937132/150‑1 (MQ=255)
                                                                                                                             cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                          <  1:5943281/150‑1 (MQ=255)
                                                                                                                              acacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCg                         >  1:3620657/1‑150 (MQ=255)
                                                                                                                               cacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                        <  2:1967658/150‑1 (MQ=255)
                                                                                                                               cacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                        >  1:5841666/1‑150 (MQ=255)
                                                                                                                               cacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                        <  2:5939406/150‑1 (MQ=255)
                                                                                                                                acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAg                       >  1:124257/1‑150 (MQ=255)
                                                                                                                                  cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCa                     >  1:4064924/1‑150 (MQ=255)
                                                                                                                                    aGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCAcc                   <  2:440308/150‑1 (MQ=255)
                                                                                                                                    aGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCAcc                   <  2:3932645/150‑1 (MQ=255)
                                                                                                                                      cACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGc                 <  2:6052237/150‑1 (MQ=255)
                                                                                                                                      cACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGc                 >  1:1665589/1‑150 (MQ=255)
                                                                                                                                         ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAtt              <  1:5230941/150‑1 (MQ=255)
                                                                                                                                         ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAtt              <  2:86360/150‑1 (MQ=255)
                                                                                                                                             gTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATg          >  2:3056114/1‑150 (MQ=255)
                                                                                                                                                gCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg       >  2:1642195/1‑150 (MQ=255)
                                                                                                                                                gCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg       <  1:1373945/150‑1 (MQ=255)
                                                                                                                                                  cACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCggg     <  1:6640952/150‑1 (MQ=255)
                                                                                                                                                  cACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCggg     <  1:6669496/150‑1 (MQ=255)
                                                                                                                                                   aCTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgc    <  2:2952808/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAG  >  NC_002947/4394352‑4394650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: