Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 925,343 G→A 5.6% A428V (GCC→GTC)  PP_0805 ← outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947925,3430GA 5.6% 214.8 / 2.2 71A428V (GCC→GTC) PP_0805outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base A (2/2);  total (33/38)
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.59e-01

GGATCTCCACACCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCTTGGCCTCCGCCAG  >  NC_002947/925197‑925491
                                                                                                                                                   |                                                                                                                                                     
ggATCTCCACACCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCg                                                                                                                                                     >  4:289668/1‑149 (MQ=255)
ggATCTCCACACCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCg                                                                                                                                                     >  7:314390/1‑149 (MQ=255)
 gATCTCCACACCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGt                                                                                                                                                    >  4:340037/1‑149 (MQ=255)
       cacaCCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAg                                                                                                                                                >  4:205730/1‑147 (MQ=255)
       cacaCCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGc                                                                                                                                               <  4:359565/148‑1 (MQ=255)
        acacCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAg                                                                                                                                             <  3:340037/149‑1 (MQ=255)
            cTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCAACGCCGAGCAggaga                                                                                                                                         <  7:165276/149‑4 (MQ=255)
             ttGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTc                                                                                                                                          >  7:332329/1‑147 (MQ=255)
             ttGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTc                                                                                                                                          >  4:340339/1‑147 (MQ=255)
             ttGAAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGAGTTGAGCAGGTc                                                                                                                                          >  1:275287/1‑147 (MQ=255)
               gCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAg                                                                                                                                       <  6:158902/148‑1 (MQ=255)
                cAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCa                                                                                                                                     <  6:172335/149‑1 (MQ=255)
                   gTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACGAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGt                                                                                                                                   >  1:278770/1‑149 (MQ=255)
                        gcCCCTCAAGGCCTAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             <  1:366099/149‑1 (MQ=255)
                        gcCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             <  2:278770/149‑1 (MQ=255)
                        gcCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             <  6:382853/149‑1 (MQ=255)
                          cccTCTAGGCCAAAGGCCGTGCGCAGGGAGCCGGTGTTGTACAGGCAGTCGACATCCTGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGTTAGATCTCCTGCTCTGCGTTGAGCATGTCGAGCAGTGGCCgt                                                                                                                             <  2:124559/147‑1 (MQ=255)
                          cccTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCCGTCGACATCCAGGCGCAGCACATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             <  2:311581/147‑1 (MQ=255)
                          cccTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgtgt                                                                                                                           >  6:50157/1‑149 (MQ=255)
                          cccTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTAAAGGCAGTCGACATAAAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgtg                                                                                                                            >  1:124559/1‑148 (MQ=255)
                           ccTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTc                                                                                                                          <  2:98707/149‑1 (MQ=255)
                               aaGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             <  3:290377/142‑1 (MQ=255)
                               aaGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCgt                                                                                                                             >  4:290377/1‑142 (MQ=255)
                                               cgcAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGAc                                                                                                      <  8:314390/148‑1 (MQ=255)
                                                   gggCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCg                                                                                                 <  5:360225/149‑1 (MQ=255)
                                                    ggCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGt                                                                                                <  2:60316/149‑1 (MQ=255)
                                                     gcgcCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCTTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTa                                                                                               <  5:232776/149‑1 (MQ=255)
                                                       gcCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTaga                                                                                             <  8:332329/149‑1 (MQ=255)
                                                        ccGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTagag                                                                                            >  4:11444/1‑149 (MQ=255)
                                                           gtgtTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGtt                                                                                                                                                   <  8:68225/92‑1 (MQ=255)
                                                           gtgtTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAg                                                                                                              <  2:14205/128‑1 (MQ=255)
                                                           gtgtTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTACTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAg                                                                                                              >  1:14205/1‑128 (MQ=255)
                                                            tgttgtACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTgg                                                                                                                                 <  6:39122/109‑1 (MQ=255)
                                                            tgttgtACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTgg                                                                                                                                 >  5:39122/1‑109 (MQ=255)
                                                              ttgtACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCg                                                                                      <  3:289668/149‑1 (MQ=255)
                                                                gtACATGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGt                                                                                     <  2:125343/148‑1 (MQ=255)
                                                                       gCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTaga                                                                                             <  5:86895/134‑1 (MQ=255)
                                                                       gCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTaga                                                                                             >  6:86895/1‑134 (MQ=255)
                                                                             gACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATg                                                                         <  2:212861/148‑1 (MQ=255)
                                                                               cATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGCAGAGCATGCGTGCCTCCTGGATGCTg                                                                      <  1:10285/149‑1 (MQ=255)
                                                                                aTCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTgg                                                                                                                                 <  8:103457/90‑1 (MQ=255)
                                                                                aTCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGc                                                                                           <  5:85392/127‑1 (MQ=255)
                                                                                aTCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGc                                                                                           >  6:85392/1‑127 (MQ=255)
                                                                                 tCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGAc                                                                    <  7:336928/149‑1 (MQ=255)
                                                                                      gcgcAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGc                                                               <  7:29738/149‑1 (MQ=255)
                                                                                       cgcAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCg                                                              <  4:193663/149‑1 (MQ=255)
                                                                                           gcagATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGctgact                                                          >  8:20049/1‑144 (MQ=255)
                                                                                           gcagATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCt                                                          >  5:144287/1‑149 (MQ=255)
                                                                                              gATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTc                                                        >  4:167560/1‑148 (MQ=255)
                                                                                              gATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTc                                                        >  5:339551/1‑148 (MQ=255)
                                                                                              gATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCATAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTc                                                        >  1:206790/1‑148 (MQ=255)
                                                                                                  ggAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGc                                                    >  8:152468/1‑148 (MQ=255)
                                                                                                   gAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCg                                                   >  3:91204/1‑148 (MQ=255)
                                                                                                   gAGCGGGTGCCGGCAAGATCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCg                                                   >  2:82248/1‑148 (MQ=255)
                                                                                                    aGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGaa                                                 <  3:340339/149‑1 (MQ=255)
                                                                                                     cggggTGCCGGCAAGTTGGAAGGGGGACTGGTACAACTCCTGGTGGGGGTTGACCAGGTCCTCCATGGGCCGGTTCCCTAGGTACCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGaa                                                 <  1:113679/147‑1 (MQ=255)
                                                                                                           gCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCg                                                              <  8:284058/129‑1 (MQ=255)
                                                                                                           gCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCg                                                              >  7:284058/1‑129 (MQ=255)
                                                                                                             cGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCcgg                                        >  4:24196/1‑149 (MQ=255)
                                                                                                                            gTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAg                         <  4:91204/149‑1 (MQ=255)
                                                                                                                               aCTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGtt                      <  5:50157/149‑1 (MQ=255)
                                                                                                                                 tGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTg                    >  2:319507/1‑149 (MQ=255)
                                                                                                                                   gTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAg                         >  1:193103/1‑142 (MQ=255)
                                                                                                                                   gTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAg                         <  2:193103/142‑1 (MQ=255)
                                                                                                                                    tAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTgcgc                 <  3:205730/149‑1 (MQ=255)
                                                                                                                                       tcgaccTGCTCAACGCCGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCAGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCtt              >  3:304542/5‑149 (MQ=255)
                                                                                                                                       aTCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGAACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCtt              >  1:289237/1‑149 (MQ=255)
                                                                                                                                           ccTGCTCAACGCCGAGCAGGTCGAGCAGTGGCCGTGTCGCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCcggcgg                                     <  5:293193/122‑1 (MQ=255)
                                                                                                                                           ccTGCTCAACGCCGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCcggcgg                                     >  6:293193/1‑122 (MQ=255)
                                                                                                                                            cTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCTTGGCCt         >  5:215637/1‑149 (MQ=255)
                                                                                                                                              gCTCGGCGTTGAGCAGGTCGAGCTGTGGCCGTGTCCGTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGGCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGACGGTTGCGCCTTGGCCTcc       >  8:147462/1‑149 (MQ=255)
                                                                                                                                                   gCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCTTGGCCTCCGCCAg  <  6:339551/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGATCTCCACACCTTGCAGGTTGCGCCCCTCAAGGCCAAAGGCCGTGCGCAGGGCGCCGGTGTTGTACAGGCAGTCGACATCCAGGCGCAGCAGATCGGAGCGGGTGCCGGCAAGTTCGAAGCGTGACTGGTAGATCTCCTGCTCGGCGTTGAGCAGGTCGAGCAGTGGCCGTGTCCCTAGGTCCAGGTACTGACGCCCGTAGAGCATGCGTGCCTCCTGGATGCTGACCTGGCGGGCTTCCAGCGAAGTCAGCCGGCGGCCCAGCCCGGAGGTTTGCGCCTTGGCCTCCGCCAG  >  NC_002947/925197‑925491

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: