Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,094,261 A→G 7.7% T420A (ACC→GCC)  PP_3602 → 2,5‑dioxovalerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,094,2610AG 7.7% 182.9 / 4.8 66T420A (ACC→GCC) PP_36022,5‑dioxovalerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (3/2);  total (34/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.19e-01

CTGTTCAAGGCCGATGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATGGTGCACGGCGGCCCG  >  NC_002947/4094120‑4094404
                                                                                                                                             |                                                                                                                                                
cTGTTCAAGGCCGATGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGAc                                                                                                                                           >  7:307522/1‑149 (MQ=255)
   ttCAAGGCCGATGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCt                                                                                                                                        >  2:179421/1‑149 (MQ=255)
       aGGCCGATGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATc                                                                                                                                    >  8:332762/1‑149 (MQ=255)
             aTGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCGGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCg                                                                                                                                >  5:313640/1‑147 (MQ=255)
             aTGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGaa                                                                                                                              <  1:179421/149‑1 (MQ=255)
               gCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAg                                                                                                                             <  4:303790/148‑1 (MQ=255)
               gCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAg                                                                                                                             >  3:303790/1‑148 (MQ=255)
                   gCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCggag                                                                                                                        <  3:293050/149‑1 (MQ=255)
                      tgctgGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGAc                                                                                                                     >  7:192430/1‑149 (MQ=255)
                         tgGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTc                                                                                                                  >  2:156140/1‑149 (MQ=255)
                            tCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGAt                                                                                                               >  7:8977/1‑149 (MQ=255)
                             cGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATg                                                                                                              >  7:54632/1‑149 (MQ=255)
                              gAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGc                                                                                                             <  4:340865/149‑1 (MQ=255)
                                    gACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGCAGCGGAGGACTTCg                                                                                                                 >  8:8709/1‑139 (MQ=255)
                                    gACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGc                                                                                                       <  3:261606/149‑1 (MQ=255)
                                       ccGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCctg                                                                                                    >  5:322374/1‑149 (MQ=255)
                                           tGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCctggctg                                                                                               <  8:307522/149‑1 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGTCAGCTGGCcgcgcag                                                                                                                                          <  8:289325/103‑4 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGTCAGCTGGCcgcgcag                                                                                                                                          >  7:289325/1‑100 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGcc                                                                                                            >  2:348914/1‑133 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGcc                                                                                                            <  1:348914/133‑1 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGc                                                                                           >  5:100726/1‑149 (MQ=255)
                                               gCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTACGCCTGGCTGGTGc                                                                                           >  8:217514/1‑149 (MQ=255)
                                                cAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGc                                                                                           <  7:332762/148‑1 (MQ=255)
                                                     agaagTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCgctgc                                                                                     <  7:68860/149‑1 (MQ=255)
                                                            tttGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGc                                                                                                             <  8:71862/119‑1 (MQ=255)
                                                            tttGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGc                                                                                                             >  7:71862/1‑119 (MQ=255)
                                                               gggCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGgaga                                                                            <  8:8977/148‑1 (MQ=255)
                                                                 gCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGc                                                                                                       >  8:237973/1‑120 (MQ=255)
                                                                 gCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGc                                                                                                       <  7:237973/120‑1 (MQ=255)
                                                                 gCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAg                                                                         >  3:341970/1‑149 (MQ=255)
                                                                   ccaccaCCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTGCATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTc                                                                       >  7:294675/1‑149 (MQ=255)
                                                                                 gTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGTCgctgct                                                                                    >  3:183461/1‑122 (MQ=255)
                                                                                 gTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGTCgctgct                                                                                    <  4:183461/122‑1 (MQ=255)
                                                                                     agGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAAt                                                     <  4:341970/149‑1 (MQ=255)
                                                                                     agGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAAt                                                     <  8:54632/149‑1 (MQ=255)
                                                                                       gTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTAGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTg                                                           <  1:237433/142‑1 (MQ=255)
                                                                                       gTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTg                                                           >  2:237433/1‑141 (MQ=255)
                                                                                       gTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATgg                                                   >  6:95728/1‑149 (MQ=255)
                                                                                           aGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTAc                                               <  1:156140/149‑1 (MQ=255)
                                                                                             gACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATgg                                                   >  1:254587/1‑143 (MQ=255)
                                                                                             gACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATgg                                                   <  2:254587/143‑1 (MQ=255)
                                                                                             gACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTAccc                                             >  8:268001/1‑149 (MQ=255)
                                                                                               cAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCt                                                            <  1:74458/132‑1 (MQ=255)
                                                                                               cAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCt                                                            >  2:74458/1‑132 (MQ=255)
                                                                                                aaCGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGAc                                          <  8:192430/149‑1 (MQ=255)
                                                                                                  cGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCg                                        <  7:268001/149‑1 (MQ=255)
                                                                                                            cgcgACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTTGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGaa                                <  6:100726/147‑1 (MQ=255)
                                                                                                            cgcgACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAg                               <  6:313640/148‑1 (MQ=255)
                                                                                                            cgcgACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGt                              <  5:95728/149‑1 (MQ=255)
                                                                                                               gACTCGCTGCGCTGCCTGCGGGGCAAGCTGACCTCGACGCTGATCGGCGAATCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGctgctg                                                                                   <  8:141888/93‑1 (MQ=255)
                                                                                                               gACGCGCTGCTCGGCCTGCCCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGCCTTCGATGCCTTCGCCTGGCTGGTGCCGctgctg                                                                                   >  7:141888/1‑93 (MQ=255)
                                                                                                               gACGCGCTCCTCGGCCTGCGCGGACAGCAGACAGCCACGCTGATCGGCGAAGCGCAGGACTTCGATGCCTTCGCCTGTCTTGTGCCGCTGCAGCAGGAGAAAGTCGCGCGAATCCTGGTCAATGGCTAACCGAACGGTGACGAAGTCtg                           >  6:178904/1‑149 (MQ=255)
                                                                                                                aCGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCtgt                          <  5:178904/149‑1 (MQ=255)
                                                                                                                         tgggCCTGCGCGGCCAGCTGTCCGCGTCGCTGTTCGGCGAATCGGAGGAATTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCAt                  <  7:217514/146‑1 (MQ=255)
                                                                                                                         tCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATg                 >  4:79284/1‑149 (MQ=255)
                                                                                                                             gcgtcgcgGTCAGCTGGCCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATGGTGc             >  3:94750/5‑149 (MQ=255)
                                                                                                                              cTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGt                              >  8:139104/1‑131 (MQ=255)
                                                                                                                              cTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGt                              <  7:139104/131‑1 (MQ=255)
                                                                                                                               gtcgcgGTCAGCTGGCCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATGGTGCAc           >  7:86046/3‑149 (MQ=255)
                                                                                                                                    ggCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAg                                                                         <  6:153951/82‑1 (MQ=255)
                                                                                                                                        aGCTGGCCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGGCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAATTCTGTGATGCCATGGTGCACGGCGGccc   <  4:94750/148‑1 (MQ=255)
                                                                                                                                         gCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATGGTGAACGGCGGCCCg  >  5:275493/1‑148 (MQ=255)
                                                                                                                                          cTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGt                              >  8:3589/1‑119 (MQ=255)
                                                                                                                                            gACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATgg                <  3:106129/131‑1 (MQ=255)
                                                                                                                                            gACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATgg                >  4:106129/1‑131 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CTGTTCAAGGCCGATGCACGCCTGCTGGTCGAGTCCGACCCGCTGTTGCAGGAAGAAGTGTTTGGGCCCACCACCGTGGCCGTAGAGGTGCAGGACAACGACCAACTGCGCGACGCGCTGCTCGGCCTGCGCGGCCAGCTGACCGCGACGCTGATCGGCGAAGCGGAGGACTTCGATGCCTTCGCCTGGCTGGTGCCGCTGCTGGAGGAGAAAGTCGGGCGAATCCTGGTCAATGGCTACCCGACCGGTGTCGAAGTCTGTGATGCCATGGTGCACGGCGGCCCG  >  NC_002947/4094120‑4094404

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: