Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,782:1 +T 7.8% coding (399/477 nt) 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,7821.T 7.8% 352.0 / 21.1 90coding (399/477 nt)ccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base . (42/41);  new base T (2/5);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

CCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911635‑4911925
                                                                                                                                                      |                                                                                                                                               
ccACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCa                                                                                                                                                >  4:333755/1‑149 (MQ=255)
ccACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCa                                                                                                                                                <  6:338404/149‑1 (MQ=255)
  aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCAc                                                                                                                                               >  1:193645/1‑148 (MQ=255)
  aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCAc                                                                                                                                               <  2:193645/148‑1 (MQ=255)
  aCGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCAATAATTAAATCATGGGTATTTTATTTTTCAAGCATTTTGCCGAGACGATCCCTTTCTACGGAAGAAAAATAAACGCCGCCATGCACg                                                                                                                                              >  6:276691/1‑149 (MQ=255)
        gATCGGGATCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTa                                                                                                                                        <  8:367004/149‑1 (MQ=255)
                 cagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGgcgc                                                                                                                               <  2:352344/149‑1 (MQ=255)
                 cagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTATTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGgcgc                                                                                                                               <  1:140065/149‑1 (MQ=255)
                    caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCga                                                                                                                            <  2:80211/149‑2 (MQ=255)
                    caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCgg                                                                                                                            >  8:88043/1‑149 (MQ=255)
                    caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCAAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGAACGGCGCCgg                                                                                                                            >  3:92525/1‑149 (MQ=255)
                     aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCgat                                                                                                                           <  1:296566/149‑3 (MQ=255)
                     aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGc                                                                                                                           <  1:375584/149‑1 (MQ=255)
                     aaaaaGATGATGGCAGCGAGCAGCAGAAGGCCCGCTCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGc                                                                                                                           <  4:92525/149‑1 (MQ=255)
                       aaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGCGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGt                                                                                                                         <  3:370183/149‑1 (MQ=255)
                        aaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTa                                                                                                                        <  7:152825/149‑1 (MQ=255)
                         aGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCATAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                       <  6:342168/149‑1 (MQ=255)
                         aGCTGAGGGCAGCGAGCACCATAAGGCCCTCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                       <  8:253415/149‑1 (MQ=255)
                          gCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACg                                                                                                                      <  8:332179/149‑1 (MQ=255)
                             gAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGGTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGAt                                                                                                                    >  4:406130/1‑148 (MQ=255)
                                          agcagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGcacca                                                                                                      >  7:91440/1‑149 (MQ=255)
                                            cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGcaccac                                                                                                     >  1:352493/1‑148 (MQ=255)
                                             agAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACg                                                                                                    >  8:354806/1‑148 (MQ=255)
                                              taaGGGCCGGGCTTAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  5:276691/148‑1 (MQ=255)
                                              gAAGGCCCGCGCTAAGCCAGAATTCAATCATGTGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  6:403787/149‑1 (MQ=255)
                                              gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  >  7:360569/1‑149 (MQ=255)
                                              gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  1:410822/149‑1 (MQ=255)
                                                    ccGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAcc                                                                                            >  2:389704/1‑149 (MQ=255)
                                                     cgcgCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCg                                                                                           >  7:328404/1‑149 (MQ=255)
                                                         cTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGTACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCa                                                                                       <  8:140606/149‑1 (MQ=255)
                                                            aGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCa                                                                                    <  7:354806/149‑1 (MQ=255)
                                                              ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTAACGGCAGACAAATCATCGGCGCCCTGCACGGCCTTACAACGCCCCCGCCGAACGATAACCGCCACCACAACTAAGCCACCGGCCAGCa                                                                                    >  1:342095/1‑148 (MQ=255)
                                                              ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCCCGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCCGc                                                                                     >  6:214954/1‑146 (MQ=255)
                                                              ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAg                                                                                  <  1:236171/149‑1 (MQ=255)
                                                                gacctacAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAg                                                                                  <  1:379314/141‑1 (MQ=255)
                                                                 gAATTCAATCATGTGTCTGTTCTTTTTCCAGCAGTTTGCTCAATACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGAt                                                                                <  1:363721/149‑1 (MQ=255)
                                                                 gAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                               >  8:33070/1‑149 (MQ=255)
                                                                      cAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                          >  3:13438/1‑149 (MQ=255)
                                                                       aaTCTTGGGTCTTTTCTTTTTGCTGCAGTTTGGCGTGACGCTCGCGTTCTTCGGCAGAGAACTCATCGGCGCCGTGCTCGGCCGTAGGAAGGCGCCGGCGTACGATCACCGCCAGCACCTCGAAGCCAGCTGGCAGCAAGATCCCCggg                                                                         <  2:342095/149‑1 (MQ=255)
                                                                       aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCggg                                                                         >  7:247166/1‑149 (MQ=255)
                                                                       aaTCATGGGTCTGTTCTTTTTCCAGCAGGTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                          <  3:406130/148‑1 (MQ=255)
                                                                        aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                        <  3:32830/149‑1 (MQ=255)
                                                                            tGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCaaa                                                                    <  7:405231/149‑1 (MQ=255)
                                                                               gTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACca                                                                 >  7:323106/1‑149 (MQ=255)
                                                                                  agTTCTTTTGCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACACCTCCTCGGCTCCCTGCACGGCCGTATCACGGCGCCGCCGGACGCACACAGACCGCACCACGAAGCCACCCGCCAGCAAGATCACCGGGCCACACCCcag                                                              >  6:434032/2‑149 (MQ=255)
                                                                                    ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                            <  8:328404/149‑1 (MQ=255)
                                                                                    ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              >  6:71061/1‑147 (MQ=255)
                                                                                      ccttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  5:266010/3‑149 (MQ=255)
                                                                                       tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGCAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  6:305978/1‑148 (MQ=255)
                                                                                            ccAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGt                                                     <  3:347058/148‑1 (MQ=255)
                                                                                             cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                         >  6:398838/1‑111 (MQ=255)
                                                                                             cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                         <  5:398838/111‑1 (MQ=255)
                                                                                              agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                >  6:201654/1‑135 (MQ=255)
                                                                                              agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                <  5:201654/135‑1 (MQ=255)
                                                                                              agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                <  4:5424/135‑1 (MQ=255)
                                                                                               gcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              >  8:136052/1‑136 (MQ=255)
                                                                                               gcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              <  7:136052/136‑1 (MQ=255)
                                                                                               gcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGc                                                 >  7:164499/1‑149 (MQ=255)
                                                                                                  gTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCg                                              >  5:297479/1‑149 (MQ=255)
                                                                                                   tttGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACgtcaca                                             >  2:155358/1‑144 (MQ=255)
                                                                                                     tGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                            <  2:352493/148‑1 (MQ=255)
                                                                                                      ggCGAGACGCTAGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTg                                           <  8:360569/148‑1 (MQ=255)
                                                                                                        gCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGa                                          <  7:33070/148‑1 (MQ=255)
                                                                                                         cGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAgc                                        <  3:333755/149‑1 (MQ=255)
                                                                                                            gaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGcc                                     <  5:71061/149‑1 (MQ=255)
                                                                                                               gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  1:189164/1‑149 (MQ=255)
                                                                                                               gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  <  3:248364/149‑1 (MQ=255)
                                                                                                               gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  7:293999/1‑149 (MQ=255)
                                                                                                               gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  4:248364/1‑149 (MQ=255)
                                                                                                                     gTTCTTAGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATAACCGCCAGCACCACGAAGCCAACGGCCAGCAAGATACCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAg                            >  4:193800/1‑149 (MQ=255)
                                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  4:133191/1‑149 (MQ=255)
                                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  7:216794/1‑149 (MQ=255)
                                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  6:86786/1‑149 (MQ=255)
                                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                        >  4:77466/1‑149 (MQ=255)
                                                                                                                             ggCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACg                      >  3:105765/1‑148 (MQ=255)
                                                                                                                              gCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                    <  7:88043/149‑1 (MQ=255)
                                                                                                                              gCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                    <  3:133191/149‑1 (MQ=255)
                                                                                                                               cAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCCCCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAg                   >  3:21307/1‑149 (MQ=255)
                                                                                                                                   cAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGc               >  2:85430/1‑149 (MQ=255)
                                                                                                                                    cgtacgccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                 >  7:317521/8‑146 (MQ=38)
                                                                                                                                    cgtacgccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                 <  8:317521/139‑1 (MQ=38)
                                                                                                                                      cTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa             >  5:423993/1‑148 (MQ=255)
                                                                                                                                       tCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa             >  7:305505/1‑147 (MQ=255)
                                                                                                                                         aTCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg          <  1:389704/148‑1 (MQ=255)
                                                                                                                                         aTCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg          <  2:340331/148‑1 (MQ=255)
                                                                                                                                          tCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg        <  5:86786/149‑1 (MQ=255)
                                                                                                                                          ccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg        >  8:102598/2‑149 (MQ=255)
                                                                                                                                           cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                     <  4:77493/135‑1 (MQ=255)
                                                                                                                                           cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                     >  3:77493/1‑135 (MQ=255)
                                                                                                                                            ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt      >  6:207971/1‑149 (MQ=255)
                                                                                                                                            ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGAGCACGAAGTCGCCATAGCgg       <  6:297479/148‑1 (MQ=255)
                                                                                                                                                cccTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc  <  8:164499/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
CCACGCAGGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911635‑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: