Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,848,920 A→T 6.0% K74* (AAG→TAG)  hemN → oxygen‑independent coproporphyrinogen III dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,848,9200AT 6.0% 204.5 / 2.3 67K74* (AAG→TAG) hemNoxygen‑independent coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

ACCCAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGCACACCGACCTTCCTCAGCCATG  >  NC_002947/4848774‑4849067
                                                                                                                                                  |                                                                                                                                                    
aCCCAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAg                                                                                                                                                    >  1:362483/1‑149 (MQ=255)
aCCCAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAg                                                                                                                                                    >  5:224947/1‑149 (MQ=255)
  ccAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGt                                                                                                                                                  >  6:145857/1‑149 (MQ=255)
   cAACCGCCGTGCCACTGCACAGCGAGGTGGGTTCGTTAGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTGCGACTTCGGCTGCAACAAGGt                                                                                                                                                  >  8:287320/1‑148 (MQ=255)
   cAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGtc                                                                                                                                                 <  8:137637/149‑1 (MQ=255)
      tcgctgTGCAGTTGCACGGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGtcatc                                                                                                                                              >  6:157321/2‑149 (MQ=255)
         ctgTGCAGTTGCACGGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCAc                                                                                                                                            >  7:288705/3‑148 (MQ=255)
         ccgTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCAc                                                                                                                                            >  5:70266/1‑148 (MQ=255)
               aaCTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGAc                                                                                                                                    >  5:215384/1‑149 (MQ=255)
                      caGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGAccgcgccc                                                                                                                             >  7:85334/1‑149 (MQ=255)
                               gggTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTgcgc                                                                                                                    >  2:126157/1‑149 (MQ=255)
                               gggTTCGTTCGACCTGCTAAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTgcgc                                                                                                                    >  2:122193/1‑149 (MQ=255)
                                ggTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTgcgc                                                                                                                    <  5:157321/148‑1 (MQ=255)
                                      ttcgACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCt                                                                                                             >  4:296713/1‑149 (MQ=255)
                                         gACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTg                                                                                                            >  6:253039/1‑147 (MQ=255)
                                          aCCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAg                                                                                                         >  6:126767/1‑149 (MQ=255)
                                              gCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGcc                                                                                                     >  7:275560/1‑149 (MQ=255)
                                               cTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCt                                                                                                    >  1:51764/1‑149 (MQ=255)
                                                  aaCGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGa                                                                                                 <  5:61536/149‑1 (MQ=255)
                                                  aaCGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGa                                                                                                 <  5:126767/149‑1 (MQ=255)
                                                  aaCGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGAACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGAGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTg                                                                                                   >  1:75325/1‑147 (MQ=255)
                                                     gCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAg                                                                                                >  4:168208/1‑147 (MQ=255)
                                                      cccTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAACAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCa                                                                                              <  1:251438/149‑1 (MQ=255)
                                                       ccTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGCCCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAg                                                                                             >  1:150166/1‑148 (MQ=255)
                                                           cgcgAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGaa                                                                                          <  7:8296/147‑1 (MQ=255)
                                                           cgcgAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGaa                                                                                          <  2:51764/147‑1 (MQ=255)
                                                            gcgAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATc                                                                                       <  8:275560/149‑1 (MQ=255)
                                                            gcgAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATc                                                                                       <  1:126157/149‑1 (MQ=255)
                                                            gcgAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATc                                                                                       <  1:122193/149‑1 (MQ=255)
                                                                   ccGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTTCAGCGCCTGGAGCAGGAAATCCAGCTgg                                                                                <  6:49527/149‑1 (MQ=255)
                                                                   ccGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTgg                                                                                <  5:145857/149‑1 (MQ=255)
                                                                          gCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTg                                                                         >  5:134858/1‑149 (MQ=255)
                                                                                   ccGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAACGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGa                                                                <  5:253039/149‑1 (MQ=255)
                                                                                   ccGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGa                                                                >  3:163452/1‑149 (MQ=255)
                                                                                    cgctgtcgctgtACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGAc                                                               <  4:352659/149‑1 (MQ=255)
                                                                                           gctgtACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTaaa                                                         <  6:70266/148‑1 (MQ=255)
                                                                                               tACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACg                                                    >  7:332555/1‑149 (MQ=255)
                                                                                               gaCGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAaacaac                                                     <  4:252399/147‑1 (MQ=255)
                                                                                                 cGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCGCCTGGACCCTAAACAACgt                                                   <  2:362483/148‑1 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACTAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAaacaa                                                      >  2:224308/1‑144 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACTAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAaacaa                                                      >  2:218901/1‑144 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACTAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAaacaa                                                      <  1:224308/144‑1 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACTAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAaacaa                                                      <  1:218901/144‑1 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACgtgt                                                 <  8:332555/149‑1 (MQ=255)
                                                                                                  gTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACgtgt                                                 <  3:168208/149‑1 (MQ=255)
                                                                                                    gCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGAcc                                                              >  7:79853/1‑134 (MQ=255)
                                                                                                    gCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGAcc                                                              <  8:79853/134‑1 (MQ=255)
                                                                                                      aTGTGGCGTTCTGCGCCATCATCTGCTACTACTGCGCCTGCTACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCGTACCTGCTGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGaa                                             <  6:224947/149‑1 (MQ=255)
                                                                                                           ccGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATcc                                                                                      >  6:195094/1‑103 (MQ=255)
                                                                                                           ccGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATcc                                                                                      <  5:195094/103‑1 (MQ=255)
                                                                                                            cGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGtt                                        >  1:271890/1‑148 (MQ=255)
                                                                                                                cTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTTGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCAt                                    <  7:287320/148‑1 (MQ=255)
                                                                                                                 tGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAg                                                                                    <  6:217600/99‑1 (MQ=255)
                                                                                                                 tGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAg                                                                                    >  5:217600/1‑99 (MQ=255)
                                                                                                                  gcgcCAACATCTGCTACTACTGCGTCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCAtttt                                 >  7:52677/1‑149 (MQ=255)
                                                                                                                   cgcCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCAtt                                   <  2:271890/146‑1 (MQ=255)
                                                                                                                      cAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACgtg                                                  >  7:128277/1‑128 (MQ=255)
                                                                                                                      cAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACgtg                                                  <  8:128277/128‑1 (MQ=255)
                                                                                                                      cAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTggcg                             >  5:315311/1‑149 (MQ=255)
                                                                                                                      cAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTggcg                             >  3:110145/1‑149 (MQ=255)
                                                                                                                            cTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCAt                                    >  7:339548/1‑136 (MQ=255)
                                                                                                                            cTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCAt                                    <  8:339548/136‑1 (MQ=255)
                                                                                                                              gCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGcaca                     >  7:281093/1‑149 (MQ=255)
                                                                                                                               ctactaCTGCGCCTGCAACACGGTCATCACCAAAGACCGCGCCCGAGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGGCTGCCACCTTGACGCTAAACAAGGTGTTGAACAGTTGCATTTTGGCGGCGGCacac                    <  2:75325/149‑1 (MQ=255)
                                                                                                                                     ctGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTggc                              >  2:343133/1‑133 (MQ=255)
                                                                                                                                     ctGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTggc                              <  1:343133/133‑1 (MQ=255)
                                                                                                                                          ccTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGCACACCGACCTTCCTc         <  6:315311/149‑1 (MQ=255)
                                                                                                                                                acaaGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGCACACCGACCTTCCTCAGCCAt   >  4:75522/1‑149 (MQ=255)
                                                                                                                                                 caaGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGACCAGTTGGATTTTGGCGGCGGCACACCGCCCTTCCTCAGCCTtg  >  4:222384/1‑149 (MQ=255)
                                                                                                                                                 caaGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGCACACCGACCTTCCTCAGCCAtg  >  8:45175/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCCAACCGCCGTGCAACTGCACAGCGAGGTGGGTTCGTTCGACCTGCTCAACGCCCTGCGCGAGAGCCGACGGGCCGTGCGCCCGCTGTCGCTGTACGTGCATGTGCCGTTCTGCGCCAACATCTGCTACTACTGCGCCTGCAACAAGGTCATCACCAAAGACCGCGCCCGCGCTGCGCCCTACCTGCAGCGCCTGGAGCAGGAAATCCAGCTGGTCGCCTGCCACCTGGACCCTAAACAACGTGTTGAACAGTTGCATTTTGGCGGCGGCACACCGACCTTCCTCAGCCATG  >  NC_002947/4848774‑4849067

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: