Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,455,567 T→G 7.6% E108D (GAA→GAC aacs ← acetoacetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,455,5670TG 7.6% 202.1 / 7.6 66E108D (GAA→GACaacsacetoacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base G (2/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

GGCCACCAGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTGGGGTATGCCAA  >  NC_002947/3455426‑3455708
                                                                                                                                              |                                                                                                                                                
ggCCACCAGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCACGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTa                                                                                                                                         >  8:320690/1‑149 (MQ=255)
  ccaccaGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACg                                                                                                                                       >  1:179704/1‑149 (MQ=255)
   caccaGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACgg                                                                                                                                      >  1:124551/1‑149 (MQ=255)
     ccaGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGc                                                                                                                                     >  3:235415/1‑148 (MQ=255)
      caGTGTTAGCCAGGTATTGGGCATGATGGCGGCTACCCGGTAGCCCGGGAAGATACCTACAGCGTTGAATGCAGTTTTCAGGCCTGCCACCTGAGCTCCCAGTTCGGCATGCGTGAACCCTTCGCGCTGGCCGTCTTCGACTACGGCGa                                                                                                                                   >  8:167901/1‑149 (MQ=255)
               ccAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACgg                                                                                                                             <  3:72300/146‑1 (MQ=255)
               ccAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACgg                                                                                                                             >  4:72300/1‑146 (MQ=255)
               ccAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGcc                                                                                                                           >  1:268555/1‑148 (MQ=255)
               ccAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGcc                                                                                                                           <  2:268555/148‑1 (MQ=255)
               ccAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCg                                                                                                                          >  3:172426/1‑149 (MQ=255)
                      ttGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGgt                                                                                                                    <  2:179704/148‑1 (MQ=255)
                                 ggcggcTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCgcggcggc                                                                                                       >  5:48123/1‑149 (MQ=255)
                                  gcggcTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCgcggcggc                                                                                                       <  6:48123/148‑1 (MQ=255)
                                   cggcTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGgcgc                                                                                                     <  3:241316/149‑1 (MQ=255)
                                      cTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCgcagc                                                                                                  <  7:320690/149‑1 (MQ=255)
                                        tcccGGGCGACCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCGCCGTGTGGGGCCGGTTCGGCTTGGGTGACCAATTCGCGGTGGCCGAATTCGCGTACGGCGCCCAGGGCCGGGCGGTCGTGgcggcggcgagg                                                                                                   >  1:91876/2‑144 (MQ=255)
                                         cccGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCgcagca                                                                                                 <  4:172426/147‑1 (MQ=255)
                                          ccGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGgt                                                                                              >  6:18971/1‑149 (MQ=255)
                                          ccGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGgt                                                                                              >  6:24147/1‑149 (MQ=255)
                                          ccGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGgt                                                                                              >  4:60762/1‑149 (MQ=255)
                                            ggTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGgtgt                                                                                            >  1:284745/1‑149 (MQ=255)
                                                gCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCgg                                                                                        <  2:124551/149‑1 (MQ=255)
                                                   cGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTGGTCGCCGTACGCGAAGACGGccag                                                                                                                         <  8:157227/114‑3 (MQ=37)
                                                   cGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTGGTCGCCGTACGCGAAGACGGccag                                                                                                                         >  7:157227/1‑112 (MQ=37)
                                                           tacctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGTCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCg                                                                                         >  6:248846/1‑137 (MQ=255)
                                                           tacctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGTCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCg                                                                                         <  5:248846/137‑1 (MQ=255)
                                                           tacctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCggg                                                                                                                        <  6:246163/107‑1 (MQ=255)
                                                           tacctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCggg                                                                                                                        >  5:246163/1‑107 (MQ=255)
                                                           tacctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAgg                                                                             <  4:149708/149‑1 (MQ=255)
                                                             cctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGt                                                                           <  4:235415/149‑1 (MQ=255)
                                                             cccgatgcacagtggttcgccaatgccaccctgaactttgccgaacacctgctgcgccgccgcgacgaccgcccGGCCGGCTTCGCGTACGGCGACCACGGCCGGGCg                                                                                                                      <  2:248563‑M1/34‑1 (MQ=255)
                                                              ctaccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCGGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTg                                                                          >  5:94523/1‑149 (MQ=255)
                                                                 ccGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCa                                                                       >  7:122729/1‑149 (MQ=255)
                                                                     cTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGcagcag                                                                                                <  8:45829/120‑1 (MQ=255)
                                                                     cTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGcagcag                                                                                                >  7:45829/1‑120 (MQ=255)
                                                                         aaTGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGa                                                                <  2:284745/148‑1 (MQ=255)
                                                                          aTGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAAc                                                              <  5:18971/149‑1 (MQ=255)
                                                                          aTGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAAc                                                              <  5:24147/149‑1 (MQ=255)
                                                                            gCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCa                                                            >  1:92479/1‑149 (MQ=255)
                                                                              cGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACt                                                          >  2:315640/1‑149 (MQ=255)
                                                                                    cAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGCCCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCAt                                                    >  3:27524/1‑149 (MQ=255)
                                                                                       gccggccACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCgg                                                 <  3:60762/149‑1 (MQ=255)
                                                                                          ggccACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCa                                              >  1:219857/1‑149 (MQ=255)
                                                                                              aCCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCGGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTg                                          <  6:94523/149‑1 (MQ=255)
                                                                                                      cGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGccc                                    <  2:219857/148‑1 (MQ=255)
                                                                                                       ggCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCtt                                  >  6:310725/1‑149 (MQ=255)
                                                                                                       ggCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCAATTGCgg                                       >  8:40263/1‑144 (MQ=255)
                                                                                                       ggCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCAATTGCgg                                       <  7:40263/144‑1 (MQ=255)
                                                                                                       ggCAAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCtt                                  >  4:157938/1‑149 (MQ=255)
                                                                                                         ccAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCa                                                                               <  4:251857/102‑1 (MQ=255)
                                                                                                         ccAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCa                                                                               >  3:251857/1‑102 (MQ=255)
                                                                                                         ccAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCg                                >  3:287068/1‑149 (MQ=255)
                                                                                                               cGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAg                           <  5:310725/148‑1 (MQ=255)
                                                                                                                        tGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCAttt                                           <  5:209673/123‑1 (MQ=255)
                                                                                                                        tGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCAttt                                           >  6:209673/1‑123 (MQ=255)
                                                                                                                        tGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACctggctg                 >  8:212819/1‑149 (MQ=255)
                                                                                                                              cTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCtgggtggg           <  2:92479/149‑1 (MQ=255)
                                                                                                                               ttCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTgg                                                                   >  1:64584/1‑92 (MQ=255)
                                                                                                                               ttCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTgg                                                                   <  2:64584/92‑1 (MQ=255)
                                                                                                                               ttCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTgggg          >  3:10060/1‑149 (MQ=255)
                                                                                                                                 cgcgCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTGGGGt         <  4:287068/148‑1 (MQ=255)
                                                                                                                                  cgcccGGCCGTGGTCGCCGTACGCGAAGACGGCCGGGCGGTCGTCGGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCggg                                                <  1:248563‑M1/104‑1 (MQ=255)
                                                                                                                                   gcccGGCCGTGGTCGCCGTACGCGAAGACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGa                            <  7:42835‑M1/123‑1 (MQ=255)
                                                                                                                                   gcccGGCCGTGGTCGCCGTACGCGAAGACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGa                            >  8:42835‑M1/5‑127 (MQ=255)
                                                                                                                                      tGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTGGGGTATGcc    >  7:4828/1‑148 (MQ=255)
                                                                                                                                        GCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTGGGGTATGCCaa  >  7:167901‑M1/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGCCACCAGTGTTTGCCAGGTATTGGGCATGATGGCGGCTACCCGGTCGCCCGGGACGATACCTACCGCCTTGAATGCCGTTTGCAGGCCGGCCACCTGTGCGGCCAGTTCGGCATGGGTGAACACTTCGCGCTGGCCGTCTTCGCGTACGGCGACCACGGCCGGGCGGTCGTCGCGGCGGCGCAGCAGGTGTTCGGCAAAGTTCAGGGTGGCATTGGCGAACCACTGTGCATCGGGCATTTGCGGCCCTTCGCTGAGCACCTGGCTGGGTGGGGTATGCCAA  >  NC_002947/3455426‑3455708

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: