Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,689,871 T→C 5.6% D8G (GAC→GGC)  pilH ← twitching motility protein PilH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,689,8710TC 5.6% 270.0 / 4.1 92D8G (GAC→GGC) pilHtwitching motility protein PilH
Reads supporting (aligned to +/- strand):  ref base T (47/38);  new base C (2/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTGCTGAAC  >  NC_002947/5689730‑5690019
                                                                                                                                              |                                                                                                                                                     
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                                                                                                                                               <  8:62978/140‑4 (MQ=255)
          ccATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGACGACTCGCCGACagag                                                                                                                                               >  7:62978/1‑137 (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                                                                                                      <  5:133205/149‑1 (MQ=255)
                                         ggccagggccaCGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGa                                                                                                      >  6:289113/1‑149 (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                                                                                   >  8:111807/1‑148 (MQ=255)
                                                             cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTc                                                                                   <  7:111807/148‑1 (MQ=255)
                                                             cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCg                                                                                  <  4:2482/149‑1 (MQ=255)
                                                                     cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                          <  7:164772/149‑1 (MQ=255)
                                                                     cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                          <  5:289113/149‑1 (MQ=255)
                                                                     cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                          >  2:309577/1‑149 (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                                                              <  4:75199/149‑1 (MQ=255)
                                                                                 aCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGa                                                              <  2:301329/149‑1 (MQ=255)
                                                                                 aCCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGa                                                              <  5:342364/149‑1 (MQ=255)
                                                                                  ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCgtgacgccgc                                                              <  1:309577/149‑8 (MQ=255)
                                                                                  ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                             <  1:183926/149‑1 (MQ=255)
                                                                                  ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                             <  3:272499/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                                                                                                                    <  4:119711/79‑1 (MQ=255)
                                                                                                   ttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCAt                                                                                                                    >  3:119711/1‑79 (MQ=255)
                                                                                                   ttCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCg                                               <  6:11386/147‑1 (MQ=255)
                                                                                                      cAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGca                                           >  6:167692/1‑148 (MQ=255)
                                                                                                       aGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGcaca                                         <  6:395118/149‑1 (MQ=255)
                                                                                                       aGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGcaca                                         >  8:94724/1‑149 (MQ=255)
                                                                                                        gCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACAt                                        >  7:324419/1‑149 (MQ=255)
                                                                                                            ttCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                                      <  4:297941/147‑1 (MQ=255)
                                                                                                            ttCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATgg                                      >  3:297941/1‑147 (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                             <  8:324419/149‑1 (MQ=255)
                                                                                                                   aaTCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATc                             >  4:4765/1‑149 (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)
                                                                                                                                        gAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGGACATGGGCGCGGATGGGGTCGAGCAGTTCTTCCTTg        >  8:347418/1‑149 (MQ=255)
                                                                                                                                           tcgtcgACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCa                                                                              >  5:224151/1‑76 (MQ=255)
                                                                                                                                           tcgtcgACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCa                                                                              <  6:224151/76‑1 (MQ=255)
                                                                                                                                           tcgtcgACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGGGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTGCt      <  7:298935/148‑1 (MQ=255)
                                                                                                                                             gtcgACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTGCTGaa   <  8:374329/149‑1 (MQ=255)
                                                                                                                                              tcgACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTGCTGAAc  <  7:174887/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
ATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTGCTGAAC  >  NC_002947/5689730‑5690019

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: