Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,689,867 G→C 5.4% D9E (GAC→GAG pilH ← twitching motility protein PilH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,689,8670GC 5.4% 300.9 / 6.6 92D9E (GAC→GAGpilHtwitching motility protein PilH
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base C (2/3);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

ATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTG  >  NC_002947/5689721‑5690013
                                                                                                                                                   |                                                                                                                                                   
aTGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAgtc                                                                                                                                                   <  4:262938/149‑1 (MQ=255)
         aTGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACaa                                                                                                                                           >  1:301329/1‑148 (MQ=255)
          tGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATc                                                                                                                                         <  4:6594/149‑1 (MQ=255)
          tGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATc                                                                                                                                         <  3:225163/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                                                                                                >  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                                                                             >  8:111807/1‑148 (MQ=255)
                                                                      cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTc                                                                             <  7:111807/148‑1 (MQ=255)
                                                                      cACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCg                                                                            <  4:2482/149‑1 (MQ=255)
                                                                              cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                    >  2:309577/1‑149 (MQ=255)
                                                                              cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                    <  5:289113/149‑1 (MQ=255)
                                                                              cTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGt                                                                    <  7:164772/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                                                        <  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:60534/1‑149 (MQ=255)
                                                                                           ccTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAg                                                       <  3:272499/149‑1 (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                                <  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                       >  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)
                                                                                                                                                 gAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGGACATGGGCGCGGATGGGGTCGAGCAGTTCTTCCTTg  >  8:347418/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGCCAGGCATGACGATGTCCATCAGCACCGCGTCGGGCTTGTCCTGCCGGGCCAGGGCCACGCCGTCGGCACCGTTGCTGGCCTTGAGCACCTGGTAGCCGTGCTTTTCCAGCCATTCGGTCAATCTGTACATCTCTGTCGGCGAGTCGTCGACAATCAGAACTCGGGCCATGCTGTTTCCCCATCAGGAAAGTAGGACCGCGCCATGGCGGGGCTCGGTCAGGGTGCGTGTTGTTGAGGTGCGGCAAACCCGGGCACATGGGCGCGGATCGCGTCGAGCAGTTCTTCCTTG  >  NC_002947/5689721‑5690013

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: