Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,202,857 G→T 6.2% R229L (CGG→CTG)  xcpY → type II secretion pathway protein XcpY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,202,8570GT 6.2% 202.8 / 3.5 64R229L (CGG→CTG) xcpYtype II secretion pathway protein XcpY
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base T (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

TGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCT  >  NC_002947/1202710‑1203001
                                                                                                                                                   |                                                                                                                                                
tGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGTCCGCCTCCCGCCAGCAGCAGCgg                                                                                                                                                 <  6:133687/149‑1 (MQ=255)
  ccGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACTCACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgct                                                                                                                                               >  1:114239/1‑149 (MQ=255)
  ccGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgct                                                                                                                                               >  5:248930/1‑149 (MQ=255)
      aGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGctgctg                                                                                                                                           <  3:219772/149‑1 (MQ=255)
       gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTgg                                                                                                                                          <  1:118286/149‑1 (MQ=255)
       gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTgg                                                                                                                                          <  2:114239/149‑1 (MQ=255)
        ccACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCgg                                                                                                                                                 <  7:135878/141‑1 (MQ=255)
        ccACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCgg                                                                                                                                                 >  8:135878/1‑141 (MQ=255)
              tGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACGGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGAGCCCGCAAGGTGAGCTCGTTGCCGGCCGCCTACCGCAAGCAGAAGCGGCTGCAGGCAGCATg                                                                                                                                   >  3:255959/1‑149 (MQ=255)
              tGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTg                                                                                                                                   <  5:407286/149‑1 (MQ=255)
               gggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGc                                                                                                                                  <  8:203398/149‑1 (MQ=255)
                ggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGc                                                                                                                                  <  4:5878/148‑1 (MQ=255)
                  cAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTg                                                                                                                               >  4:3447/1‑149 (MQ=255)
                     tgctgcACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTg                                                                                                                            >  8:27750/1‑149 (MQ=255)
                        tgcACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTg                                                                                                                         <  3:50771/149‑1 (MQ=255)
                        tgcACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTg                                                                                                                         >  3:368468/1‑149 (MQ=255)
                          cacaGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGc                                                                                                                       >  7:163043/1‑149 (MQ=255)
                          cacaGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGc                                                                                                                       >  6:132955/1‑149 (MQ=255)
                          cacaGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGc                                                                                                                       <  3:3447/149‑1 (MQ=255)
                            caGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTa                                                                                                                     >  3:88330/1‑149 (MQ=255)
                                   tCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTg                                                                                                              <  4:368468/149‑1 (MQ=255)
                                     gCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                             <  7:27750/148‑1 (MQ=255)
                                         aGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTg                                                                                                         <  6:248930/148‑1 (MQ=255)
                                          gCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGc                                                                                                       <  7:59369/149‑1 (MQ=255)
                                               aaTTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                                  <  8:163043/149‑1 (MQ=255)
                                               aaTTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                                  <  4:88330/149‑1 (MQ=255)
                                                    gCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                             <  5:132955/149‑1 (MQ=255)
                                                       aaCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                           >  1:86589/1‑148 (MQ=255)
                                                        aCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                           >  3:180265/1‑147 (MQ=255)
                                                            tATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                               >  8:319608/1‑139 (MQ=255)
                                                            tATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                               <  7:319608/139‑1 (MQ=255)
                                                                tGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGc                                                                                                       <  2:98598/127‑1 (MQ=255)
                                                                tGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGc                                                                                                       >  1:98598/1‑127 (MQ=255)
                                                                            cacacTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAAc                                                                     >  6:177277/1‑149 (MQ=255)
                                                                               acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCAGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAActct                                                                  >  2:309449/1‑149 (MQ=255)
                                                                                cTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                                 >  3:181905/1‑149 (MQ=255)
                                                                                  gCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                               >  8:209295/1‑149 (MQ=255)
                                                                                   cccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                               <  4:180265/148‑1 (MQ=255)
                                                                                       gCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATg                                                                                   >  1:149450/1‑124 (MQ=255)
                                                                                       gCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATg                                                                                   <  2:149450/124‑1 (MQ=255)
                                                                                       gCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTAAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGc                                                          >  1:301975/1‑149 (MQ=255)
                                                                                          tGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAg                                                       >  7:142166/1‑149 (MQ=255)
                                                                                            ttCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGTCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGt                                                     >  2:248810/1‑149 (MQ=255)
                                                                                            ttCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAg                                                       <  7:209295/147‑1 (MQ=255)
                                                                                              cGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCa                                                   >  3:93290/1‑149 (MQ=255)
                                                                                               gACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCAt                                                  >  2:255482/1‑149 (MQ=255)
                                                                                               gACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCAt                                                  >  2:283367/1‑149 (MQ=255)
                                                                                                 ccGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTg                                                >  2:335243/1‑149 (MQ=255)
                                                                                                  cgcgCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGc                                               >  1:358616/1‑149 (MQ=255)
                                                                                                           aaGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGCCCgg                                      >  5:11439/1‑149 (MQ=255)
                                                                                                            aGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCggg                                     >  8:65900/1‑149 (MQ=255)
                                                                                                              gTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGa                                   <  2:358616/149‑1 (MQ=255)
                                                                                                                           cggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCt                      >  8:43451/1‑149 (MQ=255)
                                                                                                                           cggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCt                      >  3:105589/1‑149 (MQ=255)
                                                                                                                             tccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg                    <  4:181905/148‑1 (MQ=255)
                                                                                                                             gccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg                    <  4:93290/149‑1 (MQ=255)
                                                                                                                              ccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTAAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCgg                   <  2:301975/149‑1 (MQ=255)
                                                                                                                                  cTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCa                           <  6:249742/137‑1 (MQ=255)
                                                                                                                                  cTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCa                           >  5:249742/1‑137 (MQ=255)
                                                                                                                                      tgCCAGCAGACGCTGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGc           >  6:52058/2‑149 (MQ=255)
                                                                                                                                          agcagcagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGTCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCTGGGAGCAGGCCAGCCCTCGGCACGCGGCACAg       <  1:248810/149‑1 (MQ=255)
                                                                                                                                           gcagccgcTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCggg                                     >  2:43739/1‑118 (MQ=255)
                                                                                                                                           gcagccgcTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCggg                                     <  1:43739/118‑1 (MQ=255)
                                                                                                                                            cagcagcGGCTGCTGGCAGCTTGCCTTGTGCTGGCTACGGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGTCCAGCCCTCGGCACGCGGCACAGGc     <  3:193485/149‑1 (MQ=255)
                                                                                                                                            cagcagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCCGGCGCAACTCTGGCGCAGCCAGGTCATTGACGTGGCAGGGGAGAATGACAGCCCTCGTCACGAGGCAAAGGc     >  8:257498/1‑149 (MQ=255)
                                                                                                                                               cagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCt  <  7:43451/149‑1 (MQ=255)
                                                                                                                                               cagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCt  <  8:142166/149‑1 (MQ=255)
                                                                                                                                                   ggCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCa                                >  7:261416/1‑115 (MQ=255)
                                                                                                                                                   ggCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCa                                <  8:261416/115‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCT  >  NC_002947/1202710‑1203001

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: