Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,756,942 A→G 5.0% H72R (CAC→CGC)  ompQ → outer membrane pyoverdine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,756,9420AG 5.0% 238.1 / 2.1 80H72R (CAC→CGC) ompQouter membrane pyoverdine efflux protein
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (2/2);  total (42/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 = 6.77e-01

ACGCCGCCAACTATGGGCATAGCCCCGCCCGCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTCGATGC  >  NC_002947/4756794‑4757089
                                                                                                                                                    |                                                                                                                                                   
aCGCCGACAACAATCGGCATAGCCCCGCCCGCCGGTAGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCa                                                                                                                                                     <  8:301821/149‑1 (MQ=255)
    cgccAACTATGGGCATAGCCCCGCCCGCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGCGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGAc                                                                                                                                                 >  3:93033/1‑149 (MQ=255)
                    aGCCCCGCCCGCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGc                                                                                                                                      >  5:222064/1‑144 (MQ=255)
                    aGCCCCGCCCGCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGc                                                                                                                                      <  6:222064/144‑1 (MQ=255)
                       cccgcccgCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGccc                                                                                                                              >  2:128937/1‑149 (MQ=255)
                        ccgcccgcCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCCACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCCGCCCg                                                                                                                             >  7:234466/1‑149 (MQ=255)
                        ccgcccgcCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCTAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCGCGGCCAGGTCGTGGCTgttgtttc                                                                                                                             <  2:267489/149‑8 (MQ=255)
                         cgcccgccGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCgg                                                                                                                            <  1:128937/149‑1 (MQ=255)
                         cgcccgccGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCgg                                                                                                                            <  7:228786/149‑1 (MQ=255)
                         cgcccgccGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCgg                                                                                                                            <  7:234777/149‑1 (MQ=255)
                             cgccgGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGcc                                                                                                                               <  8:134242/142‑1 (MQ=255)
                             cgccgGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGcc                                                                                                                               >  7:134242/1‑142 (MQ=255)
                                     ggCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGccaggc                                                                                                                 <  6:71143/148‑1 (MQ=255)
                                      gCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGccaggccc                                                                                                               >  8:82648/1‑149 (MQ=255)
                                       cAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGccaggccca                                                                                                              >  5:216690/1‑149 (MQ=255)
                                        aaaGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGccaggcccag                                                                                                             >  6:332314/1‑149 (MQ=255)
                                           gCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCgg                                                                                                                            >  3:125071/1‑131 (MQ=255)
                                           gCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCgg                                                                                                                            <  4:125071/131‑1 (MQ=255)
                                              cTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCgg                                                                                                       >  8:154628/1‑149 (MQ=255)
                                                aCTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCggg                                                                                                      >  5:187469/1‑148 (MQ=255)
                                                 cTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCggg                                                                                                      <  5:301938/147‑1 (MQ=255)
                                                 cTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGc                                                                                                     <  7:82648/148‑1 (MQ=255)
                                                 cTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCa                                                                                                    <  7:154628/149‑1 (MQ=255)
                                                 cTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGCGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCggg                                                                                                      <  4:93033/147‑1 (MQ=255)
                                                     gAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTaa                                                                                                >  3:229309/1‑149 (MQ=255)
                                                     gAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTaa                                                                                                <  8:32117/149‑1 (MQ=255)
                                                     gAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTaa                                                                                                >  1:72812/1‑149 (MQ=255)
                                                      aGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAAt                                                                                               >  8:309205/1‑149 (MQ=255)
                                                       gCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATc                                                                                              <  7:309205/149‑1 (MQ=255)
                                                       gCAAGGCACAGGCACTGCCGGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATc                                                                                              >  8:104510/1‑149 (MQ=255)
                                                            gcacaggcacTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCgg                                                                                         >  8:109193/1‑149 (MQ=255)
                                                            gcacaggcacTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCgg                                                                                         >  8:166737/1‑149 (MQ=255)
                                                                  gcacTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCCCGACCTGGCCGCGGCTGCAGCCCGGGTACGc                                                                                                                      <  4:121215/114‑1 (MQ=255)
                                                                  gcacTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGc                                                                                                                      >  3:121215/1‑114 (MQ=255)
                                                                  gcacTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGcgcc                                                                                   >  6:115563/1‑149 (MQ=255)
                                                                       gCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCgctgc                                                                              >  4:234283/1‑149 (MQ=255)
                                                                        ccTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCgctgct                                                                             >  1:199285/1‑149 (MQ=255)
                                                                            aTCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCtgctgc                                                                           <  4:48048/147‑1 (MQ=255)
                                                                            aTCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCtgctgc                                                                           <  3:234283/147‑1 (MQ=255)
                                                                            aTCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCt                                                                         <  8:35646/149‑1 (MQ=255)
                                                                               aaCGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGa                                                                       <  5:115563/148‑1 (MQ=255)
                                                                               aaCGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAg                                                                      <  7:166737/149‑1 (MQ=255)
                                                                               aaCGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAg                                                                      <  7:109193/149‑1 (MQ=255)
                                                                                aCGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAgg                                                                     >  3:169641/1‑149 (MQ=255)
                                                                                  gCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTg                                                                   >  6:12856/1‑149 (MQ=255)
                                                                                        ggAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCa                                                                 <  1:287233/145‑1 (MQ=255)
                                                                                        ggAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCa                                                                 >  2:287233/1‑145 (MQ=255)
                                                                                         gAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCg                                                            >  3:8133/1‑149 (MQ=255)
                                                                                           aGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGc                                                          >  2:180120/1‑149 (MQ=255)
                                                                                            gcccTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGcc                                                         <  5:332314/147‑1 (MQ=255)
                                                                                                ttCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGa                                                     >  7:5420/1‑149 (MQ=255)
                                                                                                  cGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCtgctgc                                                                           >  4:114281/1‑125 (MQ=255)
                                                                                                  cGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCtgctgc                                                                           <  3:114281/125‑1 (MQ=255)
                                                                                                    ccAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGGAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGcc                                                 >  4:122494/1‑149 (MQ=255)
                                                                                                     cAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCa                                                <  4:229309/149‑1 (MQ=255)
                                                                                                      aGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCGTGGATGCCAg                                               >  8:91783/1‑149 (MQ=255)
                                                                                                      aGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAg                                               <  6:187469/149‑1 (MQ=255)
                                                                                                         ccAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCAAGGCGCGGGCATAAAACCCCGCCGCGCCGCTGCTCCCTGGGGTGCAGTTCGGCCTGGATGCCAgccc                                            >  5:247698/1‑148 (MQ=255)
                                                                                                          cAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAgccgt                                           <  8:5420/149‑2 (MQ=255)
                                                                                                              gCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAgc                                       <  7:104510/149‑1 (MQ=255)
                                                                                                                ttGACCAGTTGATCGAGATTGCCCGAAACAACAGCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGcc                                    >  7:344568/1‑149 (MQ=255)
                                                                                                                   aCCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCct                                   <  4:169641/148‑1 (MQ=255)
                                                                                                                    ccAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCctgc                                 <  6:216690/149‑1 (MQ=255)
                                                                                                                       gTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGATGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCtgctgc                              >  2:340063/1‑149 (MQ=255)
                                                                                                                               agaTTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      <  4:8133/149‑1 (MQ=255)
                                                                                                                               agaTTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      >  7:118187/1‑149 (MQ=255)
                                                                                                                               agaTTGCCCGAAACAACAGCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAgg                     <  8:344568/149‑1 (MQ=255)
                                                                                                                                     cccGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCgcgg                          <  3:880/139‑1 (MQ=255)
                                                                                                                                     cccGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCgcgg                          >  4:880/1‑139 (MQ=255)
                                                                                                                                             aacaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      <  5:111545/135‑1 (MQ=255)
                                                                                                                                             aacaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      >  6:111545/1‑135 (MQ=255)
                                                                                                                                             aacaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCt        >  3:356446/1‑149 (MQ=255)
                                                                                                                                             aacaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCt        >  2:200352/1‑149 (MQ=255)
                                                                                                                                              gcaGCCGCGGCCAGGTCGTGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      <  3:105283/133‑1 (MQ=255)
                                                                                                                                              gcaGCCGCGGCCAGGTCGTGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAg                      >  4:105283/2‑134 (MQ=255)
                                                                                                                                              gcaGCCGCGGCCAGGTCGTGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTc       >  7:102481/2‑149 (MQ=255)
                                                                                                                                              acaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTc       <  3:175949/149‑1 (MQ=255)
                                                                                                                                              acaGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTc       >  4:175949/1‑149 (MQ=255)
                                                                                                                                                 gCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTCGAt    <  2:72812/149‑1 (MQ=255)
                                                                                                                                                  ccACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAgc                                       >  6:126556/1‑113 (MQ=255)
                                                                                                                                                  ccACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAgc                                       <  5:126556/113‑1 (MQ=255)
                                                                                                                                                    aCGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTCGATGc  >  4:288439/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGCCGCCAACTATGGGCATAGCCCCGCCCGCCGGTTGGCAAAGCCCTACTGCGAGCAAGGCACAGGCACTGCCTGATCAACGCTGGTGGAAGGCCTTCGCCAGCCCAGAGCTTGACCAGTTGATCGAGATTGCCCGAAACAACAGCCACGACCTGGCCGCGGCTGCAGCCCGGGTACGCCAGGCCCAGGCGCGGGCAGTAATCGCCGGCGCGCCGCTGCTGCCTGAGGTGCAGTTCGGCCTGGATGCCAGCCGCCAGCGCCTGCTGCGCGGCGAGGGCAACGATCAGCTCGATGC  >  NC_002947/4756794‑4757089

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: