Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,689,726 A→T 5.6% P56P (CCT→CCA pilH ← twitching motility protein PilH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,689,7260AT 5.6% 269.6 / 4.9 89P56P (CCT→CCApilHtwitching motility protein PilH
Reads supporting (aligned to +/- strand):  ref base A (47/37);  new base T (3/2);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CCGGTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGA  >  NC_002947/5689578‑5689874
                                                                                                                                                    |                                                                                                                                                    
ccGGTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCa                                                                                                                                                      <  8:283161/149‑1 (MQ=255)
 cGGTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAg                                                                                                                                                     >  7:262469/1‑149 (MQ=255)
   gTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTAGCACGATCACTGGAATAGCGCTGGTGTGCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGc                                                                                                                                                   <  3:238535/149‑1 (MQ=255)
   gTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAgg                                                                                                                                                    >  2:338135/1‑148 (MQ=255)
        gTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGa                                                                                                                                               <  3:353592/148‑1 (MQ=255)
         tCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAg                                                                                                                                                     <  1:127446/141‑1 (MQ=255)
         tCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAg                                                                                                                                                     >  2:127446/1‑141 (MQ=255)
         tCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGAc                                                                                                                                              >  2:220142/1‑148 (MQ=255)
           aGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCTGGCATGaac                                                                                                                                             >  5:166669/1‑145 (MQ=255)
           aGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCTGGCATGaac                                                                                                                                             <  6:166669/147‑3 (MQ=255)
           aGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGAt                                                                                                                                           >  3:126714/1‑149 (MQ=255)
             gAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCTGGCATGaacg                                                                                                                                            <  6:289675/146‑4 (MQ=255)
             gAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCTGGCATGaacg                                                                                                                                            >  5:289675/1‑143 (MQ=255)
                          cccTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACcgc                                                                                                                            <  1:338135/149‑1 (MQ=255)
                                 gcgTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAAGCAAGACGATGTCCATCAGCACCGCGTCggga                                                                                                                     >  4:38348/1‑148 (MQ=255)
                                     ggCCCAGATACCGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGt                                                                                                                 <  8:127577/149‑1 (MQ=255)
                                      gCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTc                                                                                                                >  3:1306/1‑149 (MQ=255)
                                        ccAGATACGGTCGGTTTGCTGGTCCTTGGTGGTTACCACGATCACTGTAATAGAGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCt                                                                                                              <  4:55310/149‑1 (MQ=255)
                                        ccAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCt                                                                                                              >  2:246673/1‑149 (MQ=255)
                                                  tcggtTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccag                                                                                                    <  5:43939/149‑1 (MQ=255)
                                                     gtTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccagggc                                                                                                 >  6:9035/1‑149 (MQ=255)
                                                     gtTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccagggc                                                                                                 >  4:82290/1‑149 (MQ=255)
                                                      tttCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGATGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccagggc                                                                                                 >  3:81638/1‑148 (MQ=255)
                                                      tttCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccagcaac                                                                                                >  2:155850/1‑145 (MQ=255)
                                                      tttCCTGGTCCTTGGTGGGTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCgggccagggcc                                                                                                <  6:98291/149‑1 (MQ=255)
                                                        tCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCGATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCAc                                                                                              >  5:323700/1‑149 (MQ=255)
                                                        tCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCAc                                                                                              >  4:287382/1‑149 (MQ=255)
                                                        tCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCAc                                                                                              >  7:212431/1‑149 (MQ=255)
                                                          cTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCAGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGACGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGc                                                                                            >  5:34049/1‑149 (MQ=255)
                                                           tGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGc                                                                                                                                                   <  2:28766/93‑1 (MQ=255)
                                                           tGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGc                                                                                                                                                   >  1:28766/1‑93 (MQ=255)
                                                                 ttGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCgg                                                                                      >  1:205675/1‑148 (MQ=255)
                                                                  tggtggTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCAt                                                                                                                                     <  1:141244/100‑1 (MQ=255)
                                                                  tggtggTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCAt                                                                                                                                     >  2:141244/1‑100 (MQ=255)
                                                                      ggtTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGAAGTCCATCAGCGCCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGt                                                                                >  1:158462/1‑149 (MQ=255)
                                                                       gtTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGATGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGt                                                                                <  4:209157/148‑1 (MQ=255)
                                                                           ccACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCt                                                                            <  5:9035/148‑1 (MQ=255)
                                                                           ccACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTg                                                                           <  2:38108/149‑1 (MQ=255)
                                                                           ccACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTg                                                                           >  7:264148/1‑149 (MQ=255)
                                                                           ccACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTg                                                                           >  4:50019/1‑149 (MQ=255)
                                                                           ccACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTAGGCACCGTTGCTg                                                                           >  1:84163/1‑149 (MQ=255)
                                                                              cGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCg                                                                                       >  6:47703/1‑134 (MQ=255)
                                                                              cGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCg                                                                                       <  5:47703/134‑1 (MQ=255)
                                                                               gATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCt                                                                       <  4:126714/149‑1 (MQ=255)
                                                                                 tCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCtt                                                                      >  5:176515/1‑148 (MQ=255)
                                                                                   aCTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGa                                                                    <  1:220142/148‑1 (MQ=255)
                                                                                     tGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCa                                                                 >  2:218004/1‑149 (MQ=255)
                                                                                       gAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCAcc                                                               >  4:254984/1‑149 (MQ=255)
                                                                                        aaTAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCg                                                                                       >  7:159964/1‑124 (MQ=255)
                                                                                        aaTAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCg                                                                                       <  8:159964/124‑1 (MQ=255)
                                                                                        aaTAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCt                                                              >  3:150510/1‑149 (MQ=255)
                                                                                         aTAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTg                                                             >  5:177596/1‑149 (MQ=255)
                                                                                              gcTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGGCGGCACCGTTGATGGCCTTGCGCAACTGGTAGc                                                        >  8:307260/1‑149 (MQ=255)
                                                                                                   tgtCCGTGTCCTTGCTGAGCTGGCGGGTTGCCTGGAAGCAGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGTCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGc                                                   <  1:240782/149‑1 (MQ=255)
                                                                                                    gtCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGc                                                        >  1:351721/1‑143 (MQ=255)
                                                                                                    gtCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGc                                                        <  2:351721/143‑1 (MQ=255)
                                                                                                        gggTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTc                                              <  8:264148/149‑1 (MQ=255)
                                                                                                        gggTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTc                                              <  8:212431/149‑1 (MQ=255)
                                                                                                          gTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCa                                            <  4:150510/149‑1 (MQ=255)
                                                                                                          gTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGAACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGAACAGTTGCTGGCCTTGAGAACCTGGTATCCTTGCTTTTc                                              >  7:136132/1‑147 (MQ=255)
                                                                                                           tCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAg                                           <  3:287382/149‑1 (MQ=255)
                                                                                                             cTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGcc                                         <  3:82290/149‑1 (MQ=255)
                                                                                                             cTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCACCACCGCGACGGGCTTGTCCTGCCGGGCCAGCGCCACGCCGTCGCCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGAGCTTTTCCAGc                                          >  2:84424/1‑148 (MQ=255)
                                                                                                                gCTGAGCTGGCGGGTGGCCTGGATGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCAGTCGGCACCTTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCa                                        <  4:81638/147‑1 (MQ=255)
                                                                                                                 cTGAGCTGGCGGGTGGCCTGGAAGCCGTTGATGGCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCAtt                                      <  2:125267/148‑1 (MQ=255)
                                                                                                                 cTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTc                                     <  2:84163/149‑1 (MQ=255)
                                                                                                                  tGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTc                                     >  7:111787/1‑148 (MQ=255)
                                                                                                                    aGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTc                                     <  8:111787/146‑1 (MQ=255)
                                                                                                                    aGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCg                                    <  3:254984/147‑1 (MQ=255)
                                                                                                                       tGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCaa                               >  6:356204/1‑149 (MQ=255)
                                                                                                                         gCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTc                                     <  8:63805/141‑1 (MQ=255)
                                                                                                                         gCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTc                                     >  7:63805/1‑141 (MQ=255)
                                                                                                                         gCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATc                             >  8:316203/1‑149 (MQ=255)
                                                                                                                             gTGGCCTGGAAGCCGTTCATGCCATGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTTGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTa                         <  3:38348/149‑1 (MQ=255)
                                                                                                                             gTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTa                         <  6:34049/149‑1 (MQ=255)
                                                                                                                                 ccTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACAtc                     >  3:103490/1‑149 (MQ=255)
                                                                                                                                  cTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGCGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTATCCAGCCATTCGGTCAATCTGTACAtct                    <  4:1306/149‑1 (MQ=255)
                                                                                                                                     gAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTg                 >  2:216085/1‑149 (MQ=255)
                                                                                                                                      aaGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTa                                                          >  2:257996/1‑107 (MQ=255)
                                                                                                                                       aGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTc               >  2:31662/1‑149 (MQ=255)
                                                                                                                                        gCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCg              <  1:246673/149‑1 (MQ=255)
                                                                                                                                            ttgaTGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGa          <  8:262469/146‑1 (MQ=255)
                                                                                                                                            ttCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGt                          <  2:171062/133‑1 (MQ=255)
                                                                                                                                            ttCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGt                          >  1:171062/1‑133 (MQ=255)
                                                                                                                                            ttCATGCCAGGCAGGACGATGTCCTTCAGCCCCGCTGCGGGCTTGTCCTGCCGGGCGAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCGAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGa          <  2:158462/149‑1 (MQ=255)
                                                                                                                                            gtCATGCCTGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTACTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGc            >  3:319977/2‑147 (MQ=255)
                                                                                                                                               aTGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtc       >  6:249033/1‑149 (MQ=255)
                                                                                                                                                   cAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCATGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtcgtc    >  7:201887/1‑148 (MQ=255)
                                                                                                                                                   cAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtcgtc    >  8:94570/1‑148 (MQ=255)
                                                                                                                                                   cAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtcgtc    >  4:160647/1‑148 (MQ=255)
                                                                                                                                                    aGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGa  <  5:356204/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGTTTGGTCAGGAAGTCACGGGCACCCTGGCGCGTGGCCCAGATACGGTCGGTTTCCTGGTCCTTGGTGGTTACCACGATCACTGGAATAGCGCTGGTGTCCGGGTCCTTGCTGAGCTGGCGGGTGGCCTGGAAGCCGTTCATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGA  >  NC_002947/5689578‑5689874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: