Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,689,863 G→A 5.3% P11S (CCG→TCG)  pilH ← twitching motility protein PilH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,689,8630GA 5.3% 290.5 / 2.7 93P11S (CCG→TCG) pilHtwitching motility protein PilH
Reads supporting (aligned to +/- strand):  ref base G (48/40);  new base A (2/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCT  >  NC_002947/5689716‑5690011
                                                                                                                                                    |                                                                                                                                                     
cGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGc                                                                                                                                                      <  7:97686/149‑1 (MQ=255)
cGTTCATGCCAGGAATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGc                                                                                                                                                      <  2:155240/149‑1 (MQ=255)
 gTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGc                                                                                                                                                      >  4:29762/1‑148 (MQ=255)
  ttCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGc                                                                                                                                                      >  3:6594/1‑147 (MQ=255)
     aTGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtc                                                                                                                                                 <  4:262938/149‑1 (MQ=255)
              aTGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACaa                                                                                                                                         >  1:301329/1‑148 (MQ=255)
               tGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATc                                                                                                                                       <  3:225163/149‑1 (MQ=255)
               tGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATc                                                                                                                                       <  4:6594/149‑1 (MQ=255)
                gACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAAt                                                                                                                                        >  4:272499/1‑147 (MQ=255)
                 aCGATGTCCATCTGCACCGCGTCGGGCTTGTCCTGCCGGGCGAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATGTGTACATCTCTGTCGACGACTCGCCGagagagat                                                                                                                                     <  4:274297/149‑8 (MQ=255)
                      gTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTc                                                                                                                                >  5:395118/1‑149 (MQ=255)
                       tCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTc                                                                                                                                <  3:96236/148‑1 (MQ=255)
                        ccATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGACGACTCGCCGACagag                                                                                                                                       >  7:62978/1‑137 (MQ=255)
                        ccATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGACGACTCGCCGACagag                                                                                                                                       <  8:62978/140‑4 (MQ=255)
                         cATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTTGGCGAGTCGTCGACAATCAGAACTCggg                                                                                                                             >  6:245998/1‑149 (MQ=255)
                         cATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCggg                                                                                                                             >  5:18623/1‑149 (MQ=255)
                          aTCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGc                                                                                                                            >  3:75199/1‑149 (MQ=255)
                               cACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGc                                                                                                                       >  4:298043/1‑149 (MQ=255)
                                aCCGCGTCGGGCTTGTCCTGCCGGGCAAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCt                                                                                                                      >  6:342364/1‑149 (MQ=255)
                                 ccGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTg                                                                                                                     >  6:133205/1‑149 (MQ=255)
                                          gCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCAt                                                                                                            >  6:100595/1‑149 (MQ=255)
                                            ttGTCCTGCCGGGCCAGGGCCACGCCGTGGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCa                                                                                                         <  7:232638/149‑1 (MQ=255)
                                            ttGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCa                                                                                                         <  1:236739/149‑1 (MQ=255)
                                              gTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAg                                                                                                        >  1:73748/1‑148 (MQ=255)
                                               tCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAgg                                                                                                       >  2:183926/1‑148 (MQ=255)
                                                ccTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGa                                                                                                      >  7:245224/1‑148 (MQ=255)
                                                 cTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGaaa                                                                                                    <  4:342580/149‑1 (MQ=255)
                                                 cTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGa                                                                                                      >  4:70751/1‑147 (MQ=255)
                                                 cTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACACGGGCCAGGCTGTTTCCCCATCAGGa                                                                                                      >  6:348197/1‑147 (MQ=255)
                                                   gCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAg                                                                                                   >  5:11386/1‑148 (MQ=255)
                                                    ccGGGCGAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGa                                                                                                      >  5:233935/1‑144 (MQ=255)
                                                    ccGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGa                                                                                                      <  6:233935/144‑1 (MQ=255)
                                                    ccGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTa                                                                                                 <  6:419390/149‑1 (MQ=255)
                                                      gggccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAgg                                                                                               <  5:100595/149‑1 (MQ=255)
                                                       ggccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGa                                                                                              >  6:289113/1‑149 (MQ=255)
                                                       ggccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGa                                                                                              <  5:133205/149‑1 (MQ=255)
                                                        gccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGCACTCGGGCCCTGCTGTTTGCCCAGCAGGAAAGTAGGAc                                                                                             >  6:132683/1‑149 (MQ=255)
                                                         ccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCCCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGAcc                                                                                            >  8:174887/1‑149 (MQ=255)
                                                              gccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGcc                                                                                       <  5:149040/149‑1 (MQ=255)
                                                                 aCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATg                                                                                    >  1:216361/1‑149 (MQ=255)
                                                                 aCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGGCGACAATCAGAACTCGGGCCATGCTGTTTACCCATCAGGAAAGTAGGACCGCGCCATg                                                                                    >  7:307660/1‑149 (MQ=255)
                                                                   gCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGTCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATgg                                                                                   >  1:282038/1‑149 (MQ=255)
                                                                         ggCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCt                                                                            >  2:227846/1‑149 (MQ=255)
                                                                           cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTc                                                                           <  7:111807/148‑1 (MQ=255)
                                                                           cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTc                                                                           >  8:111807/1‑148 (MQ=255)
                                                                           cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCg                                                                          <  4:2482/149‑1 (MQ=255)
                                                                                   cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                  <  7:164772/149‑1 (MQ=255)
                                                                                   cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                  >  2:309577/1‑149 (MQ=255)
                                                                                   cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                  <  5:289113/149‑1 (MQ=255)
                                                                                    tGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTg                                                                 >  7:374329/1‑149 (MQ=255)
                                                                                      gCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCg                                                               <  3:70751/149‑1 (MQ=255)
                                                                                       ccTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCgt                                                              >  6:50206/1‑149 (MQ=255)
                                                                                         ttGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCgtgt                                                            >  2:183139/1‑149 (MQ=255)
                                                                                         ttGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCgtgt                                                            >  1:298416/1‑149 (MQ=255)
                                                                                           gAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGtgtt                                                           <  1:227846/148‑1 (MQ=255)
                                                                                            aGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGtgttg                                                          <  8:370420/148‑1 (MQ=255)
                                                                                             gCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTgttgtt                                                        <  6:18623/149‑1 (MQ=255)
                                                                                               aCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGa                                                      <  5:342364/149‑1 (MQ=255)
                                                                                               aCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGa                                                      <  2:301329/149‑1 (MQ=255)
                                                                                               aCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGa                                                      <  4:75199/149‑1 (MQ=255)
                                                                                                ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCgtgacgccgc                                                      <  1:309577/149‑8 (MQ=255)
                                                                                                ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                     <  3:272499/149‑1 (MQ=255)
                                                                                                ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                     <  1:183926/149‑1 (MQ=255)
                                                                                                ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                     >  3:60534/1‑149 (MQ=255)
                                                                                                    gTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTTGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGt                                                   <  5:245998/147‑1 (MQ=255)
                                                                                                    gTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGCGTCGTCGACAATCAGTACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGc                                                 <  5:348197/149‑1 (MQ=255)
                                                                                                           tGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAAc                                          <  1:183139/149‑1 (MQ=255)
                                                                                                                tttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCgg                                      >  8:168198/1‑149 (MQ=255)
                                                                                                                tttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCgg                                      >  2:170992/1‑149 (MQ=255)
                                                                                                                 ttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCAt                                                                                                            >  3:119711/1‑79 (MQ=255)
                                                                                                                 ttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCAt                                                                                                            <  4:119711/79‑1 (MQ=255)
                                                                                                                 ttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCg                                       <  6:11386/147‑1 (MQ=255)
                                                                                                                    cAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGca                                   >  6:167692/1‑148 (MQ=255)
                                                                                                                     aGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGcaca                                 >  8:94724/1‑149 (MQ=255)
                                                                                                                     aGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGcaca                                 <  6:395118/149‑1 (MQ=255)
                                                                                                                      gCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACAt                                >  7:324419/1‑149 (MQ=255)
                                                                                                                          ttCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                              >  3:297941/1‑147 (MQ=255)
                                                                                                                          ttCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                              <  4:297941/147‑1 (MQ=255)
                                                                                                                          ttCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCGCGGGCCGGGTGCGTGGTGTTGAGGTGCGGCAAACCCGGGcacagtggc                            >  8:261924/1‑144 (MQ=255)
                                                                                                                            cGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGgcgc                          >  2:166785/1‑149 (MQ=255)
                                                                                                                              gTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGTGGCAAACCCGGGCACATGGGcgcg                         >  8:298935/1‑148 (MQ=255)
                                                                                                                              gTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCgg                        <  2:282038/149‑1 (MQ=255)
                                                                                                                              gTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCAATGCTGGTTACACATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGCTGCGTGGTGTTGAGGTGCGGCAACCCGGGGCACATGGGcgcg                         >  8:177367/1‑148 (MQ=255)
                                                                                                                                 aaTCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATc                     >  4:4765/1‑149 (MQ=255)
                                                                                                                                 aaTCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATc                     <  8:324419/149‑1 (MQ=255)
                                                                                                                                  aTCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATcg                    <  1:170992/149‑1 (MQ=255)
                                                                                                                                          ctgattGTCGACGACTCGCCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                              <  4:219377/126‑1 (MQ=255)
                                                                                                                                          ctgattGTCGACGACTCGCCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                              >  3:219377/6‑131 (MQ=255)
                                                                                                                                          aTCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAg             >  3:64576/1‑148 (MQ=255)
                                                                                                                                           tctcTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGTAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAg             >  1:315639/1‑147 (MQ=255)
                                                                                                                                            ctctGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAg          >  1:158456/1‑149 (MQ=255)
                                                                                                                                               tGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAg          <  5:150024/146‑1 (MQ=255)
                                                                                                                                                 tCGGCGAGGCGGCGACAATCAGAACTCGGGCCCTGCTGGTTCCCCATCAGGAAAGTAGGCCCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGttctt     <  7:261924/149‑1 (MQ=255)
                                                                                                                                                   ggCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGttcttc    <  3:29762/148‑1 (MQ=255)
                                                                                                                                                    gCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCCGTTCTTCCt  >  1:138817/1‑149 (MQ=255)
                                                                                                                                                    gCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCt  <  7:168198/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCT  >  NC_002947/5689716‑5690011

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: