Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,202,854 A→C 6.1% Q228P (CAG→CCG)  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,8540AC 6.1% 203.8 / 3.5 65Q228P (CAG→CCG) xcpYtype II secretion pathway protein XcpY
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base C (2/2);  total (35/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 = 9.51e-01

CGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCT  >  NC_002947/1202708‑1203001
                                                                                                                                                  |                                                                                                                                                     
cGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCagcagcagc                                                                                                                                                     >  8:91653/1‑149 (MQ=255)
  tGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCgg                                                                                                                                                   <  6:159885/149‑1 (MQ=255)
          ccACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTgg                                                                                                                                            >  4:312791/1‑148 (MQ=255)
           cACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGcg                                                                                                                                          <  8:122089/149‑2 (MQ=255)
                   gCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCGCCCGCCAGCAGCAGCGGCTGCTGGCAGCGTGCCt                                                                                                                                  >  3:201594/1‑149 (MQ=255)
                     aGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTgg                                                                                                                                <  1:193947/149‑1 (MQ=255)
                         ctgcACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCt                                                                                                                            >  7:2388/1‑149 (MQ=255)
                                 gACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCCGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCg                                                                                                                    >  4:232983/1‑149 (MQ=255)
                                     tCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTg                                                                                                                >  2:14933/1‑149 (MQ=255)
                                      cGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                               <  2:224124/149‑1 (MQ=255)
                                       gCTGAGCGGCAAGTGGGCAACCGTATTTGCGCCGCGCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTggg                                                                                                              <  4:201594/149‑1 (MQ=255)
                                            gCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGtt                                                                                                                  >  4:239479/1‑140 (MQ=255)
                                            gCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGtt                                                                                                                  <  3:239479/140‑1 (MQ=255)
                                              ggCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCt                                                                                                       >  1:219396/1‑149 (MQ=255)
                                               gCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtg                                                                                                      <  7:91653/149‑1 (MQ=255)
                                                     ggCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCGCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                                 <  3:289967/148‑1 (MQ=255)
                                                       ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGcg                                                                                                                                          <  6:67333/105‑2 (MQ=255)
                                                       ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGcg                                                                                                                                          >  5:67333/1‑104 (MQ=255)
                                                       ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAg                                                                                              >  7:36163/1‑149 (MQ=255)
                                                        cAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                             >  5:61452/1‑149 (MQ=255)
                                                        cAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                             <  8:6996/149‑1 (MQ=255)
                                                          aCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGc                                                                                                         <  1:312578/135‑1 (MQ=255)
                                                          aCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGc                                                                                                         >  2:312578/1‑135 (MQ=255)
                                                            cGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGc                                                                                         <  7:5940/149‑1 (MQ=255)
                                                               aTTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                                    <  4:101734/135‑1 (MQ=255)
                                                               aTTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                                    >  3:101734/1‑135 (MQ=255)
                                                                     gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGt                                                                                 >  6:295703/1‑148 (MQ=255)
                                                                      ctgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGAGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGt                                                                                 >  3:236230/3‑147 (MQ=255)
                                                                         cTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                                 <  6:237931/128‑1 (MQ=255)
                                                                         cTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                                 >  5:237931/1‑128 (MQ=255)
                                                                              cacacTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                               <  4:179250/125‑1 (MQ=255)
                                                                              cacacTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                               >  3:179250/1‑125 (MQ=255)
                                                                                 acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTggg                                                                                                              >  8:4708/1‑107 (MQ=255)
                                                                                 acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTggg                                                                                                              <  7:4708/107‑1 (MQ=255)
                                                                                 acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AActct                                                                  >  7:125093/1‑149 (MQ=255)
                                                                                     cccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGgc                                                               <  4:236230/148‑1 (MQ=255)
                                                                                     cccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGgc                                                               <  2:256515/148‑1 (MQ=255)
                                                                                         gCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGc                                                          >  7:42603/1‑149 (MQ=255)
                                                                                              ttCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGgc                                                               <  4:47713/139‑1 (MQ=255)
                                                                                              ttCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGgc                                                               >  3:47713/1‑139 (MQ=255)
                                                                                                 gACCGCGCCCGCAAGGTGCGCTCGTTGCCGACCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCAt                                                  <  6:220708/149‑1 (MQ=255)
                                                                                                    cgcgCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGc                                               >  5:160076/1‑149 (MQ=255)
                                                                                                    cgcgCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGc                                               >  8:152376/1‑149 (MQ=255)
                                                                                                     gcgcCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGcc                                              <  6:160076/149‑1 (MQ=255)
                                                                                                        cccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCa                                                                                        >  7:91868/1‑106 (MQ=255)
                                                                                                        cccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCa                                                                                        <  8:91868/106‑1 (MQ=255)
                                                                                                        cccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTg                                           <  2:219396/149‑1 (MQ=255)
                                                                                                            cAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCAAGCAGCAGCGGATCCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTCGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACACTGGAGCAGCCAGGTCATTGCCGTgcacg                                       >  8:259912/1‑145 (MQ=255)
                                                                                                              aGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCggg                                     <  5:295703/149‑1 (MQ=255)
                                                                                                                 tGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAg                                  <  8:2388/149‑1 (MQ=255)
                                                                                                                    gcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCCCCAAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAg                                  >  1:107806/1‑148 (MQ=255)
                                                                                                                               gccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGc                                 <  3:295748/136‑1 (MQ=255)
                                                                                                                               gccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGc                                 >  4:295748/1‑136 (MQ=255)
                                                                                                                               gccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg                    <  7:152376/149‑1 (MQ=255)
                                                                                                                                 cgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGt                                            >  4:122132/1‑123 (MQ=255)
                                                                                                                                 cgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGt                                            <  3:122132/123‑1 (MQ=255)
                                                                                                                                 cgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGc                  <  6:61452/149‑1 (MQ=255)
                                                                                                                                  gccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCa                 >  7:57689/1‑149 (MQ=255)
                                                                                                                                  gccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCa                 >  1:261014/1‑149 (MQ=255)
                                                                                                                                    cTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCAcg               >  8:154188/1‑149 (MQ=255)
                                                                                                                                    cTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCAcg               >  5:76089/1‑149 (MQ=255)
                                                                                                                                        tgCCAGCAGCCGCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCgg            >  8:48206/2‑148 (MQ=255)
                                                                                                                                        cGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGc           >  4:223830/1‑149 (MQ=255)
                                                                                                                                        cGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGc           <  7:90460/149‑1 (MQ=255)
                                                                                                                                         gCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGca          >  7:179749/1‑149 (MQ=255)
                                                                                                                                          ccAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGca          >  4:44715/1‑148 (MQ=255)
                                                                                                                                              cagcagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGc     <  8:36163/149‑1 (MQ=255)
                                                                                                                                               agcagcGCCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGcc    >  5:230704/1‑149 (MQ=255)
                                                                                                                                                 cagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCt  <  8:179749/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGC‑‑AACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCT  >  NC_002947/1202708‑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: