Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,795 C→G 7.4% R129P (CGT→CCT)  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,7950CG 7.4% 288.9 / 10.1 94R129P (CGT→CCT) ccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base G (2/5);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATC  >  NC_002947/4911652‑4911941
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                  <  1:410822/149‑1 (MQ=255)
                             gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                                  >  7:360569/1‑149 (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                                                                                                         <  5:398838/111‑1 (MQ=255)
                                                                            cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                                         >  6:398838/1‑111 (MQ=255)
                                                                             agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                                <  5:201654/135‑1 (MQ=255)
                                                                             agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTGCGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcac                                                                                >  6:201654/1‑135 (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                                                  >  7:293999/1‑149 (MQ=255)
                                                                                              gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                                  >  1:189164/1‑149 (MQ=255)
                                                                                              gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                                  >  4:248364/1‑149 (MQ=255)
                                                                                              gCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                                  <  3:248364/149‑1 (MQ=255)
                                                                                                    gTTCTTAGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATAACCGCCAGCACCACGAAGCCAACGGCCAGCAAGATACCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAg                                            >  4:193800/1‑149 (MQ=255)
                                                                                                        ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                                        >  6:86786/1‑149 (MQ=255)
                                                                                                        ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                                        >  4:133191/1‑149 (MQ=255)
                                                                                                        ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                                        >  4:77466/1‑149 (MQ=255)
                                                                                                        ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCa                                        >  7:216794/1‑149 (MQ=255)
                                                                                                            ggCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACg                                      >  3:105765/1‑148 (MQ=255)
                                                                                                             gCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                                    <  3:133191/149‑1 (MQ=255)
                                                                                                             gCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                                    <  7:88043/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                          <  2:340331/148‑1 (MQ=255)
                                                                                                                        aTCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg                          <  1:389704/148‑1 (MQ=255)
                                                                                                                         tCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg                        <  5:86786/149‑1 (MQ=255)
                                                                                                                         ccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCg                        >  8:102598/2‑149 (MQ=255)
                                                                                                                          cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                                     >  3:77493/1‑135 (MQ=255)
                                                                                                                          cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGa                                     <  4:77493/135‑1 (MQ=255)
                                                                                                                           ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt                      >  6:207971/1‑149 (MQ=255)
                                                                                                                           ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGAGCACGAAGTCGCCATAGCgg                       <  6:297479/148‑1 (MQ=255)
                                                                                                                               cccTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc                  <  8:164499/149‑1 (MQ=255)
                                                                                                                                      aCGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATg              >  5:332565/1‑147 (MQ=255)
                                                                                                                                      aCGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGaa            <  8:216794/149‑1 (MQ=255)
                                                                                                                                      aCGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGaa            <  6:266010/149‑1 (MQ=255)
                                                                                                                                        ggCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGt          >  7:362510/1‑149 (MQ=255)
                                                                                                                                        ggCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGt          <  8:323106/149‑1 (MQ=255)
                                                                                                                                            gTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCgat      <  8:362510/149‑1 (MQ=255)
                                                                                                                                               gCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATGGCGGTCGACCATGAAGTCgatgat   >  7:353813/1‑149 (MQ=255)
                                                                                                                                                cACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATc  <  6:423993/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATC  >  NC_002947/4911652‑4911941

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: