Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 740,791 T→C 7.4% intergenic (‑165/‑52) pilC ← / → pilA type IV fimbrial assembly protein PilC/fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947740,7910TC 7.4% 186.8 / 5.0 67intergenic (‑165/‑52)pilC/pilAtype IV fimbrial assembly protein PilC/fimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

CCCGGCGACAAGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATGTACACCAAC  >  NC_002947/740646‑740938
                                                                                                                                                 |                                                                                                                                                    
cccGGCGACAAGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTgg                                                                                                                                                   <  6:164443/149‑1 (MQ=255)
   ggCGACAAGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGgc                                                                                                                                                <  8:68897/149‑1 (MQ=255)
       aCAAGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCaa                                                                                                                                            <  5:30861/149‑1 (MQ=255)
          aGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCCCAACCGGCAGCGGCGCagtag                                                                                                                                         <  4:42601/149‑5 (MQ=255)
             cTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATTCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGAcc                                                                                                                                      >  2:5722/1‑149 (MQ=255)
                taaaCAGCTTGCGGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTGAGCGAAGATGTCTCCAGTGCGGGCTGGGGTGACAATAGGTGTCTGTTAAGGCCTACTGCGGCGGTGCCGGTTGGGGCGCAACTGACCt                                                                                                                                     <  2:268367/146‑1 (MQ=255)
                gAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTc                                                                                                                                    >  5:213777/1‑148 (MQ=255)
                     aGCTTACCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAAt                                                                                                                                 <  1:286226/146‑1 (MQ=255)
                     aGCTTACCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAAt                                                                                                                                 >  2:286226/1‑146 (MQ=255)
                             gcgTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGaa                                                                                                                      <  1:166454/149‑1 (MQ=255)
                               gTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAAt                                                                                                                     >  2:327824/1‑148 (MQ=255)
                                     gCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAg                                                                                                              >  8:22901/1‑149 (MQ=255)
                                     gCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAg                                                                                                              >  8:55306/1‑149 (MQ=255)
                                           gCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTa                                                                                                       <  2:148506/149‑1 (MQ=255)
                                             aCGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACg                                                                                                     >  7:294865/1‑149 (MQ=255)
                                               gcgTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGc                                                                                                                         >  6:19091/1‑128 (MQ=255)
                                               gcgTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGc                                                                                                                         <  5:19091/128‑1 (MQ=255)
                                               gcgTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGtt                                                                                                   <  8:267651/149‑1 (MQ=255)
                                               gcgTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCCCAACCGGCAGCGGCGCagtag                                                                                                                                         >  6:79762/1‑108 (MQ=38)
                                               gcgTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCCCAACCGGCAGCGGCGCagtag                                                                                                                                         <  5:79762/112‑5 (MQ=38)
                                                 gTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCa                                                                                                 >  5:157839/1‑149 (MQ=255)
                                                  tCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCAt                                                                                                <  7:22901/149‑1 (MQ=255)
                                                     gTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCACAATGTGTAAGCGAATCAAAAAGGAACTACGTTCAtca                                                                                              >  7:94285/1‑146 (MQ=255)
                                                        aTGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAgg                                                                                          >  2:232019/1‑149 (MQ=255)
                                                          gCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAggg                                                                                         >  1:140120/1‑148 (MQ=255)
                                                          gCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAggg                                                                                         >  3:132772/1‑148 (MQ=255)
                                                          gCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCCCAACCGGCAGCGGCGCagtag                                                                                                                                         <  6:237747/101‑5 (MQ=37)
                                                          gCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCCCAACCGGCAGCGGCGCagtag                                                                                                                                         >  5:237747/1‑97 (MQ=37)
                                                                tAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGtt                                                                                                                                                     <  5:101365/83‑1 (MQ=255)
                                                                tAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGtt                                                                                                                                                     >  6:101365/1‑83 (MQ=255)
                                                                  gagTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCgg                                                                                 <  7:55306/148‑1 (MQ=255)
                                                                     tagtTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATc                                                                             <  1:5722/149‑1 (MQ=255)
                                                                       gtTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAAACTACGTTCATGAAAGGGCAACGCgg                                                                                 <  1:14662/144‑1 (MQ=255)
                                                                       gtTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCgg                                                                                 >  2:14662/1‑143 (MQ=255)
                                                                         tCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCAcc                                                                          <  8:294865/148‑1 (MQ=255)
                                                                           aGCGAAGATGACGCCCGTGCGGGCATGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCt                                                                        >  6:321729/1‑148 (MQ=255)
                                                                           aGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCt                                                                        >  6:311967/1‑148 (MQ=255)
                                                                              gAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGAt                                                                     <  5:311967/148‑1 (MQ=255)
                                                                              gAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATc                                                                    >  3:95308/1‑149 (MQ=255)
                                                                                  cccctgCCCGCACGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTa                                                                               <  7:297838/129‑1 (MQ=255)
                                                                                  cccctgCCCGCACGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTa                                                                               >  8:297838/6‑134 (MQ=255)
                                                                                     aCGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACtga                                                             >  5:180858/1‑149 (MQ=255)
                                                                                     aCGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACtga                                                             <  2:140120/149‑1 (MQ=255)
                                                                                     aCGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACtga                                                             <  1:232019/149‑1 (MQ=255)
                                                                                          cGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGAtc                                                        <  6:213777/149‑1 (MQ=255)
                                                                                                 gCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGAtcgtcgtcg                                                 <  1:327824/149‑1 (MQ=255)
                                                                                                       gTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGatca                                           >  6:215661/1‑149 (MQ=255)
                                                                                                           cAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGACCTGATGAACGTCGTCGCGatcatccg                                       >  8:124548/1‑147 (MQ=255)
                                                                                                            aaaaGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCgg                                       >  2:345138/1‑148 (MQ=255)
                                                                                                                 gtgtCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCt                                 <  8:221315/149‑1 (MQ=255)
                                                                                                                        ttCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGa                            <  6:157839/147‑1 (MQ=255)
                                                                                                                           aGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATc                       <  4:132772/149‑1 (MQ=255)
                                                                                                                           aGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATc                       >  3:253031/1‑149 (MQ=255)
                                                                                                                            ggCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCg                      <  4:253031/149‑1 (MQ=255)
                                                                                                                            ggCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCg                      >  7:167303/1‑149 (MQ=255)
                                                                                                                               cTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGAGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCt                   >  4:298063/1‑149 (MQ=255)
                                                                                                                                 aCTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTa                 >  2:70259/1‑149 (MQ=255)
                                                                                                                                  cTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTa                 <  6:180858/148‑1 (MQ=255)
                                                                                                                                  cTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTAc                <  5:215661/149‑1 (MQ=255)
                                                                                                                                   tGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTAcc               >  6:115032/1‑149 (MQ=255)
                                                                                                                                     cgccgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGACCTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGa             >  7:78179/1‑149 (MQ=255)
                                                                                                                                     cgccgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGTCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCg              <  5:321729/148‑1 (MQ=255)
                                                                                                                                      gccgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGAt            <  1:70259/149‑1 (MQ=255)
                                                                                                                                       ccgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATg           <  8:167303/149‑1 (MQ=255)
                                                                                                                                        cgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGc                              <  7:247521/129‑1 (MQ=255)
                                                                                                                                        cgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGc                              >  8:247521/1‑129 (MQ=255)
                                                                                                                                        cgcTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATg           >  1:132489/1‑148 (MQ=255)
                                                                                                                                             ccGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATGTacac      <  4:95308/148‑1 (MQ=255)
                                                                                                                                                gTTGGGGCGCAACTGACCTCAATGAGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATGTACaccaac  <  3:298063/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCGGCGACAAGCCTGGAAACAGCTTGCCGCGTCCATGCAAGGGCACGCGTCCGTGATGCTCACTAGAGTAGTTCAGCGAAGATGACGCCCGTGCGGGCAGGGGTGACAAAAGGTGTCTGTTCAGGCCTACTGCGCCGCTGCCGGTTGGGGCGCAACTGACCTCAATGTGTAAGCGAATCAAAAAGGAACTACGTTCATGAAAGGGCAACGCGGTATCACCCTGATCGAACTGATGATCGTCGTCGCGATCATCGGCATTCTGGCGACCATCGCCTTACCGATGTACACCAAC  >  NC_002947/740646‑740938

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: