Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,735 T→A 5.2% L166Q (CTG→CAG)  pilF → type IV pili biogenesis protein PilF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947985,7350TA 5.2% 232.6 / 2.2 77L166Q (CTG→CAG) pilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base T (31/42);  new base A (2/2);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

CTTTATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGC  >  NC_002947/985587‑985881
                                                                                                                                                    |                                                                                                                                                  
cTTTATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCt                                                                                                                                                    >  3:289793/1‑149 (MQ=255)
 tttATGCTCAGGGGTGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTg                                                                                                                                                   >  5:252361/1‑149 (MQ=255)
  ttATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGt                                                                                                                                                  >  2:314198/1‑149 (MQ=255)
  ttATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGt                                                                                                                                                  >  3:55053/1‑149 (MQ=255)
  ttATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGACTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGTCCTGACCGCCCTGAAGCTCGAACGCAGCGATCAGGCGCACGCGTATCTGt                                                                                                                                                  >  7:70503/1‑149 (MQ=255)
     tGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGa                                                                                                                                               <  6:379501/149‑1 (MQ=255)
       cTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGaaa                                                                                                                                             <  1:339514/149‑1 (MQ=255)
        tCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAg                                                                                                                                            >  5:199204/1‑149 (MQ=255)
        tCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAg                                                                                                                                            >  5:238439/1‑149 (MQ=255)
        tCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAg                                                                                                                                            >  5:35892/1‑149 (MQ=255)
          aGGGGCGATTTGCCGAGGCCGAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCg                                                                                                                                          <  6:184285/149‑1 (MQ=255)
          aGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCg                                                                                                                                          <  3:179291/149‑1 (MQ=255)
          aGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCg                                                                                                                                          <  6:19305/149‑1 (MQ=255)
            gggCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCAATCAGGCGCACGCGTATCTGTTGAAAGCGtt                                                                                                                                        <  6:15302/149‑1 (MQ=255)
             ggCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTa                                                                                                                                       <  1:366513/149‑1 (MQ=255)
              gCGATTTGCCGAGGCCTAGCAGATGGTTCGCCTGGCCAGGGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTAc                                                                                                                                      <  8:213403/149‑1 (MQ=255)
                 aTTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAAc                                                                                                                                   <  6:342874/149‑1 (MQ=255)
                    tGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTc                                                                                                                                 <  5:334126/148‑1 (MQ=255)
                      ccGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCaa                                                                                                                               <  1:363287/148‑1 (MQ=255)
                      ccGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAAc                                                                                                                              <  6:238439/149‑1 (MQ=255)
                           gagaaGCAGATGTTTCGCCTTGGCTGTTCCGATACGCTTTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCg                                                                                                                         <  8:70503/146‑1 (MQ=255)
                           gCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCg                                                                                                                         <  6:2969/149‑1 (MQ=255)
                             cAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCgg                                                                                                                        >  5:136258/1‑148 (MQ=255)
                              aaGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGc                                                                                                                       >  3:201873/1‑148 (MQ=255)
                              aaGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCa                                                                                                                      <  4:289793/149‑1 (MQ=255)
                                 cAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAAc                                                                                                                    >  3:183427/1‑148 (MQ=255)
                                  aGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCATGAGCGCTCACGCGTGTACGAGAATATGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAAcc                                                                                                                   >  2:31722/1‑148 (MQ=255)
                                   gATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCa                                                                                                                  >  5:346660/1‑148 (MQ=255)
                                      gTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAg                                                                                                                           >  2:304054/1‑136 (MQ=255)
                                      gTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAg                                                                                                                           <  1:304054/136‑1 (MQ=255)
                                       tttCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGc                                                                                                             >  1:199697/1‑149 (MQ=255)
                                           gCCTGGCCAGTGACGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCTCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTg                                                                                                         >  8:198923/1‑149 (MQ=255)
                                            ccTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGc                                                                                                        <  1:37607/149‑1 (MQ=255)
                                                  cAGTGCCGATACCGTGTATCGTGTGCGCTCACGCGTGTAGGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCTACCAAAAGCGTTGCTGGaaa                                                                                                  <  4:55053/149‑1 (MQ=255)
                                                     tGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATgg                                                                                               <  1:114829/149‑1 (MQ=255)
                                                     tGCAGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATgg                                                                                               <  4:201873/149‑1 (MQ=255)
                                                       ccGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGc                                                                                              <  6:35892/148‑1 (MQ=255)
                                                        cGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGaaa                                                                                                  <  2:335838/143‑1 (MQ=255)
                                                        cGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGaaa                                                                                                  >  1:335838/1‑143 (MQ=255)
                                                               cTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTc                                                                                     <  3:361257/149‑1 (MQ=255)
                                                               cTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTc                                                                                     <  6:346660/149‑1 (MQ=255)
                                                               cTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCGGAAGCTCGAACGCCGCGATCAGGCGCACGCGGAGCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAACGCGTTGCGGGAAATGGCTgctgtgta                                                                                     >  5:389514/1‑142 (MQ=255)
                                                                  tATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTa                                                                                  <  1:33039/149‑1 (MQ=255)
                                                                  tATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCCTTACAACTCACCCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTa                                                                                  <  6:199204/149‑1 (MQ=255)
                                                                    tCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTAc                                                                                 >  4:67665/1‑148 (MQ=255)
                                                                       tGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAACGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGaa                                                                              >  6:354431/1‑148 (MQ=255)
                                                                                aCGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCAt                                                                     >  4:29925/1‑148 (MQ=255)
                                                                                aCGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCAt                                                                     >  7:216625/1‑148 (MQ=255)
                                                                                  gcgTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTAt                                                                  <  4:351497/149‑1 (MQ=255)
                                                                                    gtgtACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTAtgt                                                                <  1:314198/149‑1 (MQ=255)
                                                                                         cGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGATGCTGCAAATGGCTCAGTTGTCCTACGAAAACAGGCATAATGTGCCgg                                                           >  4:46611/1‑149 (MQ=255)
                                                                                         cGAGAATCTGGGCCTGAACGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACCCGTATATGTTGAAAGCGTTACAACTCAACCAGCGGCACCCAAAAGCGTTGCTGGCAATGGCTGAGTTGTCCTACGAACACAGGCATTATGTGCCgg                                                           >  6:299944/1‑149 (MQ=255)
                                                                                           agaATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGcc                                                         >  2:66032/1‑149 (MQ=255)
                                                                                            gaATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGccc                                                        >  7:350516/1‑149 (MQ=255)
                                                                                            gaATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGccc                                                        <  4:183427/149‑1 (MQ=255)
                                                                                              aTCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGTCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCgg                                                      <  6:332413/149‑1 (MQ=255)
                                                                                                 ttggCCTGTCCGCCGTGAAGCTCGAACGCGGCGATCAGGCGCACGCGTATCTGTTGAACGCGTTACAACTCAACCATCGGCAACCATAAGCGTTGCTGGAAATGGCTTAGTTGTCCTACGAAAACAGGCTTTATGTGCCGGCCCGGGa                                                    <  8:325791/146‑1 (MQ=255)
                                                                                                  gggCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCAGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTcc                                                                                    >  5:102317/1‑115 (MQ=255)
                                                                                                  gggCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCAGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTcc                                                                                    >  1:40544/1‑115 (MQ=255)
                                                                                                  gggCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCAGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTcc                                                                                    <  6:102317/115‑1 (MQ=255)
                                                                                                  gggCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCAGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTcc                                                                                    <  2:40544/115‑1 (MQ=255)
                                                                                                   ggCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGacta                                                 >  7:189455/1‑149 (MQ=255)
                                                                                                     ccTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGactact                                               <  2:199697/149‑1 (MQ=255)
                                                                                                      cTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGActacta                                              <  2:368255/149‑1 (MQ=255)
                                                                                                                  aaGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTc                                    <  8:350516/147‑1 (MQ=255)
                                                                                                                   aGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGc                                 <  1:87013/149‑1 (MQ=255)
                                                                                                                   aGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTGGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGGTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCacc                                 <  1:87010/149‑3 (MQ=255)
                                                                                                                     cTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCa                               <  4:82929/149‑1 (MQ=255)
                                                                                                                       cGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGt                             >  1:78619/1‑149 (MQ=255)
                                                                                                                       cGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGt                             >  1:213560/1‑149 (MQ=255)
                                                                                                                            gccgcGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCCACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGc                        >  7:297183/1‑149 (MQ=255)
                                                                                                                            gccgcGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGAGCGTTTCAGCCAGTTGAGc                        >  6:155452/1‑149 (MQ=255)
                                                                                                                             ccgcGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCg                       <  6:252361/149‑1 (MQ=255)
                                                                                                                               gcgATCAGGCGCACGCGTATCTGTTGAAAGGGTTACAACTCAACCAGCGGCAACCTAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGAc                     <  3:46611/149‑1 (MQ=255)
                                                                                                                                         gcACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTGGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGccc           <  6:136258/149‑1 (MQ=255)
                                                                                                                                            cgcgTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGAGGACGACGCCCGTa        <  1:162528/149‑1 (MQ=255)
                                                                                                                                             gcgTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAg       <  3:29925/149‑1 (MQ=255)
                                                                                                                                                  tCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCctgc  >  5:127767/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTTATGCTCAGGGGCGATTTGCCGAGGCCAAGCAGATGTTTCGCCTGGCCAGTGCCGATACCCTGTATCCTGAGCGCTCACGCGTGTACGAGAATCTGGGCCTGACCGCCCTGAAGCTCGAACGCCGCGATCAGGCGCACGCGTATCTGTTGAAAGCGTTACAACTCAACCAGCGGCAACCAAAAGCGTTGCTGGAAATGGCTGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGC  >  NC_002947/985587‑985881

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: ‑