Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 857,884 G→T 8.7% P314T (CCG→ACG)  phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947857,8840GT 8.7% 202.2 / 8.4 69P314T (CCG→ACG) phrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base G (25/38);  new base T (2/4);  total (27/42)
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.79e-01

CAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAAC  >  NC_002947/857736‑858032
                                                                                                                                                    |                                                                                                                                                     
cAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                                                       >  8:194483/1‑149 (MQ=255)
cAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCACCAATGATCGGGACGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                                                       >  8:130652/1‑149 (MQ=255)
cAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGAGCATCCACCCGCTGAGCAGAAATTGGCGCATGGCGGCACCTATCATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgg                                                                                                                                                       >  8:134334/1‑149 (MQ=255)
   tCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGc                                                                                                                                                    <  6:119434/149‑1 (MQ=255)
        cagcagGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgggg                                                                                                                                                     <  5:186394/143‑1 (MQ=255)
        cagcagGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCgggg                                                                                                                                                     >  6:186394/1‑143 (MQ=255)
            agGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCca                                                                                                                                           <  5:144063/149‑1 (MQ=255)
             ggTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCcag                                                                                                                                          <  5:140475/149‑1 (MQ=255)
             ggTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCcag                                                                                                                                          <  8:150523/149‑1 (MQ=255)
              gTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCcag                                                                                                                                          <  4:99864/148‑1 (MQ=255)
                 cTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCcagggc                                                                                                                                       <  6:257472/148‑1 (MQ=255)
                    gCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCcagggcaggg                                                                                                                                   >  1:55297/1‑149 (MQ=255)
                    gCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGGGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCCTGGCGTGACGCcccggca                                                                                                                                      <  6:30713/146‑1 (MQ=255)
                     cTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCCTGGCGTGACGCcccggcagacc                                                                                                                                  <  2:158594/149‑4 (MQ=255)
                      tGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCt                                                                                                                                 >  1:266184/1‑149 (MQ=255)
                         ggAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCg                                                                                                                              <  3:89360/149‑1 (MQ=255)
                           aaCATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTCCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGt                                                                                                                            >  5:137669/1‑149 (MQ=255)
                            aCATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTc                                                                                                                           >  6:71880/1‑149 (MQ=255)
                               tGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTgg                                                                                                                        <  2:36163/149‑1 (MQ=255)
                                 gCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGGCACGCCAGGGCAGGGCTTCGGTCTggg                                                                                                                       <  6:266271/148‑1 (MQ=255)
                                        ccATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAg                                                                                                               >  7:179620/1‑149 (MQ=255)
                                           tGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCAc                                                                                                                                               >  5:144824/1‑114 (MQ=255)
                                           tGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCAc                                                                                                                                               <  6:144824/114‑1 (MQ=255)
                                              gcAAGCGGTTGTGCATCCACCCGTTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCgg                                                                                                         >  7:169213/1‑149 (MQ=255)
                                              gcAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGACTTCGGTCTGGGCGCGGAAGGCGCgg                                                                                                         >  8:247441/1‑149 (MQ=255)
                                                aaGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGCAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTg                                                                                                       <  5:216356/149‑1 (MQ=255)
                                                aaGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTg                                                                                                       <  7:22564/149‑1 (MQ=255)
                                                  gCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGgc                                                                                                     >  4:155479/1‑149 (MQ=255)
                                                  gCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCACCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGgc                                                                                                     >  2:257128/1‑149 (MQ=255)
                                                      ttGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGcgcg                                                                                                  >  1:230044/1‑148 (MQ=255)
                                                        gtgCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCgaga                                                                                               >  6:86504/1‑149 (MQ=255)
                                                             tCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                          <  2:266184/149‑1 (MQ=255)
                                                             tCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgt                                                                                          <  8:194242/149‑1 (MQ=255)
                                                               cACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCtgtg                                                                                         <  8:169213/148‑1 (MQ=255)
                                                                    ggTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAAc                                                                                   <  3:155479/149‑1 (MQ=255)
                                                                         gcagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCt                                                                               >  6:124726/1‑148 (MQ=255)
                                                                          cagcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCa                                                                             <  2:55297/149‑1 (MQ=255)
                                                                           agcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCACGGTGGCGCGAGACCTGTGGATAACCGCTCAg                                                                            >  8:85466/1‑149 (MQ=255)
                                                                            gcaATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGACGCGGTGGCGCGAGACCTGTGGATAACCGCTCAg                                                                            >  4:66240/1‑148 (MQ=255)
                                                                               aTTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCACGGTGGCGCGAGACCTGTGGATAACCGCTCAGGAtg                                                                        <  7:85466/149‑1 (MQ=255)
                                                                                   gcgcATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGc                                                                     <  5:71880/148‑1 (MQ=255)
                                                                                     gcATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTg                                                                  <  8:179620/149‑1 (MQ=255)
                                                                                      cATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGACTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGt                                                                 <  7:247441/149‑1 (MQ=255)
                                                                                          gcggcATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGaa                                                            >  7:231071/1‑149 (MQ=255)
                                                                                           cggcATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAAc                                                           <  7:194483/149‑1 (MQ=255)
                                                                                                ccAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACg                                                      <  1:257128/148‑1 (MQ=255)
                                                                                                  aaTGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGc                                                     <  2:230044/148‑1 (MQ=255)
                                                                                                     gATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCa                                                   >  6:30271/1‑147 (MQ=255)
                                                                                                       tCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAgg                                                                                                              <  7:123989/87‑1 (MQ=255)
                                                                                                       tCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGACAGGGCAGGGCTTCGGTCTGGGCGCGGAAgg                                                                                                              >  8:123989/1‑87 (MQ=255)
                                                                                                                     cGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATcc                                  <  5:124726/148‑1 (MQ=255)
                                                                                                                        cccTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGt                              >  7:126595/1‑149 (MQ=255)
                                                                                                                        cccTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGCCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGt                              <  3:66240/149‑1 (MQ=255)
                                                                                                                          cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAg                                <  2:148012/145‑1 (MQ=255)
                                                                                                                          cTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAg                                >  1:148012/1‑145 (MQ=255)
                                                                                                                            gCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCAAGCCAGGGCAGGGCTTCGGTCTGGCCGCGGAAGGCGCGGTGGCGCGAGACCTGAGGATAACCGCTCAGGATGTGCTTGAAGAACTCACGCCAGAGAAGCTCATTGATCCAAGTCTg                           >  2:209137/1‑148 (MQ=255)
                                                                                                                             cTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAATGCGCGGTGGCGCGAGACCTGTTGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGc                          <  6:137669/148‑1 (MQ=255)
                                                                                                                              tCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCAc                        >  2:118216/1‑149 (MQ=255)
                                                                                                                              tCCCAGGCCGCCAGGTCTGCCTGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCaga                                                 >  6:40478/1‑123 (MQ=255)
                                                                                                                              tCCCAGGCCGCCAGGTCTGCCTGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCaga                                                 <  5:40478/123‑1 (MQ=255)
                                                                                                                                     ccgccAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCtgtg                 <  8:231071/149‑1 (MQ=255)
                                                                                                                                      cgccAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGtgc                >  1:98835/1‑149 (MQ=255)
                                                                                                                                        ccAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGtgctg              >  4:20068/1‑149 (MQ=255)
                                                                                                                                           ggTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCa           <  7:134334/149‑1 (MQ=255)
                                                                                                                                           ggTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCa           <  7:130652/149‑1 (MQ=255)
                                                                                                                                           ggTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCa           <  5:108222/149‑1 (MQ=255)
                                                                                                                                            gTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGc                      <  5:238821/137‑1 (MQ=255)
                                                                                                                                            gTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGc                      >  6:238821/1‑137 (MQ=255)
                                                                                                                                             ccTGCCCTGGCGTGACGCCCCGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCa                         >  6:134156/2‑133 (MQ=255)
                                                                                                                                             ccTGCCCTGGCGTGACGCCCCGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCa                         <  5:134156/132‑1 (MQ=255)
                                                                                                                                                 ccGGGTCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTc      <  6:271364/148‑1 (MQ=255)
                                                                                                                                                 ccGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTg        <  5:30271/146‑1 (MQ=255)
                                                                                                                                                 ccGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGt       <  8:126595/147‑1 (MQ=255)
                                                                                                                                                 ccGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGGGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTATGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGt       <  1:209137/148‑1 (MQ=255)
                                                                                                                                                    gggCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAAc  >  8:152119/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAGCGGTTGTGCATCCACCCGGTGTGCAGCAATTGGCGCATGGCGGCATCAATGATCGGGAAGCCGGTGCGGCCCTGCTCCCAGGCCGCCAGGTCTGCCGGGGCGTCACGCCAGGGCAGGGCTTCGGTCTGGGCGCGGAAGGCGCGGTGGCGCGAGACCTGTGGATAACCGCTCAGGATGTGCTTGTAGAACTCACGCCAGAGTAGCTCATTGATCCAGGTCTGCACGCCTGTGCTGCCACTGTCGAAC  >  NC_002947/857736‑858032

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: