Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,017 G→T 5.2% L744M (CTG→ATG) ‡ cadA‑II ← cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947684,0170GT 5.2% 253.5 / 3.1 77L744M (CTG→ATG) ‡cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base T (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

CTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGTCGGCCGCG  >  NC_002947/683869‑684162
                                                                                                                                                    |                                                                                                                                                  
ctgGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAg                                                                                                                                                    <  4:43270/149‑1 (MQ=255)
   gtggGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAgg                                                                                                                                                   <  5:80248/147‑1 (MQ=255)
   gtggGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGtt                                                                                                                                                 >  8:117048/1‑149 (MQ=255)
     gggTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCt                                                                                                                                               >  8:20254/1‑149 (MQ=255)
      ggTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTg                                                                                                                                              >  2:169490/1‑149 (MQ=255)
        tttCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGc                                                                                                                                            >  7:300290/1‑149 (MQ=255)
          tCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACGGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACCTGTTCTgggc                                                                                                                                           <  5:116229/148‑3 (MQ=255)
            cACGTTTTAAGCCACAGTGCATTGCTGAAGACGCTGACACTGGACACGGCCATGGCGGCTCCTGCCAGTACCGCGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGGTCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGAt                                                                                                                                        >  5:268481/1‑149 (MQ=255)
                      gCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAg                                                                                                                             <  6:121883/149‑1 (MQ=255)
                        cacaGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCAAGAACAGGTTCTGGCGGATCTTCGCGTAGGt                                                                                                                           >  8:175031/1‑149 (MQ=255)
                         acaGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTc                                                                                                                          <  2:198606/149‑1 (MQ=255)
                             tGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGc                                                                                                                      >  6:208000/1‑149 (MQ=255)
                                aTTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTTTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGc                                                                                                                   <  3:154227/149‑1 (MQ=255)
                                  tGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCt                                                                                                                  >  5:382445/1‑148 (MQ=255)
                                      gacgacGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGAttt                                                                                                             >  5:330904/1‑149 (MQ=255)
                                        cgacgCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGGTCAGGTTGTCGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTc                                                                                                            <  6:268481/148‑1 (MQ=255)
                                                acacTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCgcggc                                                                                                   <  2:369788/149‑1 (MQ=255)
                                                   cTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCggg                                                                                                >  3:29620/1‑149 (MQ=255)
                                                            gCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGgcg                                                                                       >  7:115730/1‑149 (MQ=255)
                                                             ccATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTTGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGgcgc                                                                                      <  2:92365/149‑1 (MQ=255)
                                                               aTGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCgg                                                                                    >  4:312432/1‑149 (MQ=255)
                                                               aTGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCgg                                                                                    <  8:300290/149‑1 (MQ=255)
                                                                 ggcggcTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGt                                                                                  >  5:357732/1‑149 (MQ=255)
                                                                  gcggcTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGt                                                                                  <  7:117048/148‑1 (MQ=255)
                                                                   cggcTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCg                                                                                >  7:159653/1‑149 (MQ=255)
                                                                      cTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTTTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCa                                                                             <  6:357732/149‑1 (MQ=255)
                                                                      cTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCa                                                                             <  6:382445/149‑1 (MQ=255)
                                                                         cTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTc                                                                                                                                   <  6:26125/92‑1 (MQ=255)
                                                                         cTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTc                                                                                                                                   >  5:26125/1‑92 (MQ=255)
                                                                         cTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTc                                                                                                            <  8:78283/115‑1 (MQ=255)
                                                                         cTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTc                                                                                                            >  7:78283/1‑115 (MQ=255)
                                                                            ccAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCt                                                                                                                                     >  5:306116/1‑87 (MQ=255)
                                                                            ccAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCt                                                                                                                                     <  6:306116/87‑1 (MQ=255)
                                                                                tACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAggg                                                                   >  6:348669/1‑149 (MQ=255)
                                                                                   cGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGCACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGa                                                                >  6:53360/1‑149 (MQ=255)
                                                                                      gTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGc                                                             >  4:16582/1‑149 (MQ=255)
                                                                                         gAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCcg                                                           <  6:330904/148‑1 (MQ=255)
                                                                                         gAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCACCCTGATGcgt                                                           <  7:20254/148‑3 (MQ=255)
                                                                                          aGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCACCCTGATGcgtg                                                          <  4:261920/148‑4 (MQ=255)
                                                                                            gTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCcggcgg                                                       >  3:322257/1‑149 (MQ=255)
                                                                                                 ccAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTg                                                   >  7:305420/1‑148 (MQ=255)
                                                                                                       cggcgAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGcc                                            <  2:177391/149‑1 (MQ=255)
                                                                                                         gcgAGCGGAATGCCGATCAGGTTGAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCAc                                                                            >  4:217912/1‑114 (MQ=255)
                                                                                                         gcgAGCGGAATGCCGATCAGGTTGAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCAc                                                                            <  3:217912/114‑1 (MQ=255)
                                                                                                            aGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTc                                       <  1:169490/149‑1 (MQ=255)
                                                                                                              cGGAATGCCGATCAGGTTGTTGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCgg                                     >  4:134344/1‑149 (MQ=255)
                                                                                                              cGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCgg                                     <  3:16582/149‑1 (MQ=255)
                                                                                                              cGGAATGCCGATCAGATTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCa                                                                                           <  2:125373/95‑1 (MQ=255)
                                                                                                                gAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCggg                                                                                                >  4:245332/1‑88 (MQ=255)
                                                                                                                gAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCggg                                                                                                <  3:245332/88‑1 (MQ=255)
                                                                                                                  aTGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCcggcg                                                        >  1:114977/1‑126 (MQ=255)
                                                                                                                  aTGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCcggcg                                                        <  2:114977/126‑1 (MQ=255)
                                                                                                                    gCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTCCCGc                               >  8:341718/1‑149 (MQ=255)
                                                                                                                     ccGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGc                                                  <  5:217223/129‑1 (MQ=255)
                                                                                                                     ccGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGc                                                  >  6:217223/1‑129 (MQ=255)
                                                                                                                      cGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCa                             >  6:42223/1‑149 (MQ=255)
                                                                                                                      cGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCa                             >  7:1647/1‑149 (MQ=255)
                                                                                                                       gATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCAc                            >  8:285863/1‑149 (MQ=255)
                                                                                                                       gATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCAc                            >  6:17891/1‑149 (MQ=255)
                                                                                                                        aTCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCAcc                           <  8:1647/149‑1 (MQ=255)
                                                                                                                            ggTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGAGCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATg                       <  7:166080/149‑1 (MQ=255)
                                                                                                                                    aTGAAGGCCCAGAACAGGTTCTGGCGGATCGTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGc                >  1:33897/1‑149 (MQ=255)
                                                                                                                                    aTGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCAt                        >  8:112558/1‑140 (MQ=255)
                                                                                                                                    aTGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCAt                        <  7:112558/140‑1 (MQ=255)
                                                                                                                                    aTGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGcc               <  3:312432/149‑1 (MQ=255)
                                                                                                                                      gAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGa             >  2:246960/1‑149 (MQ=255)
                                                                                                                                       atccgccAGAACCTGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGc                                                             >  8:248897/6‑100 (MQ=255)
                                                                                                                                       atccgccAGAACCTGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGc                                                             <  7:248897/95‑1 (MQ=255)
                                                                                                                                       aaGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGa             <  5:208000/148‑1 (MQ=255)
                                                                                                                                       aaGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGAt            <  5:53360/149‑1 (MQ=255)
                                                                                                                                       aaGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGAt            >  8:311197/1‑149 (MQ=255)
                                                                                                                                       aaGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGGGATGCCGGCGGCCTGCATGGCCACGTCGGTACCACCACCCATGGCAATGCCGAt            <  4:29620/149‑1 (MQ=255)
                                                                                                                                         ccgccAGAACCTGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGt          >  8:128419/4‑149 (MQ=255)
                                                                                                                                           cccAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGTCg        >  7:176459/1‑149 (MQ=255)
                                                                                                                                           cccAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGTCg        >  6:83123/1‑149 (MQ=255)
                                                                                                                                            ccAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCgg                                     >  8:144577/1‑119 (MQ=255)
                                                                                                                                            ccAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCgg                                     <  7:144577/119‑1 (MQ=255)
                                                                                                                                                aaCAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGACCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGTCGGCcgc   <  8:159653/149‑1 (MQ=255)
                                                                                                                                                 aCAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATACCGATGTCGGCcgcg  <  5:348669/149‑1 (MQ=255)
                                                                                                                                                    ggTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCa                                           >  7:64141/1‑105 (MQ=255)
                                                                                                                                                    ggTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCa                                           <  8:64141/105‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGTGGGTTTCCACGTTTTAAGCCACAGTGCATTGCTGACGACGCTGACACTGGACAGGGCCATGGCGGCTCCTGCCAGTACCGGGTTGAGGTAACCCAGCGCGGCGAGCGGAATGCCGATCAGGTTGTAGATGAAGGCCCAGAACAGGTTCTGGCGGATCTTCGCGTAGGTCTTGCGGCTGATTTCCAGCGCGGCGGGCACCAGGCGCGGGTCGCCACGCATCAGGGTGATGCCGGCGGCCTGCATGGCCACGTCGGTACCGCCACCCATGGCAATGCCGATGTCGGCCGCG  >  NC_002947/683869‑684162

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: