Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,305,604 G→A 5.8% R151C (CGC→TGC)  PP_3777 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,305,6040GA 5.8% 260.1 / 3.2 86R151C (CGC→TGC) PP_3777hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/47);  new base A (3/2);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305458‑4305745
                                                                                                                                                  |                                                                                                                                             
acCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCt                                                                                                                                             >  2:217955/1‑149 (MQ=255)
  cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCGGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTg                                                                                                                                            >  1:14952/1‑148 (MQ=255)
   ggCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                          <  1:226445/149‑1 (MQ=255)
    gCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAg                                                                                                                                         >  4:243012/1‑149 (MQ=255)
    gCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAg                                                                                                                                         >  2:208030/1‑149 (MQ=255)
     cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  5:177549/149‑1 (MQ=255)
     cAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGt                                                                                                                                        <  6:254388/149‑1 (MQ=255)
        tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     >  5:158067/1‑149 (MQ=255)
        tGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:122668/149‑1 (MQ=255)
          ccGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATc                                                                                                                                   >  4:110458/1‑149 (MQ=255)
            gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGa                                                                                                                                 <  1:208030/149‑1 (MQ=255)
            gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGa                                                                                                                                 <  8:36439/149‑1 (MQ=255)
            gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 <  3:110458/149‑1 (MQ=255)
                aaCGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGAGTCCCGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGCggca                                                                                                                             >  3:156154/1‑147 (MQ=255)
                      ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     >  1:141429/1‑135 (MQ=255)
                      ttCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCa                                                                                                                                     <  2:141429/135‑1 (MQ=255)
                          gCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAg                                                                                                                   >  7:181066/1‑149 (MQ=255)
                           cccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGTTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTTTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                  <  8:227887/149‑1 (MQ=255)
                             cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  2:226422/1‑149 (MQ=255)
                             cGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                >  7:114967/1‑149 (MQ=255)
                                  atCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCt                                                                                                            <  3:248400/148‑1 (MQ=255)
                                   tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               <  5:206334/128‑1 (MQ=255)
                                   tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAgg                                                                                                                               >  6:206334/1‑128 (MQ=255)
                                   tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGaa                                                                                                                             <  4:206581/130‑1 (MQ=255)
                                   tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGaa                                                                                                                             >  3:206581/1‑130 (MQ=255)
                                   tCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGc                                                                                                          <  3:243012/149‑1 (MQ=255)
                                     aGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCa                                                                                                        >  5:179031/1‑149 (MQ=255)
                                      gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAc                                                                                                       >  4:108657/1‑149 (MQ=255)
                                       cTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCAcc                                                                                                      >  3:51380/1‑149 (MQ=255)
                                          gCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTg                                                                                                   >  1:238869/1‑149 (MQ=255)
                                            cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGTCAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 >  7:222715/1‑117 (MQ=255)
                                            cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGTCAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                 <  8:222715/117‑1 (MQ=255)
                                            cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  6:243160/149‑1 (MQ=255)
                                            cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  5:64318/149‑1 (MQ=255)
                                            cTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGc                                                                                                 <  3:108657/149‑1 (MQ=255)
                                             tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                >  4:76771/1‑149 (MQ=255)
                                             tCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                >  8:240028/1‑149 (MQ=255)
                                                gggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGACGTCTGCCACCCTGGCGGAt                                                                                             <  3:74607/149‑1 (MQ=255)
                                                 ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  3:215019/149‑1 (MQ=255)
                                                 ggCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  4:238290/149‑1 (MQ=255)
                                                    ttCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATACAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  2:62614/1‑149 (MQ=255)
                                                      cAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  4:51380/147‑1 (MQ=255)
                                                      cAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAAt                                                                                        <  6:179031/148‑1 (MQ=255)
                                                       aGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTNAATa                                                                                       <  5:3743/148‑1 (MQ=255)
                                                         gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  2:54119/1‑144 (MQ=255)
                                                         gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  1:54119/144‑1 (MQ=255)
                                                         gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  3:205536/149‑1 (MQ=255)
                                                         gTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  1:226422/149‑1 (MQ=255)
                                                         gTAGCGCGGCACCCGTTCGAAACACCGGCGCACGGTCCGCGCGCCATATAACAGCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGa                                                                                    <  4:156154/149‑1 (MQ=255)
                                                             cTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATca                                                                                 >  2:112429/1‑148 (MQ=255)
                                                                ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATc                                                                                                                                   >  4:22279/1‑95 (MQ=255)
                                                                ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATc                                                                                                                                   <  3:22279/95‑1 (MQ=255)
                                                                ggCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGAGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATcaca                                                                               <  2:148622/147‑1 (MQ=255)
                                                                  cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            <  1:61490/132‑1 (MQ=255)
                                                                  cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                            >  2:61490/1‑132 (MQ=255)
                                                                  cACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCg                                                                            >  5:230136/1‑148 (MQ=255)
                                                                      cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTa                                                                          >  8:13777/1‑146 (MQ=255)
                                                                      cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTa                                                                          <  7:13777/146‑1 (MQ=255)
                                                                            aaaaaCCTGCGCAAGGTCGGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCctc                                                                 <  2:14952/149‑1 (MQ=255)
                                                                                        aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAgg                                                                                     >  8:178548/1‑117 (MQ=255)
                                                                                        aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAgg                                                                                     <  7:178548/117‑1 (MQ=255)
                                                                                        aaGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAt                                                     >  7:85237/1‑149 (MQ=255)
                                                                                         aGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  8:114967/149‑1 (MQ=255)
                                                                                         aGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  3:76771/149‑1 (MQ=255)
                                                                                         aGGTCCGCTCGCCATACAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  1:62614/149‑1 (MQ=255)
                                                                                            tCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCg                                                 >  6:58637/1‑149 (MQ=255)
                                                                                               gatagCCATAGAACATCGCTTACAGCGACAGCCCTTGGGCAGCCAAGATGCCCTGTAGATCATCGAGGAACCCATCCGAGCGCGTCTGCCACCCTCGCGTATTCAATAGGATCAAATCGAATACGACCTCTGCCCAAGGCATAGccag                                               >  4:82074/5‑145 (MQ=255)
                                                                                               gCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGc                                                  <  3:82074/145‑1 (MQ=255)
                                                                                                   gCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    <  8:69860/139‑1 (MQ=255)
                                                                                                   gCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCAtt                                                    >  7:69860/1‑139 (MQ=255)
                                                                                                     cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                             >  3:106454/1‑112 (MQ=255)
                                                                                                     cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATc                                                                             <  4:106454/112‑1 (MQ=255)
                                                                                                           aCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTccc                                  <  6:158067/149‑1 (MQ=255)
                                                                                                            cATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCt                                 >  3:59260/1‑149 (MQ=255)
                                                                                                            cATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCt                                 <  7:234851/149‑1 (MQ=255)
                                                                                                                      aGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTc                       >  5:228510/1‑149 (MQ=255)
                                                                                                                      aGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTc                       <  8:85237/149‑1 (MQ=255)
                                                                                                                       gCGACAGCCATTGGTCAGCCTCGATGCGCTTTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCg                      <  2:238869/149‑1 (MQ=255)
                                                                                                                        cGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGa                     <  7:240028/149‑1 (MQ=255)
                                                                                                                             gCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAg                <  4:30496/149‑1 (MQ=255)
                                                                                                                                   tggCAGCCTCGTTGCGCTGTTGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGATTCCGTCCTCTGCGGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTc           <  5:58637/147‑1 (MQ=255)
                                                                                                                                   gggCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGAtt            <  4:59260/147‑1 (MQ=255)
                                                                                                                                        tccTCGATGCACTACAGCGCATCGAGGCAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAg     >  3:137409/2‑149 (MQ=255)
                                                                                                                                          cTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt   <  6:228510/149‑1 (MQ=255)
                                                                                                                                          cTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt   <  1:217955/149‑1 (MQ=255)
                                                                                                                                           tCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGtt  >  4:94440/1‑149 (MQ=255)
                                                                                                                                                  gcTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         >  3:248503/1‑55 (MQ=255)
                                                                                                                                                  gcTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                         <  4:248503/55‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTT  >  NC_002947/4305458‑4305745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: