Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 590,711 T→A 6.5% W279R (TGG→AGG)  PP_0500 → dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947590,7110TA 6.5% 231.9 / 3.5 76W279R (TGG→AGG) PP_0500dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTCGA  >  NC_002947/590569‑590854
                                                                                                                                               |                                                                                                                                               
ggccggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGgcgcg                                                                                                                                           >  5:251482/1‑149 (MQ=255)
   cggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTg                                                                                                                                        >  7:196804/1‑149 (MQ=255)
     gccAATCCCGCCACCGGGCCGTGCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGACGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGgcacgcggc                                                                                                                                       <  6:33307/149‑4 (MQ=255)
       cAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCtt                                                                                                                                    <  6:251482/149‑1 (MQ=255)
         aTACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGc                                                                                                                                  <  8:196804/149‑1 (MQ=255)
           aCCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCa                                                                                                                                >  5:336375/1‑149 (MQ=255)
            ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGGCACGC‑GGCTTgatac                                                                                                                               <  5:61158/149‑5 (MQ=255)
            ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                               >  5:164647/1‑149 (MQ=255)
            ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                               <  2:117017/149‑1 (MQ=255)
            ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                               >  8:65066/1‑149 (MQ=255)
            ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                               <  6:311394/149‑1 (MQ=255)
             cGCAAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCTGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTg                                                                                                                                        <  2:200397/139‑1 (MQ=255)
             cGCAAGCTGGCCGTTCAGGCCCCCACGACGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTg                                                                                                                                        >  1:200397/1‑139 (MQ=255)
               ccAGCTGGCCGTTCAGG‑CCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACctc                                                                                                                            >  3:328069/1‑148 (MQ=255)
                 aGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTa                                                                                                                          <  7:65066/149‑1 (MQ=255)
                    tGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                         <  7:184150/147‑1 (MQ=255)
                    tGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                         >  8:184150/1‑147 (MQ=255)
                       ccGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAc                                                                                                                    >  6:112802/1‑149 (MQ=255)
                        cGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAcc                                                                                                                   >  7:25157/1‑149 (MQ=255)
                               gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGgcacgc                                                                                                                                          >  5:304731/1‑120 (MQ=255)
                               gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGgcacgc                                                                                                                                          <  6:304731/120‑1 (MQ=255)
                               gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGCCTTGCAACCTCTACTGGACCGGTTCTa                                                                                                            >  6:136396/1‑149 (MQ=255)
                                    cACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGcc                                                                                                       <  8:25157/149‑1 (MQ=255)
                                      cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCAt                                                                                                     <  8:277413/149‑1 (MQ=255)
                                       gccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGTCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATg                                                                                                    >  2:64377/1‑149 (MQ=255)
                                         cgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGc                                                                                                  >  7:28137/1‑149 (MQ=255)
                                          gcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTg                                                                                                                                                >  1:15394/1‑102 (MQ=255)
                                          gcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTg                                                                                                                                                <  2:15394/102‑1 (MQ=255)
                                             ggCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGAc                                                                                              <  6:164647/149‑1 (MQ=255)
                                                ccATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGcc                                                                                           >  7:227158/1‑149 (MQ=255)
                                                 cATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATgg                                                                                                   <  6:221884/140‑1 (MQ=255)
                                                 cATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATgg                                                                                                   >  5:221884/1‑140 (MQ=255)
                                                   tGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGTTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                         >  4:77247/1‑116 (MQ=255)
                                                   tGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGTTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                         <  3:77247/116‑1 (MQ=255)
                                                    gCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACAGCCATGGCTGACGccccc                                                                                        >  6:192588/1‑148 (MQ=255)
                                                                gATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACcg                                                                           >  2:65152/1‑149 (MQ=255)
                                                                    ccAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACcggcg                                                                        <  1:119398/148‑1 (MQ=255)
                                                                     cAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGt                                                                      <  6:49803/149‑1 (MQ=255)
                                                                     cAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACAGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGt                                                                      >  2:304980/1‑149 (MQ=255)
                                                                     cAGAGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGt                                                                      >  4:298393/1‑149 (MQ=255)
                                                                       gCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGc                                                                    >  6:257611/1‑149 (MQ=255)
                                                                       gCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGc                                                                    <  7:36708/149‑1 (MQ=255)
                                                                        cGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGcc                                                                   >  5:108575/1‑149 (MQ=255)
                                                                           agagCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGc                                                                <  5:136396/149‑1 (MQ=255)
                                                                             agCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCtt                                                              >  4:200370/1‑149 (MQ=255)
                                                                                 acaGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                          >  6:14369/1‑149 (MQ=255)
                                                                                  caGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTg                                                         <  5:192588/149‑1 (MQ=255)
                                                                                      aCGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGccccc                                                                                        >  1:128773/1‑114 (MQ=255)
                                                                                      aCGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGccccc                                                                                        <  2:128773/114‑1 (MQ=255)
                                                                                      aCGCGGTACTCGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGCCGCCCCAATCCTGATCACCGGCGGTGCCGGCTTCATTGGcac                                                     >  2:321050/1‑147 (MQ=255)
                                                                                                 gCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCtgtg                                            <  8:227158/147‑1 (MQ=255)
                                                                                                  ccTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGGGCGCTGGCTTGCCACCTCTACTGGTCCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGa                                         <  4:328069/149‑1 (MQ=255)
                                                                                                      cAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCg                                     <  5:112802/149‑1 (MQ=255)
                                                                                                         gAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCgttg                                  >  7:123458/1‑149 (MQ=255)
                                                                                                            aaTCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGttgttg                               <  5:14369/149‑1 (MQ=255)
                                                                                                            aaTCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGttgttg                               <  1:65152/149‑1 (MQ=255)
                                                                                                             aTCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGttgttg                               >  1:226683/1‑148 (MQ=255)
                                                                                                                 ttCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACa                           >  8:86125/1‑148 (MQ=255)
                                                                                                                 ttCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACa                           >  5:5168/1‑148 (MQ=255)
                                                                                                                   cATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                         <  6:336375/148‑1 (MQ=255)
                                                                                                                        cTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTAc                   <  3:298393/149‑1 (MQ=255)
                                                                                                                         ttCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATg                                       <  4:305284/128‑1 (MQ=255)
                                                                                                                         ttCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATg                                       >  3:305284/1‑128 (MQ=255)
                                                                                                                         ttCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACg                  >  8:183223/1‑149 (MQ=255)
                                                                                                                            ggTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAGGGCTAc                   <  3:60860/145‑1 (MQ=255)
                                                                                                                            ggTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAGGGCTAc                   >  4:60860/1‑145 (MQ=255)
                                                                                                                            ggTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGcc                <  5:257611/148‑1 (MQ=255)
                                                                                                                             gTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCCTCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCg               >  7:3124/1‑148 (MQ=255)
                                                                                                                              tATGAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTa             <  7:187925/149‑1 (MQ=255)
                                                                                                                                    gccagcgcgcCAGGCACGC‑GGCTTGATACCACTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCAtt       >  1:203247/5‑149 (MQ=255)
                                                                                                                                    gCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCAtt       >  6:33688/1‑149 (MQ=255)
                                                                                                                                     cc‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCAc                                                 >  5:300967/1‑106 (MQ=255)
                                                                                                                                     cc‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCAc                                                 <  6:300967/106‑1 (MQ=255)
                                                                                                                                     cc‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATtc      <  1:136216/149‑1 (MQ=255)
                                                                                                                                          gTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTcga  <  7:183223/149‑1 (MQ=255)
                                                                                                                                          gTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTcga  <  6:5168/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATTCTCGA  >  NC_002947/590569‑590854

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