Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,779 C→G 8.3% Q134H (CAG→CAC)  ccmH ← cytochrome c heme lyase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,911,7790CG 8.3% 281.4 / 10.5 93Q134H (CAG→CAC) ccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base G (2/6);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911632‑4911925
                                                                                                                                                     |                                                                                                                                                   
cggcCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCt                                                                                                                                                    >  7:253415/1‑149 (MQ=255)
 ggcCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACACATCGGCGCCCTg                                                                                                                                                   >  2:363721/1‑149 (MQ=255)
   ccACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑a                                                                                                                                                >  4:333755/1‑149 (MQ=255)
   ccACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑a                                                                                                                                                <  6:338404/149‑1 (MQ=255)
     aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑Ac                                                                                                                                               >  1:193645/1‑148 (MQ=255)
     aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑Ac                                                                                                                                               <  2:193645/148‑1 (MQ=255)
     aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCAATAATTAAATCATGGGTATTTTATTTTTCAAGCATTTTG‑CCGAGACGATCCCTTTCT‑ACGGAAGAAAAATAAACGCCGCCATGC‑ACg                                                                                                                                              >  6:276691/1‑149 (MQ=255)
           gATCGGGATCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTa                                                                                                                                        <  8:367004/149‑1 (MQ=255)
                    cagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGgcgc                                                                                                                               <  2:352344/149‑1 (MQ=255)
                    cagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTATTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGgcgc                                                                                                                               <  1:140065/149‑1 (MQ=255)
                       caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCga                                                                                                                            <  2:80211/149‑2 (MQ=255)
                       caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCgg                                                                                                                            >  8:88043/1‑149 (MQ=255)
                       caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCAAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGAACGGCGCCgg                                                                                                                            >  3:92525/1‑149 (MQ=255)
                        aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCgat                                                                                                                           <  1:296566/149‑3 (MQ=255)
                        aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGc                                                                                                                           <  1:375584/149‑1 (MQ=255)
                        aaaaaGATGATGGCAGCGAGCAGCAGAAGGCCCGCTCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGc                                                                                                                           <  4:92525/149‑1 (MQ=255)
                          aaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGCGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGt                                                                                                                         <  3:370183/149‑1 (MQ=255)
                           aaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTa                                                                                                                        <  7:152825/149‑1 (MQ=255)
                            aGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCATAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                       <  6:342168/149‑1 (MQ=255)
                            aGCTGAGGGCAGCGAGCACCATAAGGCCCTCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                       <  8:253415/149‑1 (MQ=255)
                             gCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACg                                                                                                                      <  8:332179/149‑1 (MQ=255)
                                gAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGGTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGAt                                                                                                                    >  4:406130/1‑148 (MQ=255)
                                             agcagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGcacca                                                                                                      >  7:91440/1‑149 (MQ=255)
                                               cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGcaccac                                                                                                     >  1:352493/1‑148 (MQ=255)
                                                agAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACg                                                                                                    >  8:354806/1‑148 (MQ=255)
                                                 taaGGGCCGGGCTTAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  5:276691/148‑1 (MQ=255)
                                                 gAAGGCCCGCGCTAAGCCAGAATTCAATCATGTGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  6:403787/149‑1 (MQ=255)
                                                 gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  >  7:360569/1‑149 (MQ=255)
                                                 gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  1:410822/149‑1 (MQ=255)
                                                       ccGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAcc                                                                                            >  2:389704/1‑149 (MQ=255)
                                                        cgcgCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCg                                                                                           >  7:328404/1‑149 (MQ=255)
                                                            cTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGTACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCa                                                                                       <  8:140606/149‑1 (MQ=255)
                                                               aGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCa                                                                                    <  7:354806/149‑1 (MQ=255)
                                                                 ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCTAACGGCAGACAAATCATCGGCGCCCTGC‑ACGGCCTTACAACGCCCCCGCCGAACGATAACCGCCACCACAACTAAGCCACCGGCCAGCa                                                                                    >  1:342095/1‑148 (MQ=255)
                                                                 ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑CCGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCCGc                                                                                     >  6:214954/1‑146 (MQ=255)
                                                                 ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAg                                                                                  <  1:236171/149‑1 (MQ=255)
                                                                   gacctacAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAg                                                                                  <  1:379314/141‑1 (MQ=255)
                                                                    gAATTCAATCATGTGTCTGTTCTTTTTCCAGCAGTTTGCTCAATACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGAt                                                                                <  1:363721/149‑1 (MQ=255)
                                                                    gAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                               >  8:33070/1‑149 (MQ=255)
                                                                         cAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                          >  3:13438/1‑149 (MQ=255)
                                                                          aaTCTTGGGTCTTTTCTTTTTGCTGCAGTTTG‑GCGTGACGCTCGCGTTCT‑TCGGCAGAGAACTCATCGGCGCCGTGC‑TCGGCCGTAGGAAGGCGCCGGCGTACGATCACCGCCAGCACCTCGAAGCCAGCTGGCAGCAAGATCCCCggg                                                                         <  2:342095/149‑1 (MQ=255)
                                                                          aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCggg                                                                         >  7:247166/1‑149 (MQ=255)
                                                                          aaTCATGGGTCTGTTCTTTTTCCAGCAGGTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                          <  3:406130/148‑1 (MQ=255)
                                                                           aTCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                        <  3:32830/149‑1 (MQ=255)
                                                                               tGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCaaa                                                                    <  7:405231/149‑1 (MQ=255)
                                                                                  gTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACca                                                                 >  7:323106/1‑149 (MQ=255)
                                                                                     agTTCTTTTGCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACACCTCCTCGGCTCCCTGC‑ACGGCCGTATCACGGCGCCGCCGGACGCACACAGACCGCACCACGAAGCCACCCGCCAGCAAGATCACCGGGCCACACCCcag                                                              >  6:434032/2‑149 (MQ=255)
                                                                                       ttCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                            <  8:328404/149‑1 (MQ=255)
                                                                                       ttCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              >  6:71061/1‑147 (MQ=255)
                                                                                         ccttttCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  5:266010/3‑149 (MQ=255)
                                                                                          tttttCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGCAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  6:305978/1‑148 (MQ=255)
                                                                                               ccAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGt                                                     <  3:347058/148‑1 (MQ=255)
                                                                                                cagcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                         >  6:398838/1‑111 (MQ=255)
                                                                                                cagcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                         <  5:398838/111‑1 (MQ=255)
                                                                                                 agcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                >  6:201654/1‑135 (MQ=255)
                                                                                                 agcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                <  5:201654/135‑1 (MQ=255)
                                                                                                 agcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                <  4:5424/135‑1 (MQ=255)
                                                                                                  gcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              >  8:136052/1‑136 (MQ=255)
                                                                                                  gcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              <  7:136052/136‑1 (MQ=255)
                                                                                                  gcagTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGc                                                 >  7:164499/1‑149 (MQ=255)
                                                                                                     gTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCg                                              >  5:297479/1‑149 (MQ=255)
                                                                                                      tttG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACgtcaca                                             >  2:155358/1‑144 (MQ=255)
                                                                                                        tG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                            <  2:352493/148‑1 (MQ=255)
                                                                                                         g‑gCGAGACGCTAGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTg                                           <  8:360569/148‑1 (MQ=255)
                                                                                                           gCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGa                                          <  7:33070/148‑1 (MQ=255)
                                                                                                            cGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAgc                                        <  3:333755/149‑1 (MQ=255)
                                                                                                               gaCGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGcc                                     <  5:71061/149‑1 (MQ=255)
                                                                                                                  gCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  1:189164/1‑149 (MQ=255)
                                                                                                                  gCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  7:293999/1‑149 (MQ=255)
                                                                                                                  gCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  4:248364/1‑149 (MQ=255)
                                                                                                                  gCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  <  3:248364/149‑1 (MQ=255)
                                                                                                                        gTTCT‑TAGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATAACCGCCAGCACCACGAAGCCAACGGCCAGCAAGATACCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAg                            >  4:193800/1‑149 (MQ=255)
                                                                                                                            t‑tcGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  4:133191/1‑149 (MQ=255)
                                                                                                                            t‑tcGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  7:216794/1‑149 (MQ=255)
                                                                                                                            t‑tcGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  4:77466/1‑149 (MQ=255)
                                                                                                                            t‑tcGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  6:86786/1‑149 (MQ=255)
                                                                                                                                ggCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACg                      >  3:105765/1‑148 (MQ=255)
                                                                                                                                 gCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                    <  7:88043/149‑1 (MQ=255)
                                                                                                                                 gCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                    <  3:133191/149‑1 (MQ=255)
                                                                                                                                  cAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCCCCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAg                   >  3:21307/1‑149 (MQ=255)
                                                                                                                                      cAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGc               >  2:85430/1‑149 (MQ=255)
                                                                                                                                       cgtacgccGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                 >  7:317521/8‑146 (MQ=38)
                                                                                                                                       cgtacgccGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                 <  8:317521/139‑1 (MQ=38)
                                                                                                                                         cTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa             >  5:423993/1‑148 (MQ=255)
                                                                                                                                          tCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa             >  7:305505/1‑147 (MQ=255)
                                                                                                                                            aTCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg          <  1:389704/148‑1 (MQ=255)
                                                                                                                                            aTCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg          <  2:340331/148‑1 (MQ=255)
                                                                                                                                             tCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg        <  5:86786/149‑1 (MQ=255)
                                                                                                                                             ccGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg        >  8:102598/2‑149 (MQ=255)
                                                                                                                                              cGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                     >  3:77493/1‑135 (MQ=255)
                                                                                                                                              cGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                     <  4:77493/135‑1 (MQ=255)
                                                                                                                                               ggCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt      >  6:207971/1‑149 (MQ=255)
                                                                                                                                               ggCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGAGCACGAAGTCGCCATAGCgg       <  6:297479/148‑1 (MQ=255)
                                                                                                                                                   cccTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc  <  8:164499/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGCCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTG‑GCGAGACGCTCGCGTTCT‑TCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911632‑4911925

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