Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,260,749 T→C 5.9% V446A (GTC→GCC)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,260,7490TC 5.9% 197.7 / 2.3 67V446A (GTC→GCC) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCCGCGCTGCTG  >  NC_002947/2260606‑2260897
                                                                                                                                               |                                                                                                                                                    
ccTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAg                                                                                                                                                 >  3:268862/1‑149 (MQ=255)
          gCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGcc                                                                                                                                       >  7:138134/1‑149 (MQ=255)
          gCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGCTCTCGat                                                                                                                                         <  3:186331/147‑2 (MQ=255)
            aCTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCgg                                                                                                                                     <  2:117066/149‑1 (MQ=255)
              tgctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCCAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGgt                                                                                                                                   >  3:144056/1‑149 (MQ=255)
                ctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGgt                                                                                                                                   <  1:11618/147‑1 (MQ=255)
                ctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGgt                                                                                                                                   <  2:390477/147‑1 (MQ=255)
                   ctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAc                                                                                                                              <  8:266512/149‑1 (MQ=255)
                     gtgGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAccc                                                                                                                            >  4:66919/1‑149 (MQ=255)
                      tgGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCt                                                                                                                           <  7:41443/149‑1 (MQ=255)
                      tgGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCt                                                                                                                           >  3:288662/1‑149 (MQ=255)
                         ctgctgGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTc                                                                                                                        <  7:119080/149‑1 (MQ=255)
                                     aaaCGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGt                                                                                                            >  2:263558/1‑149 (MQ=255)
                                          cAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCg                                                                                                       <  1:162461/149‑1 (MQ=255)
                                           aaGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCtgt                                                                                                               >  1:196242/1‑140 (MQ=255)
                                           aaGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCtgt                                                                                                               <  2:196242/140‑1 (MQ=255)
                                              gCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTgg                                                                                                             >  6:124169/1‑139 (MQ=255)
                                              gCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTgg                                                                                                             <  5:124169/139‑1 (MQ=255)
                                                ccAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCgg                                                                                                 >  5:72352/1‑149 (MQ=255)
                                                ccAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCgg                                                                                                 >  6:312692/1‑149 (MQ=255)
                                                  aGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCt                                                                                               <  4:268862/149‑1 (MQ=255)
                                                  aGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCt                                                                                               >  8:384666/1‑149 (MQ=255)
                                                     aagaagCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCa                                                                                            <  8:138134/149‑1 (MQ=255)
                                                       gaagCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGc                                                                                          >  7:172333/1‑149 (MQ=255)
                                                        aagCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctct                                                                                                                                             >  4:8066/1‑94 (MQ=255)
                                                        aagCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctct                                                                                                                                             <  3:8066/97‑4 (MQ=255)
                                                         agCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTg                                                                                        >  3:56091/1‑149 (MQ=255)
                                                          gCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGc                                                                                       >  2:15943/1‑149 (MQ=255)
                                                             gagaAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGt                                                                                    >  4:154993/1‑149 (MQ=255)
                                                             gagaAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGt                                                                                    <  7:384666/149‑1 (MQ=255)
                                                                aaGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGc                                                                                 <  3:66919/149‑1 (MQ=255)
                                                                          gcATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGaaa                                                                         >  5:300596/1‑147 (MQ=255)
                                                                           cATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAg                                                                      <  4:56091/149‑1 (MQ=255)
                                                                            aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGt                                                                                                                                                      >  7:25708/1‑68 (MQ=255)
                                                                            aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGt                                                                                                                                                      <  8:25708/68‑1 (MQ=255)
                                                                            aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTc                                                                                   >  1:19706/1‑135 (MQ=255)
                                                                            aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTc                                                                                   <  2:19706/135‑1 (MQ=255)
                                                                            aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGt                                                                     >  8:320510/1‑149 (MQ=255)
                                                                             tGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTg                                                                    <  7:172646/149‑1 (MQ=255)
                                                                               gcgcGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                                  <  4:288662/149‑1 (MQ=255)
                                                                                  cgGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGa                                                               <  6:300596/149‑1 (MQ=255)
                                                                                   gggCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGt                                                                     >  1:29481/1‑142 (MQ=255)
                                                                                   gggCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGt                                                                     <  2:29481/142‑1 (MQ=255)
                                                                                        cTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGt                                                         <  8:172333/149‑1 (MQ=255)
                                                                                           gaagaaGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCgg                                                                            >  7:34871/1‑127 (MQ=255)
                                                                                           gaagaaGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCgg                                                                            <  8:34871/127‑1 (MQ=255)
                                                                                           gaagaaGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGc                                                      >  2:187069/1‑149 (MQ=255)
                                                                                              gaagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGcc                                                     >  1:66970/1‑147 (MQ=255)
                                                                                               aagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAg                                                  <  7:320510/149‑1 (MQ=255)
                                                                                                  agCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGACAGAACCGGTCGCCGAGgca                                               <  1:15943/149‑1 (MQ=255)
                                                                                                          gcAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGc                                                                                                         <  1:360384/83‑1 (MQ=255)
                                                                                                          gcAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGc                                                                                                         >  2:360384/1‑83 (MQ=255)
                                                                                                            aGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGATGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGc                                     >  5:208836/1‑149 (MQ=255)
                                                                                                              gTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCt                                   >  5:189414/1‑149 (MQ=255)
                                                                                                                  cGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGaa                                <  6:189414/148‑1 (MQ=255)
                                                                                                                   gACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                                  <  7:92384/113‑1 (MQ=255)
                                                                                                                   gACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                                  >  8:92384/1‑113 (MQ=255)
                                                                                                                     cAGCGATATCGAGAGCTTCGGAGTCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACcac                            <  4:144056/149‑1 (MQ=255)
                                                                                                                         gATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTg                                                                    <  5:233537/105‑1 (MQ=255)
                                                                                                                         gATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTg                                                                    >  6:233537/1‑105 (MQ=255)
                                                                                                                           tatCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTgg                                                                   <  2:388830/104‑1 (MQ=255)
                                                                                                                             tcgggggCTTCGACGCCGCCGAAGCTCTCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGAc                    >  8:10519/7‑149 (MQ=255)
                                                                                                                               gagaGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGAccc                  <  3:154993/149‑1 (MQ=255)
                                                                                                                                        ggcggcGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCAcgccgc         <  6:150554/149‑1 (MQ=255)
                                                                                                                                         gcggcgTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCcgcg        >  4:30650/1‑149 (MQ=255)
                                                                                                                                            gcgTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCCGCgctg     >  6:342469/1‑149 (MQ=255)
                                                                                                                                              gTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCCGCgctgc    <  2:398066/148‑1 (MQ=255)
                                                                                                                                              gTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCCGCgctgc    >  1:398066/1‑148 (MQ=255)
                                                                                                                                              gTCGAAGCCCCCGAGCCGGGTCTAACCCTGTCGCCGGCTCTCGTTCACGCTTCGGCTTCAGATGCCGTCGCTTCGGAAAAAGTGGCAGAACCGGTCGCCGAGCCAGCACCCGAGCCTGAACAACAAGCAGACCCGCACGCCGCgctgct   >  2:66836/1‑149 (MQ=255)
                                                                                                                                               tCGAAGCCCCCGAGCCGGGTGTAACCCTTTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCTCAAGCAGACCCGCACGCCGCGctgctg  <  1:66836/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCCGCACGCCGCGCTGCTG  >  NC_002947/2260606‑2260897

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