Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,314 T→C 8.0% F284S (TTT→TCT)  PP_1478 → xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,684,3140TC 8.0% 183.5 / 4.4 62F284S (TTT→TCT) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base C (3/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GACTGGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACT  >  NC_002947/1684166‑1684457
                                                                                                                                                    |                                                                                                                                                
gactgGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGACCCTGGCTATCAGACGCGCttt                                                                                                                                                 >  2:1329985/1‑150 (MQ=255)
  ctgGGTCGAGGACGGCTGGAACCCGGATGAATCCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGc                                                                                                                                               >  1:299796/1‑150 (MQ=255)
  ctgGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGc                                                                                                                                               <  1:1383862/150‑1 (MQ=255)
  ctgGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGGGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGc                                                                                                                                               >  2:1066333/1‑150 (MQ=255)
                gCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCa                                                                                                                                 >  1:156092/1‑150 (MQ=255)
                 cTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                <  2:156092/150‑1 (MQ=255)
                   ggAACCCGGATGAAACCGTCGAACTGGCGCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAgg                                                                                                                              <  1:290917/150‑1 (MQ=255)
                    gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                     <  1:852795/142‑2 (MQ=255)
                    gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                     >  2:852795/1‑141 (MQ=255)
                      aCCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAg                                                                                                                           >  1:286712/1‑150 (MQ=255)
                        ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCttt                                                                                                                                                 >  2:1357743/1‑126 (MQ=255)
                        ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCttt                                                                                                                                                 <  1:1357743/126‑1 (MQ=255)
                        ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                         >  2:434049/1‑150 (MQ=255)
                                   cGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAc                                                                                                             >  1:1588469/1‑150 (MQ=255)
                                     tCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCg                                                                                                           <  1:1216496/150‑1 (MQ=255)
                                       gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGCGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGAAGAGCGCGTGCGCAAGGAAGCaga                                                                                                                       >  1:1514969/1‑137 (MQ=255)
                                       gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGCGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGAAGAGCGCGTGCGCAAGGAAGCaga                                                                                                                       <  2:1514969/137‑1 (MQ=255)
                                         aCTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCAc                                                                                                       >  1:979049/1‑150 (MQ=255)
                                              cTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTAGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCg                                                                                                  >  2:117305/1‑150 (MQ=255)
                                               tcgtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGCGTTTGAGGAAGCAGAGATCGCCACCGGCACCGTCgg                                                                                                 <  1:1365475/150‑1 (MQ=255)
                                                 gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCggaa                                                                                               <  2:1588469/150‑3 (MQ=255)
                                                    gCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGa                                                                                            >  1:567024/1‑150 (MQ=255)
                                                     ccTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGAt                                                                                           >  1:636279/1‑150 (MQ=255)
                                                      cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAg                                                                                                                               <  2:1350601/114‑1 (MQ=255)
                                                      cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGCCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAg                                                                                                                               >  1:1350601/1‑114 (MQ=255)
                                                           ggTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGa                                                                                     <  2:567024/150‑1 (MQ=255)
                                                                          ccTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAAGCCCGGGGGAGCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGa                                                                      <  2:631835/150‑1 (MQ=255)
                                                                             gATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                   >  1:363237/1‑150 (MQ=255)
                                                                                  aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAAc                                                                                   <  1:377084/129‑1 (MQ=255)
                                                                                  aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAAc                                                                                   >  2:377084/1‑129 (MQ=255)
                                                                                       tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGt                                                         >  1:1462299/1‑150 (MQ=255)
                                                                                         ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGAACGGAAACAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                       >  1:690816/1‑150 (MQ=255)
                                                                                         ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACACGCTTTGCAGAACGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATAATCACTGAACCTGCCCAGACCGAGCACATTCTGCGTAc                                                       >  2:241651/1‑150 (MQ=255)
                                                                                           ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTg                                                            >  2:952571/1‑143 (MQ=255)
                                                                                           ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTg                                                            <  1:952571/143‑1 (MQ=255)
                                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                >  2:1632932/1‑66 (MQ=255)
                                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                <  1:1632932/66‑1 (MQ=255)
                                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCag                                                                                                                        >  1:792612/1‑74 (MQ=255)
                                                                                                     cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCag                                                                                                                        <  2:792612/74‑1 (MQ=255)
                                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGAc                                          <  1:434049/150‑1 (MQ=255)
                                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTACCCAGGCCGa                                                                      >  1:479532/1‑122 (MQ=255)
                                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTACCCAGGCCGa                                                                      <  2:479532/122‑1 (MQ=255)
                                                                                                      tCGGTCAACGCCGAAATTCCTACAGGCCATGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGGAAAGATCGAACACCGCACCGGCGGTATTATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACGGGGCAGGCCGAc                                          >  1:372207/1‑150 (MQ=255)
                                                                                                                ccGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTcc                                >  2:1576205/1‑150 (MQ=255)
                                                                                                                 cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                   >  2:13242/1‑67 (MQ=255)
                                                                                                                 cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                   <  1:13242/67‑1 (MQ=255)
                                                                                                                  gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                        <  1:1582311/92‑1 (MQ=255)
                                                                                                                  gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                        >  2:1582311/1‑92 (MQ=255)
                                                                                                                     aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCt                           >  1:1240766/1‑150 (MQ=255)
                                                                                                                         cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                           <  1:995414/114‑1 (MQ=255)
                                                                                                                         cTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                           >  2:995414/1‑114 (MQ=255)
                                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                         >  2:300911/1‑131 (MQ=255)
                                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                         <  1:300911/131‑1 (MQ=255)
                                                                                                                          tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                      >  2:794607/1‑150 (MQ=255)
                                                                                                                                cccTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACCTTCTGCTTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctgct                >  1:351064/1‑150 (MQ=255)
                                                                                                                                      cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTc                                                              <  2:126694/98‑1 (MQ=255)
                                                                                                                                      cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTc                                                              >  1:126694/1‑98 (MQ=255)
                                                                                                                                         tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                       >  1:605531/1‑102 (MQ=255)
                                                                                                                                         tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                       <  2:605531/102‑1 (MQ=255)
                                                                                                                                            gACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGTAGCTGCTGCGCGATCCTTa    >  2:1076474/1‑150 (MQ=255)
                                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                    >  2:385555/2‑70 (MQ=38)
                                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                    <  1:385555/69‑1 (MQ=38)
                                                                                                                                            caCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTa    >  1:481807/2‑150 (MQ=255)
                                                                                                                                             aCGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTAc   <  2:979049/150‑1 (MQ=255)
                                                                                                                                              cgcgCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTCATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACt  <  1:1329985/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GACTGGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACT  >  NC_002947/1684166‑1684457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: