Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,503,149 T→G 6.9% P33P (CCT→CCG) ‡ oprC → copper receptor OprC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,503,1490TG 6.9% 204.7 / 7.5 73P33P (CCT→CCG) ‡oprCcopper receptor OprC
Reads supporting (aligned to +/- strand):  ref base T (27/39);  new base G (2/3);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

AGCCAGGAGGCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGA  >  NC_002947/5503004‑5503297
                                                                                                                                                  |                                                                                                                                                          
aGCCAGGAGGCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                                         <  3:237038/149‑1 (MQ=255)
aGCCAGGAGGCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                                         <  5:405919/149‑1 (MQ=255)
 gCCAGGAGGCTTCGCGCTAACTCATGACTGGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                                         >  1:242238/1‑149 (MQ=255)
  ccAGGAGGCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                                       <  7:365004/149‑1 (MQ=255)
    aggaggCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATg                                                                                                                                                     <  4:255493/149‑1 (MQ=255)
     ggagcCTTCGCGCTAACTCCTGTCT‑GATCGACTGGAATCTTTATGTCCGGCTGCACCGCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCTCCTCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGa                                                                                                                                                    <  8:374151/149‑1 (MQ=255)
       aggCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgggc                                                                                                                                                  <  1:385720/149‑4 (MQ=255)
                    aCTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt‑g                                                                                                                                         <  2:87283/144‑1 (MQ=255)
                    aCTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt‑g                                                                                                                                         >  1:87283/1‑144 (MQ=255)
                      tCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCCCCCATGTTTTTGCCTTGTAACTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGAGCCCATGGCCCTGGCAGCAGGCCCTGCCCATGAAGACCATGT‑GCa‑‑cgacg                                                                                                                                >  6:215659/1‑149 (MQ=255)
                      tCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCa‑‑cgacg                                                                                                                                >  3:178858/1‑149 (MQ=255)
                        aTGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                             <  7:219300/137‑1 (MQ=255)
                        aTGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                             >  8:219300/1‑137 (MQ=255)
                        aTGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCa                                                                                                                              >  6:65969/1‑149 (MQ=255)
                                aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                                           <  7:205131/116‑1 (MQ=255)
                                aTCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                                           >  8:205131/1‑116 (MQ=255)
                                   gACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTg                                                                                                                     <  5:141534/148‑1 (MQ=255)
                                   gACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGa                                                                                                                    <  7:142671/149‑1 (MQ=255)
                                     cTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGc                                                                                                                  >  2:314241/1‑149 (MQ=255)
                                        gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCg                                                                                                               >  3:185562/1‑149 (MQ=255)
                                          ataaTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTAGCGCCCATGGCCCTGGCCCCCCGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGCCGCACCCGAGTTGAGCCCGAc                                                                                                             >  1:184419/4‑149 (MQ=255)
                                           tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACg                                                                                                            <  2:242238/149‑1 (MQ=255)
                                           tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACg                                                                                                            <  4:203506/149‑1 (MQ=255)
                                            cgccttGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGTCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACg                                                                                                            <  1:347771/143‑1 (MQ=255)
                                             aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGt                                                                                                          >  4:265831/1‑149 (MQ=255)
                                              ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGt                                                                                                          >  6:381652/1‑148 (MQ=255)
                                              ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCAGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGCCGGt                                                                                                          >  4:243151/1‑148 (MQ=255)
                                                aTGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                                        <  6:115949/148‑1 (MQ=255)
                                                      ggCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCAccgc                                                                                                 >  5:273116/1‑149 (MQ=255)
                                                       gCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCAccgcc                                                                                                >  4:325410/1‑149 (MQ=255)
                                                           cACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTg                                                                                             <  4:185562/148‑1 (MQ=255)
                                                             ccccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCa‑‑cgacg                                                                                                                                >  3:12372/1‑111 (MQ=255)
                                                             ccccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCa‑‑cgacg                                                                                                                                <  4:12372/111‑1 (MQ=255)
                                                             ccccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑Gcc                                                                                         >  5:103242/1‑149 (MQ=255)
                                                              cccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑Gccc                                                                                        <  5:205298/149‑1 (MQ=255)
                                                              cccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑Gccc                                                                                        >  4:306690/1‑149 (MQ=255)
                                                               ccTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑Gcccc                                                                                       >  8:31637/1‑149 (MQ=255)
                                                                cTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATCTATTAGGCGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑Gcc                                                                                         >  1:135096/1‑149 (MQ=255)
                                                                     tttGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTc                                                                                 >  7:157600/1‑149 (MQ=255)
                                                                        gCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTcgc                                                                               <  2:190652/148‑1 (MQ=255)
                                                                         ccTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTcgccg                                                                             <  5:65969/149‑1 (MQ=255)
                                                                         ccTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTcgccg                                                                             <  3:265831/149‑1 (MQ=255)
                                                                          cTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGCCGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTcgccgc                                                                            >  8:304438/1‑149 (MQ=255)
                                                                              tCCCTGCGTGGGACCATCGCCGCCCTGTGCGTCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGa                                                                       <  2:389427/148‑1 (MQ=255)
                                                                                    cGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACtgtg                                                                  <  6:273116/147‑1 (MQ=255)
                                                                                     gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTgg                                                                 <  8:365588/147‑1 (MQ=255)
                                                                                     gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTgg                                                                 >  7:365588/1‑147 (MQ=255)
                                                                                     gTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTc                                                               <  4:178858/149‑1 (MQ=255)
                                                                                      tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCAGGCCATGAAGCCCATGT‑GCA‑‑CGCCGCACCCGAGTTGAGCCCGCCGGTCATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCa                                                              >  3:25160/1‑149 (MQ=255)
                                                                                           ccATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAAc                                                         <  5:343743/149‑1 (MQ=255)
                                                                                           ccATCGCCGCCCTGTGCGGATCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGAACCCGAGTTGAGCCCGACTGTGATCAACGCCGTT‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAAc                                                         <  5:215659/149‑1 (MQ=255)
                                                                                                gccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGa                                                     >  2:250381/1‑148 (MQ=255)
                                                                                                 ccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGaaa                                                   <  3:323073/149‑1 (MQ=255)
                                                                                                     ccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGAcc                                               <  7:31637/149‑1 (MQ=255)
                                                                                                         tgCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCcg                                             <  5:48651/147‑1 (MQ=255)
                                                                                                            ggCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCgcca                                        <  5:105150/149‑1 (MQ=255)
                                                                                                            ggCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCgcca                                        <  3:243151/149‑1 (MQ=255)
                                                                                                              cTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCgccagc                                      >  3:125496/1‑149 (MQ=255)
                                                                                                              cTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCgccagc                                      <  5:381652/149‑1 (MQ=255)
                                                                                                                     ctCGCGCCTATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTg                                 <  2:184419/147‑1 (MQ=255)
                                                                                                                     ctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTGCGCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTg                                 <  1:314241/148‑1 (MQ=255)
                                                                                                                        gcgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGAGCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                            <  1:250381/149‑1 (MQ=255)
                                                                                                                         cgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGc                            >  2:87378/1‑148 (MQ=255)
                                                                                                                          gcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCa                          <  2:109452/149‑1 (MQ=255)
                                                                                                                           ttcATGGCCAGGGCCGGCGGCCAGGGCCATGGAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGaaa                                                   >  4:186484/3‑123 (MQ=38)
                                                                                                                           ttcATGGCCAGGGCCGGCGGCCAGGGCCATGGAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGaaa                                                   <  3:186484/121‑1 (MQ=38)
                                                                                                                            tcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGc                        >  4:15132/2‑149 (MQ=255)
                                                                                                                                  cccTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCAAAATCGCCGCAATGA‑TGTGGTCACCAACCCGAAAGACCCGCGCCACCCAGTGCCGCCCAGCGACgg                   >  1:334284/1‑149 (MQ=255)
                                                                                                                                       gccgccGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGa             <  8:155474/149‑1 (MQ=255)
                                                                                                                                             ggCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCt       <  2:135096/149‑1 (MQ=255)
                                                                                                                                              gCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTg      <  4:125496/149‑1 (MQ=255)
                                                                                                                                                ccTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc         <  3:8852/144‑1 (MQ=255)
                                                                                                                                                ccTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc         >  4:8852/1‑144 (MQ=255)
                                                                                                                                                  tGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGa  >  2:329767/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                          
AGCCAGGAGGCTTCGCGCTAACTCATGACT‑GATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTCCCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGT‑GCA‑‑CGACGCACCCGAGTTGAGCCCGACGGTGATCACCGCCGTG‑GCCCCGAGTTCGCCGC‑‑TGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTACCTGAAGA  >  NC_002947/5503004‑5503297

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