Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,584,898 G→T 7.4% intergenic (+45/‑90) benC → / → benD benzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,584,8980GT 7.4% 258.5 / 7.7 81intergenic (+45/‑90)benC/benDbenzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base T (2/4);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCGG  >  NC_002947/3584752‑3585046
                                                                                                                                                  |                                                                                                                                                     
tgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa                                                                                                                                                      <  1:323525/148‑1 (MQ=255)
       cGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTc                                                                                                                                              <  6:295264/149‑1 (MQ=255)
        gCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                             <  6:328561/149‑1 (MQ=255)
        gCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTcc                                                                                                                                             <  6:202778/149‑1 (MQ=255)
            cccATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                         <  4:66602/149‑1 (MQ=255)
                     gAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATcgc                                                                                                                                 <  6:241631/148‑1 (MQ=255)
                              aGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTc                                                                                                                       <  2:232208/149‑1 (MQ=255)
                                 cAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAgg                                                                                                                    <  1:359432/149‑1 (MQ=255)
                                      cgtgcgtgATCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                               <  3:5476/148‑1 (MQ=255)
                                      cgggcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                               <  6:297392/145‑1 (MQ=255)
                                         gcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAg                                                                                                            >  8:47906/1‑148 (MQ=255)
                                                    gCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGtt                                                                                                  <  4:310174/147‑1 (MQ=255)
                                                    gCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGtt                                                                                                  >  3:310174/1‑147 (MQ=255)
                                                    gCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgc                                                                                                >  5:358061/1‑149 (MQ=255)
                                                           cccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTAGCACAGTGATCGCGCAGCCGTTAAGGCTTTGAGCAAGAGAGTTGCTGCca                                                                                           >  5:204959/1‑147 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAg                                                                                        >  4:13538/1‑149 (MQ=255)
                                                               gcaaaCTTCTCGTAGTAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGt                                                                                     >  5:7071/4‑149 (MQ=39)
                                                               gCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGt                                                                                     <  7:145689/149‑1 (MQ=255)
                                                                ccAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtc                                                                                    >  5:343307/1‑149 (MQ=255)
                                                                  aaCTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctc                                                                                  >  2:357666/1‑149 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcg                                                                               <  3:5277/149‑1 (MQ=255)
                                                                          ctacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTCGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCACTCGCGAAc                                                                          >  7:160867/1‑149 (MQ=255)
                                                                              gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                        >  2:74905/1‑147 (MQ=255)
                                                                              gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                        >  1:200960/1‑147 (MQ=255)
                                                                               agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                        >  7:45612/1‑146 (MQ=255)
                                                                                 aaGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTCTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGgaga                                                                    <  1:74905/149‑1 (MQ=255)
                                                                                 aaGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGAt                                                                   <  6:358061/149‑1 (MQ=255)
                                                                                 aaGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGAt                                                                   >  4:136546/1‑149 (MQ=255)
                                                                                   gTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGAttt                                                                 <  7:47906/149‑1 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                              <  3:95388/135‑1 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                              >  4:95388/1‑135 (MQ=255)
                                                                                         ggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACatat                                                           >  6:72966/1‑149 (MQ=255)
                                                                                         ggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACatat                                                           <  3:136546/149‑1 (MQ=255)
                                                                                           ctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATg                                                          >  6:244088/1‑148 (MQ=255)
                                                                                            tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaa                                                        >  8:45370/1‑149 (MQ=255)
                                                                                              cggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaa                                                        <  6:64911/147‑1 (MQ=255)
                                                                                                gcCTGACGCTCCCTGATTCAATGCAATACGCCGTCGCTGTGTGAGCGGGCGAGCCCTCTCACACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTGTCGCGAACTGGAGATTTACATATGaacaac                                                    <  8:160867/149‑1 (MQ=255)
                                                                                                 ccTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca                                                   <  6:95449/149‑1 (MQ=255)
                                                                                                       gCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTc                                              >  7:13371/1‑148 (MQ=255)
                                                                                                         tCCCTGATTCATATCAATACGAGGACCCTGTGGGAGCGGGCGAGCCGGCTCCGACAGTGATCGCGCAGCGTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtc                                                                                    <  2:21600/108‑1 (MQ=255)
                                                                                                         tCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGAACCCCCAACCTTTCAGGCAATGACCAAGAGAGTTGCAGCCACAgggac                                                                                    >  1:21600/1‑106 (MQ=37)
                                                                                                           ccTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaa                                           <  8:13371/147‑1 (MQ=255)
                                                                                                            cTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaagg                                         >  1:295399/1‑148 (MQ=255)
                                                                                                             tGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                       >  5:369517/1‑101 (MQ=255)
                                                                                                             tGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                       <  6:369517/101‑1 (MQ=255)
                                                                                                                ttCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaagg                                    >  6:356651/1‑149 (MQ=255)
                                                                                                                   aaaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTg                                 >  7:8587/1‑149 (MQ=255)
                                                                                                                   aaaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTg                                 >  2:199300/1‑149 (MQ=255)
                                                                                                                   aaaGCAATAAGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAATGTTg                                 >  5:194042/1‑149 (MQ=255)
                                                                                                                     aGCATTACTCCGACCCTGTTGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGTCTTTGAGCAAGACAGTTGCTGCCACAGGGTCTCTCGCTTACTGGAGATTTAGATATGATCAACAGATTCCTAGGCAATGTTGcc                               <  6:204959/149‑1 (MQ=37)
                                                                                                                     aGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGcc                               <  2:295399/149‑1 (MQ=255)
                                                                                                                      gCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCGACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCATGAGAGTTGCTGCCACAGGGTCTCTCGCGAAGTGGAGATTTACTTATGAACAACAGATTCCAAGGCAAGGTTGccc                              >  5:277802/1‑149 (MQ=255)
                                                                                                                        aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                            <  5:356651/149‑1 (MQ=255)
                                                                                                                        aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                            <  1:357666/149‑1 (MQ=255)
                                                                                                                        aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                            <  6:7071/149‑1 (MQ=255)
                                                                                                                         aTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTgg                           >  7:199326/1‑149 (MQ=255)
                                                                                                                          tACGCCGACACTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTcc                                             <  3:167811/130‑1 (MQ=255)
                                                                                                                          tACGCCGACACTGTGGGAGCGGGCGAGCCCGCTCCCACACTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTcc                                             >  4:167811/1‑130 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                          >  6:101061/1‑82 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                          >  6:101057/1‑82 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                          <  5:101057/82‑1 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                          <  5:101061/82‑1 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca                                                   >  1:251441/1‑121 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca                                                   <  2:251441/121‑1 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCAc                       <  6:343307/149‑1 (MQ=255)
                                                                                                                                atcacTGTGGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGATTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCgg                    >  1:273586/5‑149 (MQ=255)
                                                                                                                                   ccTGTGGGAGCGGGCGAGCCCGCTCCCCCAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCGCTGGTCACCGGTgc                 >  5:244437/1‑149 (MQ=255)
                                                                                                                                    cTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgcc                >  1:143924/1‑149 (MQ=255)
                                                                                                                                    cTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgcc                >  6:176259/1‑149 (MQ=255)
                                                                                                                                     tgtgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGa                                                 >  8:249359/1‑115 (MQ=255)
                                                                                                                                     tgtgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGa                                                 <  7:249359/115‑1 (MQ=255)
                                                                                                                                      gtgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgc              <  1:195248/149‑1 (MQ=255)
                                                                                                                                       tgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgcc             <  4:131075/149‑1 (MQ=255)
                                                                                                                                        gggAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgcc             >  5:270057/1‑148 (MQ=255)
                                                                                                                                          gAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGccc            >  4:330126/1‑147 (MQ=255)
                                                                                                                                           aGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaagg                                    <  5:210505/122‑1 (MQ=255)
                                                                                                                                           aGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAg         <  3:63412/149‑1 (MQ=255)
                                                                                                                                           aGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAg         <  2:200960/149‑1 (MQ=255)
                                                                                                                                              gggCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaagg                                         <  1:247056/114‑1 (MQ=255)
                                                                                                                                              gggCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaagg                                         >  2:247056/1‑114 (MQ=255)
                                                                                                                                                 cGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACTGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCg   <  3:13538/149‑1 (MQ=255)
                                                                                                                                                  gAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCgg  >  1:41340/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCGG  >  NC_002947/3584752‑3585046

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: