Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 563,007 T→G 5.7% P30P (CCT→CCG) ‡ rpoA → DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947563,0070TG 5.7% 222.6 / 4.1 70P30P (CCT→CCG) ‡rpoADNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base G (2/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

AACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAATT  >  NC_002947/562860‑563154
                                                                                                                                                    |                                                                                                                                                   
aaCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCct                                                                                                                                                     <  1:275073/148‑1 (MQ=255)
 aCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctct                                                                                                                                                   >  1:288515/1‑149 (MQ=255)
    aaaaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGa                                                                                                                                                <  6:203492/149‑1 (MQ=255)
      aaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGc                                                                                                                                              <  7:86147/149‑1 (MQ=255)
         cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  6:359257/149‑1 (MQ=255)
         cTGATTGGCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGATCCTGACTCCCCGCCACAT‑TGATCTGCAGGAAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  7:295977/149‑1 (MQ=255)
            aTTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGCGtaatt                                                                                                                                        <  7:289014/149‑5 (MQ=255)
             ttGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGttt                                                                                                                                       <  6:101634/149‑1 (MQ=255)
                tCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAG‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTCGTTTCgg                                                                                                                                    >  6:86700/1‑149 (MQ=255)
                  gAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                   <  5:266003/148‑1 (MQ=255)
                      tctACTCCAAGTAAGGTCTAGAAAATAGGTGCATCCATGCAGATTTCTGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  1:123147/136‑1 (MQ=255)
                      tctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAc                                                                                                                              >  3:232876/1‑149 (MQ=255)
                      tctACTACAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTTCAGGTAGTCAGTCCGACAAGCGCCAAAATTACGATAGAGCCTCTCGAGCGTg                                                                                                                                           >  2:123147/1‑136 (MQ=255)
                          cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTg                                                                                                                          <  5:387292/149‑1 (MQ=255)
                            ccAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTggg                                                                                                                        >  3:243046/1‑149 (MQ=255)
                             cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                        <  6:202687/132‑1 (MQ=255)
                             cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                        >  5:202687/1‑132 (MQ=255)
                                      gCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                           >  7:101821/1‑136 (MQ=255)
                                      gCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                           <  8:101821/136‑1 (MQ=255)
                                        tAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTgcg                                                                                                            >  2:333536/1‑149 (MQ=255)
                                                  gTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCa                                                                                                                                 >  2:50960/1‑118 (MQ=255)
                                                  gTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCa                                                                                                                                 <  1:50960/118‑1 (MQ=255)
                                                     cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCa                                                                                                       <  2:142135/141‑1 (MQ=255)
                                                     cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCa                                                                                                       >  1:142135/1‑141 (MQ=255)
                                                     cATCCATGCAGATTTCGGTAAATGAGGTCCTGACTCCCCGACACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCTCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                                <  8:153405/148‑1 (MQ=255)
                                                      aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATGTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCtctc                                                                                                                                                  >  1:78877/1‑98 (MQ=255)
                                                      aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCtctc                                                                                                                                                  <  2:78877/97‑1 (MQ=255)
                                                      aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                              >  6:42261/1‑149 (MQ=255)
                                                           tGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcctc                                                                                          >  6:178064/1‑148 (MQ=255)
                                                             cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttg                                                                                               >  7:257405/1‑141 (MQ=255)
                                                             cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttg                                                                                               <  8:257405/141‑1 (MQ=255)
                                                             cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCAt                                                                                       <  2:158959/149‑1 (MQ=255)
                                                              aGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATg                                                                                      <  5:365998/149‑1 (MQ=255)
                                                                 tttCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                                   <  4:232876/149‑1 (MQ=255)
                                                                     ggTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCt                                                                               >  5:35981/1‑149 (MQ=255)
                                                                         aaTGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGc                                                                           <  2:288515/149‑1 (MQ=255)
                                                                                 ccTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                    >  2:378572/1‑148 (MQ=255)
                                                                                  cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                    >  6:126757/1‑147 (MQ=255)
                                                                                  cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                    >  2:278407/1‑147 (MQ=255)
                                                                                  cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                  <  6:35981/149‑1 (MQ=255)
                                                                                    gACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                                >  2:300956/1‑149 (MQ=255)
                                                                                     aCTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGc                                                               <  1:300956/149‑1 (MQ=255)
                                                                                       tCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                              >  1:278330/1‑148 (MQ=255)
                                                                                       tCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                              >  7:180426/1‑148 (MQ=255)
                                                                                         cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                                >  6:205978/1‑144 (MQ=255)
                                                                                         cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                                <  5:205978/144‑1 (MQ=255)
                                                                                         cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                              >  3:370249/1‑146 (MQ=255)
                                                                                         cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCga                                                            <  2:59201/148‑1 (MQ=255)
                                                                                          ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           >  5:916/1‑68 (MQ=255)
                                                                                          ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                           <  6:916/68‑1 (MQ=255)
                                                                                          ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCgaga                                                          <  4:370249/149‑1 (MQ=255)
                                                                                           agCCAGAT‑TGATGTGCAGGTTGTCAGTCCGACCGGGTCCAAAATTACGCTCGAGCCTCTCGAGCGTGGATTCGGCCTTACCCTGGGCAACGCGCTGCGCCGCATCCTGGTGTCCTGCATGCCTGGGTGTGCAGTAGTCGAGGCCGAGAt                                                         <  5:222323/148‑1 (MQ=255)
                                                                                             ccACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCg                                                       <  1:378572/149‑1 (MQ=255)
                                                                                              cacaT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGa                                                      <  8:180426/149‑1 (MQ=255)
                                                                                                 aT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATgg                                                   >  8:222013/1‑149 (MQ=255)
                                                                                                  t‑tGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTg                                                                                                              <  5:231095/89‑1 (MQ=255)
                                                                                                  t‑tGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTg                                                                                                              >  6:231095/1‑89 (MQ=255)
                                                                                                         tgCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACt                                             <  4:243046/148‑1 (MQ=255)
                                                                                                            aGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCAt                                         <  5:42261/149‑1 (MQ=255)
                                                                                                            aGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCAt                                         <  5:86700/149‑1 (MQ=255)
                                                                                                                   cagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCTGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                  <  8:336499/117‑1 (MQ=255)
                                                                                                                   cagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCTGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                  >  7:336499/1‑117 (MQ=255)
                                                                                                                     gtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTcc                                >  3:219669/1‑149 (MQ=255)
                                                                                                                      tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCAACCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcg                               >  1:121165/1‑149 (MQ=255)
                                                                                                                           cccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATc                          <  1:333536/149‑1 (MQ=255)
                                                                                                                            ccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATc                          >  5:19778/1‑148 (MQ=255)
                                                                                                                                gcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAgg                     >  5:22129/1‑149 (MQ=255)
                                                                                                                                    aaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGtt                 <  6:192652/149‑1 (MQ=255)
                                                                                                                                        caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg              >  7:171114/2‑148 (MQ=255)
                                                                                                                                         aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAgg             >  2:175829/1‑148 (MQ=255)
                                                                                                                                             tCGAGCCTCTCGAGCGTGGTTTCGGCCAGACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGAc        >  3:216491/1‑149 (MQ=255)
                                                                                                                                               gAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGAt                           <  3:233710/128‑1 (MQ=255)
                                                                                                                                               gAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGAt                           >  4:233710/1‑128 (MQ=255)
                                                                                                                                                 gCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTaa    <  2:121165/149‑1 (MQ=255)
                                                                                                                                                    tctcGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAAtt  >  3:2565/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAATT  >  NC_002947/562860‑563154

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