Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,704,190 C→A 5.8% S604S (TCG→TCT paaZ ← bifunctional oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,704,1900CA 5.8% 221.0 / 3.6 69S604S (TCG→TCTpaaZbifunctional oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base A (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCAAGTCAGCTTCT  >  NC_002947/3704042‑3704337
                                                                                                                                                    |                                                                                                                                                   
ccTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCGGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGc                                                                                                                                                     >  5:201019/1‑149 (MQ=255)
    cGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCGTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATATCGGTCCGGCGCCGGGCGa                                                                                                                                                   <  3:240888/147‑1 (MQ=255)
     ggCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGGTGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGacac                                                                                                                                                <  3:288422/149‑1 (MQ=255)
     ggCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCTCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGacac                                                                                                                                                <  6:255436/149‑1 (MQ=255)
           cccAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAg                                                                                                                                          >  3:312706/1‑149 (MQ=255)
            ccAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGc                                                                                                                                         >  4:13254/1‑149 (MQ=255)
             cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGcc                                                                                                                                        <  2:131877/149‑1 (MQ=255)
             cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGcc                                                                                                                                        >  6:313098/1‑149 (MQ=255)
                 gggCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCGGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGc                                                                                                                                    <  6:201019/149‑1 (MQ=255)
                     tGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCt                                                                                                                                >  8:197240/1‑149 (MQ=255)
                            gCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGa                                                                                                                           >  8:256683/1‑147 (MQ=255)
                            gCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGa                                                                                                                           <  3:13254/147‑1 (MQ=255)
                             cccTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACa                                                                                                                         <  5:313098/148‑1 (MQ=255)
                               cTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTa                                                                                                                      >  6:162946/1‑149 (MQ=255)
                                 ggTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTAc                                                                                                                     >  4:215282/1‑148 (MQ=255)
                                   tCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCTCCAATGCCAACCGGGTTGATGAAGCTCAGGGTGTCCAGCCCATAGTTGTCAAATACCGGTCCGGCGCCGGGAGCCACGAACAGCCCGGCCGCTGCCGACAGTACg                                                                                                                    >  1:276427/1‑147 (MQ=255)
                                   tCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGaa                                                                                                                  <  3:201759/149‑1 (MQ=255)
                                   tCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGAGACACGAACAGCCCGGCCGCTGCCGACAGTACg                                                                                                                    >  1:272488/1‑147 (MQ=255)
                                   tCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGAGACACGAACAGCCCGGCCGCTGCCGACAGTACg                                                                                                                    <  2:272488/147‑1 (MQ=255)
                                   tCGATCTTGAGCTTGCAGGTCAGGCTCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGAGACACGAACAGCCCGGCCGCTGCCGACAGTACg                                                                                                                    <  2:276427/147‑1 (MQ=255)
                                         ttGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAAcc                                                                                                            >  4:328217/1‑149 (MQ=255)
                                         ttGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAAcc                                                                                                            <  2:328422/149‑1 (MQ=255)
                                            cgcTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGt                                                                                                          <  7:197240/148‑1 (MQ=255)
                                              cTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGAAGGGTGTCCAGCCCATAGTTGGCCAAAACCGGTCCGGCGCCGGGCGACACGAAAAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTggg                                                                                                       >  8:72064/1‑149 (MQ=255)
                                                  cAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAAt                                                                                                   <  5:162946/149‑1 (MQ=255)
                                                     gTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTaa                                                                                                              >  8:190528/1‑135 (MQ=255)
                                                     gTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCg                                                                                                >  4:297345/1‑149 (MQ=255)
                                                     gTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCg                                                                                                <  6:268953/149‑1 (MQ=255)
                                                     gTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCCGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTaa                                                                                                              <  7:190528/135‑1 (MQ=255)
                                                         ggCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTgg                                                                                                        >  7:27018/1‑137 (MQ=255)
                                                         ggCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTgg                                                                                                        <  8:27018/137‑1 (MQ=255)
                                                                 tGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTg                                                                                    >  2:234017/1‑149 (MQ=255)
                                                                   gATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACCCGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAAGGCGTTTACCGAACGGc                                                                                   >  8:189016/1‑148 (MQ=255)
                                                                   gATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTg                                                                                    <  3:215282/147‑1 (MQ=255)
                                                                   gATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGAACGAAGTAACCGTGGGCAATGCGTTTACCGAACtcc                                                                                   >  3:294329/1‑146 (MQ=255)
                                                                       gTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGc                                                                              >  7:60606/1‑149 (MQ=255)
                                                                                    ccACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGAtt                                                                 <  7:256683/149‑1 (MQ=255)
                                                                                    ccACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGAtt                                                                 <  2:132051/149‑1 (MQ=255)
                                                                                      aCCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTtc                                                               >  6:330751/1‑149 (MQ=255)
                                                                                          ggTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTtcgtc                                                            >  2:347064/1‑148 (MQ=255)
                                                                                              gatgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGt                                                                                                                       >  6:318312/1‑85 (MQ=255)
                                                                                              gatgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGt                                                                                                                       <  5:318312/85‑1 (MQ=255)
                                                                                              gatgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTtcgtc                                                            >  4:300820/1‑144 (MQ=255)
                                                                                              gatgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTtcgtc                                                            <  3:300820/144‑1 (MQ=255)
                                                                                              gatgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAg                                                        <  3:328217/148‑1 (MQ=255)
                                                                                                tgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGc                                                     <  5:330751/149‑1 (MQ=255)
                                                                                                tgaAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGc                                                     <  3:297345/149‑1 (MQ=255)
                                                                                                        aGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGaa                                             <  1:234017/149‑1 (MQ=255)
                                                                                                          ggTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGc                                                     >  4:169148/1‑139 (MQ=255)
                                                                                                          ggTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGc                                                     <  3:169148/139‑1 (MQ=255)
                                                                                                                ctgCCCATAGTTGGCCAATACCGGTGCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGc                                     <  4:312706/147‑1 (MQ=255)
                                                                                                                 aGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGcc                                    >  6:298442/1‑149 (MQ=255)
                                                                                                                    ccATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGc                                     >  7:215759/1‑145 (MQ=255)
                                                                                                                    ccATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGc                                     <  8:215759/145‑1 (MQ=255)
                                                                                                                      aTAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGa                                                                   <  1:43871/113‑1 (MQ=255)
                                                                                                                      aTAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGa                                                                   >  2:43871/1‑113 (MQ=255)
                                                                                                                      aTAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGa                                 >  3:37326/1‑147 (MQ=255)
                                                                                                                         gTTGGACAATACCGGTCCGGCGCCGGGCGACACGATGAGCCCGGCCGGTGCCGACAGTTCGAAGTAACGGTTGGCAATGCGTTTACAGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAgg                             <  7:72064/148‑1 (MQ=255)
                                                                                                                           tGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAg                          >  1:288967/1‑149 (MQ=255)
                                                                                                                             gCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGc                         <  2:288967/148‑1 (MQ=255)
                                                                                                                                 aTACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGaaa                    <  8:60606/149‑1 (MQ=255)
                                                                                                                                      ggTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGGGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGAc               >  1:70156/1‑149 (MQ=255)
                                                                                                                                          cGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGAc               >  5:229482/1‑145 (MQ=255)
                                                                                                                                          cGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGAc               <  6:229482/145‑1 (MQ=255)
                                                                                                                                          cGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCaa            <  4:294329/148‑1 (MQ=255)
                                                                                                                                           ggCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCTTGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCAAGt          >  8:8017/1‑149 (MQ=255)
                                                                                                                                              gcCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGAc                                >  3:244348/1‑124 (MQ=255)
                                                                                                                                              gcCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGAc                                <  4:244348/124‑1 (MQ=255)
                                                                                                                                              gcCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCAAGTCAg       <  1:347064/149‑1 (MQ=255)
                                                                                                                                                   gCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCAAGTCAGCTTCt  >  8:253124/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTGGCGGCGATTTCGTCGAAGTGCATGTAGAAATGGTCGCCTGACAGGCAGCCGAAATTGACCAAGTCAGCTTCT  >  NC_002947/3704042‑3704337

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: