Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,063,976 G→C 7.4% L87L (CTG→CTC)  PP_2675 → cytochrome c‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,063,9760GC 7.4% 205.6 / 4.8 64L87L (CTG→CTC) PP_2675cytochrome c‑type protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base C (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGACCCCGATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTACCTCAGCCA  >  NC_002947/3063828‑3064122
                                                                                                                                                    |                                                                                                                                                  
gCCTGACCCCGATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTg                                                                                                                                                    <  7:229866/149‑1 (MQ=255)
     aCCCCGATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGc                                                                                                                                               <  8:36521/149‑1 (MQ=255)
       cccGATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTCCAGGCCa                                                                                                                                             >  3:97848/1‑149 (MQ=255)
          gATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTCCAGGCCatga                                                                                                                                          <  4:97848/149‑4 (MQ=255)
            tCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTCCAGGCCatgaca                                                                                                                                        <  7:69414/149‑6 (MQ=255)
                    acTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAgg                                                                                                                                                >  7:174559/1‑133 (MQ=255)
                    acTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAgg                                                                                                                                                <  8:174559/133‑1 (MQ=255)
                    acTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTAt                                                                                                                                >  4:72216/1‑149 (MQ=255)
                      tGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATcg                                                                                                                              >  4:286413/1‑149 (MQ=255)
                        gCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATcgcg                                                                                                                            >  3:93167/1‑149 (MQ=255)
                          gTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGcc                                                                                                                          >  6:240311/1‑149 (MQ=255)
                              gCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGAc                                                                                                                      <  7:254603/149‑1 (MQ=255)
                                cTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACtt                                                                                                                    >  2:242912/1‑149 (MQ=255)
                                 tGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTg                                                                                                                   >  2:301787/1‑149 (MQ=255)
                                           cggcTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCg                                                                                                         <  1:242912/149‑1 (MQ=255)
                                             gcTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGAt                                                                                                       >  1:97094/1‑149 (MQ=255)
                                               tGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGc                                                                                                     >  6:287788/1‑149 (MQ=255)
                                                 ggCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGc                                                                                                     <  5:240311/147‑1 (MQ=255)
                                                  gCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCggg                                                                                                  <  4:93167/149‑1 (MQ=255)
                                                   ccGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGt                                                                                                 >  7:269992/1‑149 (MQ=255)
                                                        gTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGc                                                                                            >  6:156075/1‑149 (MQ=255)
                                                         tCAACCCCTACCGCAGCTCGCCCGAGCACGACAAGCCCGTGGAAATCGGCGCCTCCGCCTACAACCAGAACTGCGCTGCCTGTCATCACCTCGAGGCCAAGTCAGCTGGTATCGCGACTGAAATGCTCATGCTCCAAGCGGGAGATGcc                                                                                           >  4:234695/1‑149 (MQ=255)
                                                             cccGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGc                                                                                        >  7:119466/1‑148 (MQ=255)
                                                             cccGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGc                                                                                        >  7:251211/1‑148 (MQ=255)
                                                               cGTACCGCAGCTCGCCCTAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGAt                                                                                     <  7:225987/149‑1 (MQ=255)
                                                               cGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGAt                                                                                     <  3:234695/149‑1 (MQ=255)
                                                                gTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGAt                                                                                     >  3:91229/1‑148 (MQ=255)
                                                                   ccGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAAt                                                                                 <  1:12536/149‑1 (MQ=255)
                                                                    cGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATg                                                                                >  8:239437/1‑149 (MQ=255)
                                                                     gCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATgg                                                                               >  4:218854/1‑149 (MQ=255)
                                                                      cAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGt                                                                              >  5:224888/1‑149 (MQ=255)
                                                                        gCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTtc                                                                            <  3:15520/149‑1 (MQ=255)
                                                                         cTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGc                                                                                            <  7:71615/132‑1 (MQ=255)
                                                                         cTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGc                                                                                            >  8:71615/1‑132 (MQ=255)
                                                                         cTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTtcg                                                                           >  7:84035/1‑149 (MQ=255)
                                                                          tCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTtcgt                                                                          >  6:79226/1‑149 (MQ=255)
                                                                               cGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGc                                                                     <  5:156075/149‑1 (MQ=255)
                                                                                gAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCg                                                                    >  3:164952/1‑149 (MQ=255)
                                                                                 aGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCgg                                                                   <  8:251211/149‑1 (MQ=255)
                                                                                  gCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCggg                                                                  <  4:164952/149‑1 (MQ=255)
                                                                                    aCGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTa                                                                <  3:93494/149‑1 (MQ=255)
                                                                                      gataaGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGt                                                                                                 >  8:256393/4‑114 (MQ=255)
                                                                                      gataaGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGt                                                                                                 <  7:256393/111‑1 (MQ=255)
                                                                                      gACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACg                                                              >  3:86803/1‑149 (MQ=255)
                                                                                            gCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACgg                                                        >  6:199237/1‑149 (MQ=255)
                                                                                                tGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGcgcg                                                    <  4:91229/149‑1 (MQ=255)
                                                                                                    aaTCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTc                                                <  5:287788/149‑1 (MQ=255)
                                                                                                      tCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTcgc                                              <  8:84035/149‑1 (MQ=255)
                                                                                                         gcgcCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGAc                                           <  3:218854/149‑1 (MQ=255)
                                                                                                          cgcCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACg                                          <  5:79226/149‑1 (MQ=255)
                                                                                                               cGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCcgc                                     <  5:199237/149‑1 (MQ=255)
                                                                                                                 gCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCgt                                   <  8:119466/149‑1 (MQ=255)
                                                                                                                      cAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTa                                <  1:301787/147‑1 (MQ=255)
                                                                                                                       aaCCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACAt                             <  3:146578/149‑1 (MQ=255)
                                                                                                                        aCCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATg                            <  3:286413/149‑1 (MQ=255)
                                                                                                                                tGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACg                                          >  7:263984/1‑127 (MQ=255)
                                                                                                                                tGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACg                                          <  8:263984/127‑1 (MQ=255)
                                                                                                                                 gcgcTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCGACGGCGCGGTTCGCGACGGCCGCGTGTACATGc                           >  5:280718/1‑141 (MQ=255)
                                                                                                                                 gcgcTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGc                           <  6:280718/141‑1 (MQ=255)
                                                                                                                                  cgcTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATgg                  >  7:251225/1‑149 (MQ=255)
                                                                                                                                           tCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACTGCCGCGTGTACATGCCGAAGATGGCCGATTa           <  6:224888/147‑1 (MQ=255)
                                                                                                                                           tCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTa           <  8:269992/147‑1 (MQ=255)
                                                                                                                                           tCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTAcc         <  3:72216/149‑1 (MQ=255)
                                                                                                                                           acttGGCCTCCAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTAcc         >  7:251573/4‑149 (MQ=255)
                                                                                                                                                  cTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTACCTCAGCCa  <  2:97094/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCTGACCCCGATCAAGGACACTGGCGTGAGCCTGGACGGCGACGGCTGGGCCGCGGTCAACCCGTACCGCAGCTCGCCCGAGCACGACAAGGCCGTGGAAATCGGCGCCTCGGCCTACAACCAGAACTGCGCTGCCTGTCATGGCCTGGAGGCCAAGTCCGGTGGTATCGCGCCTGACTTGCGCATGCTCGATGCGGGTGATGCCGGCGATGAATGGTTCGTCGAGCGGGTACGCCACGGCGCGGTTCGCGACGGCCGCGTGTACATGCCGAAGATGGCCGATTACCTCAGCCA  >  NC_002947/3063828‑3064122

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