Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 857,577 C→T 16.9% G416D (GGT→GAT)  phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947857,5770CT16.9% 160.0 / 22.8 72G416D (GGT→GAT) phrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base C (34/25);  new base T (6/6);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

CCGCTGGCGACTGCTGCTGAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCC  >  NC_002947/857438‑857721
                                                                                                                                           |                                                                                                                                                 
ccgaaggCCACTGCTGCAGAGATCTACGATTGCACTTAGGTACGAAGTAATACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATAGCTTTTTAAACCAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACAGAACCAGAGGa                                                                                                                                           >  2:333573/6‑148 (MQ=255)
      gCGACTGCTGCTGAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGaa                                                                                                                                      <  6:254950/147‑1 (MQ=255)
             ctgctgAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGa                                                                                                                             <  6:146788/149‑1 (MQ=255)
                 tgAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACtgt                                                                                                                         >  6:79144/1‑149 (MQ=255)
                 tgAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATAGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACtgt                                                                                                                         >  5:48108/1‑149 (MQ=255)
                   agaTCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACtgtg                                                                                                                        >  8:268737/1‑148 (MQ=255)
                   agaTCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACtgt                                                                                                                         <  8:14996/147‑1 (MQ=255)
                   agaTCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCATGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACtgtg                                                                                                                        >  3:344010/1‑148 (MQ=255)
                        tACGAATGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATCGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCACGTAACCAGTGTCGAATGAACCGACCCTGTGGa                                                                                                                                           >  5:273585/1‑124 (MQ=255)
                         aCGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCg                                                                                                                  <  6:48108/148‑1 (MQ=255)
                               ggACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGa                                                                                                            <  3:44744/148‑1 (MQ=255)
                                     aGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCggg                                                                                                                >  4:217107/1‑138 (MQ=255)
                                     aGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCggg                                                                                                                <  3:217107/138‑1 (MQ=255)
                                         aCGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGaaa                                                                                                 >  6:7382/1‑149 (MQ=255)
                                           gAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATa                                                                                               <  7:268737/149‑1 (MQ=255)
                                                    ccGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACcgc                                                                                     >  1:303743/1‑149 (MQ=255)
                                                      gAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGt                                                                                                               >  1:80382/1‑122 (MQ=255)
                                                      gAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGt                                                                                                               <  2:80382/122‑1 (MQ=255)
                                                       aaaaGGTCCGTCGATTTAATGGGTTGATGGATGGATTTTTCATCGAGCTCCTGTACTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGCATCGAAGCGCTGCGACTTTGAAACCGCGTTGAAGATa                                                                                                       >  3:135782/1‑129 (MQ=255)
                                                            gTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGc                                                                             >  6:209261/1‑149 (MQ=255)
                                                            gTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGc                                                                             >  6:27050/1‑149 (MQ=255)
                                                            gTCCGTCGATTTAATGGGTAGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGc                                                                             <  4:358885/149‑1 (MQ=255)
                                                               cGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGAGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCg                                                                           <  2:303743/148‑1 (MQ=255)
                                                                    atatAATGGGTTGATGGATGGATTTTTCATCGAGCTCCTGTAGTTACGGTACCCAGAGGCGAATGAACCGACCAGGTGGATGGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTAGGTGCCGGTGGag                                                                     >  8:228149/4‑148 (MQ=255)
                                                                     tttAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAg                                                                    <  5:79144/149‑1 (MQ=255)
                                                                       tAATGGGTGGATGGATGGCTTTTTCAGCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGc                                                                   <  4:344010/148‑1 (MQ=255)
                                                                          tGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTc                                                               <  5:9137/149‑1 (MQ=255)
                                                                          tGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTc                                                               >  7:18857/1‑149 (MQ=255)
                                                                           gggTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTc                                                               >  5:260269/1‑148 (MQ=255)
                                                                           gggTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTc                                                               >  6:94261/1‑148 (MQ=255)
                                                                             gTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCAc                                                            <  5:209261/149‑1 (MQ=255)
                                                                              ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGGGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                           >  5:335674/1‑149 (MQ=255)
                                                                              ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTc                                                                                  >  5:325178/1‑126 (MQ=255)
                                                                              ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTc                                                                                  <  6:325178/126‑1 (MQ=255)
                                                                              ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTACGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCTCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                           >  4:89297/1‑149 (MQ=255)
                                                                                   ggatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCaa                                                      >  6:335001/1‑149 (MQ=255)
                                                                                   ggatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCaa                                                      >  2:84279/1‑149 (MQ=255)
                                                                                    gatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAAc                                                     >  7:49442/1‑149 (MQ=255)
                                                                                    gatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGAGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTTCCAAc                                                     <  8:18857/149‑1 (MQ=255)
                                                                                     atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTg                                                                                                             <  3:184805/93‑1 (MQ=255)
                                                                                     atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTg                                                                                                             >  4:184805/1‑93 (MQ=255)
                                                                                      tggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAAccc                                                   >  4:3436/1‑149 (MQ=255)
                                                                                      tggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGt                                                                               >  1:97283/1‑121 (MQ=255)
                                                                                      tggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGt                                                                               <  2:97283/121‑1 (MQ=255)
                                                                                        gCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAAccccc                                                 >  5:276622/1‑149 (MQ=255)
                                                                                               cATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGTACCGATCCTGTGGATCGAAGCTCTGCGACTGTGAAAGGGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCgttgtt                                           <  1:333573/148‑1 (MQ=255)
                                                                                               cATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCgttgt                                            <  1:84279/148‑1 (MQ=255)
                                                                                                aTCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGttgttg                                          <  8:49442/148‑1 (MQ=255)
                                                                                                 tCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggc                                        <  6:276622/149‑1 (MQ=255)
                                                                                                        ccTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       <  1:228180/143‑1 (MQ=255)
                                                                                                        ccTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       >  2:228180/1‑143 (MQ=255)
                                                                                                           gTAGTTCCGGTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCg                              >  6:355925/1‑149 (MQ=255)
                                                                                                           gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       >  1:167673/1‑140 (MQ=255)
                                                                                                           gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       <  2:167673/140‑1 (MQ=255)
                                                                                                           gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       >  3:194614/1‑140 (MQ=255)
                                                                                                           gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcg                                       <  4:194614/140‑1 (MQ=255)
                                                                                                               ttCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGt                          >  5:95904/1‑149 (MQ=255)
                                                                                                                  cGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCt                        <  6:95904/148‑1 (MQ=255)
                                                                                                                  cGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTa                       >  1:47053/1‑149 (MQ=255)
                                                                                                                    gTACCCAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATc                     <  5:355925/149‑1 (MQ=255)
                                                                                                                       cccAGTGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAgg                  <  5:7382/149‑1 (MQ=255)
                                                                                                                         cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGTGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTg                <  3:89297/149‑1 (MQ=255)
                                                                                                                           gTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGCGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGATGTTGGCGGACAGGTCGCCGTCTATCAGGTGgc              >  5:338655/1‑149 (MQ=255)
                                                                                                                            tGGCGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGgcg             >  7:306142/1‑149 (MQ=255)
                                                                                                                               cGAATGAACCGATCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt          <  3:3436/149‑1 (MQ=255)
                                                                                                                               cGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt          >  5:164129/1‑149 (MQ=255)
                                                                                                                                gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt          >  3:352439/1‑148 (MQ=255)
                                                                                                                                gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt          >  8:137683/1‑148 (MQ=255)
                                                                                                                                  aTGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGac       <  5:94261/149‑2 (MQ=255)
                                                                                                                                    gAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGt     >  6:308230/1‑149 (MQ=255)
                                                                                                                                      aCCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGc   >  5:16427/1‑149 (MQ=255)
                                                                                                                                       ccGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCACCCCCCGTGGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGcc  >  4:201974/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CCGCTGGCGACTGCTGCTGAGATCTACGATTGGACTTAGGTACGAAGTATTACCGAAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCC  >  NC_002947/857438‑857721

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: