Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 961,720 A→T 5.8% L232Q (CTG→CAG)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947961,7200AT 5.8% 252.2 / 3.7 85L232Q (CTG→CAG) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base T (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAG  >  NC_002947/961572‑961865
                                                                                                                                                     |                                                                                                                                                 
cccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCa                                                                                                                                                   >  4:298020/1‑149 (MQ=255)
   gTAGCGGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACTGCGGCGTAGTTCAGGTGCGCGTTGAGGCCGCTGCCTTCGGCAAGCTCGT‑CCAGGTGTCCGAGCAGGc                                                                                                                                                <  6:294904/149‑1 (MQ=255)
    tAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGG‑CCACCTGGTCGAGCTGGcc                                                                                                                                               <  5:164204/149‑1 (MQ=255)
        tGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAg                                                                                                                                           >  4:254480/1‑149 (MQ=255)
          tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGCTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGCGCAGGCCGAAGcc                                                                                                                                         >  8:165371/1‑149 (MQ=255)
            gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCgg                                                                                                                                       <  1:146738/149‑1 (MQ=255)
             tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                      <  6:179493/149‑1 (MQ=255)
             tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGCTAATGGTCGACACTAGGCAGGGGTGGCTCGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCTGCTCGG‑CCACCTGGTCGAGCTGGCCGAAggcagc                                                                                                                                      <  4:277654/149‑6 (MQ=37)
                   cgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                                <  2:171823/149‑1 (MQ=255)
                   cgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                                >  3:20314/1‑149 (MQ=255)
                    gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                               >  1:103527/1‑149 (MQ=255)
                    gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                               >  6:148549/1‑149 (MQ=255)
                    gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                               >  7:78404/1‑149 (MQ=255)
                     cAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                              <  4:224109/149‑1 (MQ=255)
                      aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                             >  6:371626/1‑149 (MQ=255)
                         gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCAt                                                                                                                          <  8:138257/149‑1 (MQ=255)
                          ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                         <  3:298020/149‑1 (MQ=255)
                              aCCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCt                                                                                                                     >  6:84481/1‑149 (MQ=255)
                               ccGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                                    >  3:316137/1‑149 (MQ=255)
                                         aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                           >  5:26345/1‑148 (MQ=255)
                                          gCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                         <  2:150549/149‑1 (MQ=255)
                                           ccTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAg                                                                                                        >  1:377156/1‑149 (MQ=255)
                                            cTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGt                                                                                                                                 >  4:298016/1‑123 (MQ=255)
                                            cTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGt                                                                                                                                 <  3:298016/123‑1 (MQ=255)
                                               cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                                     <  8:78404/148‑1 (MQ=255)
                                               cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGtt                                                                                                      <  2:103527/147‑1 (MQ=255)
                                               cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGtt                                                                                                      <  7:85552/147‑1 (MQ=255)
                                               cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTAGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGtt                                                                                                      <  1:104957/147‑1 (MQ=255)
                                                  ccAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGcc                                                                                                  >  2:386931/1‑148 (MQ=255)
                                                         aTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                           <  8:54696/148‑1 (MQ=255)
                                                             tCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCCGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                                      >  6:12492/1‑149 (MQ=255)
                                                                 cacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAACAAGA‑CCAGGTTGCCGACCAGGCCCCAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTTGCGAAGAGGGTGCCCGGc                                                                                  >  2:204719/1‑149 (MQ=255)
                                                                    tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTGGGCCACCACGG‑CCAGGTGGCCGCGCAGGCCGAAGGCGGTTGCGTCGGTCATGTGCTTTTCACCCGGCAGGTTGCCGAAGCGCCTTCcgtgcgtg                                                                               >  6:236065/1‑149 (MQ=255)
                                                                          ggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                         <  7:272858/149‑1 (MQ=255)
                                                                            cGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                       <  1:386931/149‑1 (MQ=255)
                                                                             ggggCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGTGCGTGTTGAGGGTg                                                                      <  1:404771/147‑1 (MQ=255)
                                                                              tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGc                                                                     <  5:371626/149‑1 (MQ=255)
                                                                                 cACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGcaca                                                                  >  3:60187/1‑149 (MQ=255)
                                                                                   cggcggCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATc                                                                <  4:20314/149‑1 (MQ=255)
                                                                                    ggcggcGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                                    >  5:83548/1‑96 (MQ=255)
                                                                                    ggcggcGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                                    <  6:83548/96‑1 (MQ=255)
                                                                                      cggcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAg                                                             >  2:225845/1‑149 (MQ=255)
                                                                                       ggcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGt                                                            >  7:251347/1‑149 (MQ=255)
                                                                                          gTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGAGCCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcc                                                                                                                >  1:83595/1‑95 (MQ=255)
                                                                                          gTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcc                                                                                                                <  2:83595/94‑1 (MQ=255)
                                                                                              ttCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                           <  5:357401/111‑1 (MQ=255)
                                                                                              ttCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                           >  6:357401/1‑111 (MQ=255)
                                                                                                cAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGcc                                                    <  1:161664/148‑1 (MQ=255)
                                                                                                cAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGc                                                     >  1:384222/1‑147 (MQ=255)
                                                                                                 aGGTGCGCGTTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGGTGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCa                                                   <  5:236065/148‑1 (MQ=255)
                                                                                                 aGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                                  >  2:249477/1‑149 (MQ=255)
                                                                                                   gTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGc                                                >  3:120167/1‑149 (MQ=255)
                                                                                                      cgcgGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCtt                                             <  6:26345/149‑1 (MQ=255)
                                                                                                        cgGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                           <  5:84481/149‑1 (MQ=255)
                                                                                                        cgGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                           <  4:60187/149‑1 (MQ=255)
                                                                                                          gTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTc                                         <  5:148549/149‑1 (MQ=255)
                                                                                                              ggCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGc                                     <  8:251347/149‑1 (MQ=255)
                                                                                                                 cGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                         <  1:155919/78‑1 (MQ=255)
                                                                                                                 cGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                         >  2:155919/1‑78 (MQ=255)
                                                                                                                   cTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACg                                <  1:204719/149‑1 (MQ=255)
                                                                                                                    tGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGc                               <  4:316137/149‑1 (MQ=255)
                                                                                                                       cTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGc                            >  7:105319/1‑149 (MQ=255)
                                                                                                                        ttCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTTCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCt                           >  7:9426/1‑149 (MQ=255)
                                                                                                                          cGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                          <  3:254480/148‑1 (MQ=255)
                                                                                                                           ggCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTgg                        >  5:115708/1‑149 (MQ=255)
                                                                                                                            gCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                       >  7:210246/1‑149 (MQ=255)
                                                                                                                             cctgCTCGG‑CCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCACGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCATGGTCATGCACACGCAGCTTGGcc                      >  3:313821/4‑149 (MQ=255)
                                                                                                                                  tCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                          <  6:51994/140‑1 (MQ=255)
                                                                                                                                  tCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                          >  5:51994/1‑140 (MQ=255)
                                                                                                                                   cGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCttt                >  1:394448/1‑149 (MQ=255)
                                                                                                                                          ggtggtCGAGCTGGCCGAAGGCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCtttt               >  6:272418/1‑144 (MQ=255)
                                                                                                                                          ggtggtCGAGCTGGCCGAAGGCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCtttt               <  5:272418/144‑1 (MQ=255)
                                                                                                                                           gTGGCCGAGCAGGACGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAg         >  4:147590/1‑149 (MQ=255)
                                                                                                                                            tGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgt        >  7:243496/1‑149 (MQ=255)
                                                                                                                                              gCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtg       <  1:225845/148‑1 (MQ=255)
                                                                                                                                               ccGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                            >  6:273136/1‑110 (MQ=255)
                                                                                                                                               ccGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                            <  5:273136/110‑1 (MQ=255)
                                                                                                                                               ccGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCGTTCTTTTCGGCAgtggt     >  4:47966/1‑149 (MQ=255)
                                                                                                                                                cGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtggtg    <  1:249477/149‑1 (MQ=255)
                                                                                                                                                cGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtggt     >  7:269931/1‑148 (MQ=255)
                                                                                                                                                cGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACAAGCAGATTGGCCTTATTTTCGGAAgtggtg    >  6:276729/1‑149 (MQ=255)
                                                                                                                                                 gAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGCAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtggt     >  6:101496/1‑147 (MQ=255)
                                                                                                                                                 gAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGa   <  5:12492/149‑1 (MQ=255)
                                                                                                                                                 gAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGa   <  6:115708/149‑1 (MQ=255)
                                                                                                                                                   gCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAg  >  1:406781/1‑148 (MQ=255)
                                                                                                                                                   gCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAg  >  3:65859/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGA‑CCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAG  >  NC_002947/961572‑961865

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