Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,581,991 A→C 5.7% Y21S (TAC→TCC)  benA → benzoate 1,2‑dioxygenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,581,9910AC 5.7% 226.5 / 3.9 70Y21S (TAC→TCC) benAbenzoate 1,2‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base C (2/2);  total (39/31)
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.73e-01

CGTTCTGCACAATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAG  >  NC_002947/3581847‑3582136
                                                                                                                                                 |                                                                                                                                                 
cGTTCTGCACAATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGc                                                                                                                                               >  7:97467/1‑149 (MQ=255)
  ttCTGCACAATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTg                                                                                                                                             >  5:335676/1‑149 (MQ=255)
      gCACAATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAg                                                                                                                                         >  7:311587/1‑149 (MQ=255)
       cacaATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGc                                                                                                                                        <  1:142952/149‑1 (MQ=255)
           aTCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCgg                                                                                                                                      >  1:359592/1‑147 (MQ=255)
                    gCCCCCCGCCAGCCGGCTCCCTAACCTGACCACGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCAc                                                                                                                           <  4:103598/149‑1 (MQ=255)
                      cccccGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCAcc                                                                                                                          >  8:283959/1‑148 (MQ=255)
                       ccccGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGa                                                                                                                        >  6:187435/1‑149 (MQ=255)
                            ccAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCAc                                                                                                                   <  8:97467/149‑1 (MQ=255)
                             cAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACg                                                                                                                  >  6:77200/1‑149 (MQ=255)
                             cAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGACATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACg                                                                                                                  >  7:167206/1‑149 (MQ=255)
                                  gTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGt                                                                                                             >  2:169830/1‑149 (MQ=255)
                                       ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCa                                                                                                                            <  4:48936/129‑1 (MQ=255)
                                       ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCa                                                                                                                            >  3:48936/1‑129 (MQ=255)
                                       ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGAATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGAc                                                                                                        >  8:285487/1‑149 (MQ=255)
                                        cTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGAcc                                                                                                       >  3:142867/1‑149 (MQ=255)
                                             cTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGATATGTTCACCGACCCACGGCTGTTCGACCTg                                                                                                     >  5:89420/1‑146 (MQ=255)
                                              tGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGtt                                                                                                            >  5:241131/1‑138 (MQ=255)
                                              tGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGtt                                                                                                            <  6:241131/138‑1 (MQ=255)
                                                aCCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGgaga                                                                                                 <  6:335676/147‑1 (MQ=255)
                                                aCCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATg                                                                                               <  4:142867/149‑1 (MQ=255)
                                                 ccAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGa                                                                                              <  2:359592/149‑1 (MQ=255)
                                                  cAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGCCCTGGAGATGa                                                                                              >  8:294129/1‑148 (MQ=255)
                                                   aGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGaaa                                                                                            <  8:283922/149‑1 (MQ=255)
                                                     gTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTCCCGCttgcag                                                                                                                                         <  6:299893/102‑6 (MQ=255)
                                                     gTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTCCCGCttgcag                                                                                                                                         >  5:299893/1‑97 (MQ=255)
                                                        tAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAAcaca                                                                                        >  7:310863/1‑148 (MQ=255)
                                                          aaCAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGgaga                                                                                                                                  >  8:52694/1‑104 (MQ=255)
                                                          aaCAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCt                                                                                     >  3:93954/1‑149 (MQ=255)
                                                           aCAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCtt                                                                                    >  1:165286/1‑149 (MQ=255)
                                                               ttaCAAGGGATACCCTGGCCATGTCCCTGGGATTCTACTACATCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAg                                                                                <  8:167206/147‑1 (MQ=255)
                                                               tAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAg                                                                                <  7:283959/149‑1 (MQ=255)
                                                                 aCAAGGGAGACCCTGGCCATGTCCCTGGGAGTTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCAACGGCTGTTCGACCTGGAGATGAAaca                                                                                          >  1:33690/1‑138 (MQ=255)
                                                                 aCAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGTCATCTACCGCTGCAAGCGGGATATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAgg                                                                               <  6:89420/148‑1 (MQ=255)
                                                                 aCAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAaca                                                                                          <  2:33690/137‑1 (MQ=255)
                                                                 aCAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAgg                                                                               <  4:93954/148‑1 (MQ=255)
                                                                    aGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGAGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCaa                                                                           >  1:292348/1‑149 (MQ=255)
                                                                            ccTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATc                                                                    >  7:299579/1‑148 (MQ=255)
                                                                                 ccATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTg                                                              <  8:311587/149‑1 (MQ=255)
                                                                                        ccTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCAc                                                        <  4:124931/148‑1 (MQ=255)
                                                                                        ccTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCAc                                                        >  3:124931/1‑148 (MQ=255)
                                                                                         cTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAgg                                                                               >  7:44224/1‑124 (MQ=255)
                                                                                         cTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAgg                                                                               <  8:44224/124‑1 (MQ=255)
                                                                                          tGGGATTCGACTACGTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGa                                                      >  7:323891/1‑148 (MQ=255)
                                                                                          tGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGa                                                                                 >  6:111113/1‑121 (MQ=255)
                                                                                          tGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGa                                                                                 <  5:111113/121‑1 (MQ=255)
                                                                                           gggATTCGACTGCCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAg                                                                                >  3:63555/1‑121 (MQ=255)
                                                                                           gggATTCGACTGCCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAg                                                                                <  4:63555/121‑1 (MQ=255)
                                                                                            ggATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGaaa                                                    <  7:285487/148‑1 (MQ=255)
                                                                                                ttCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCa                                                <  5:77200/149‑1 (MQ=255)
                                                                                                    aCTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATc                                            <  1:169830/149‑1 (MQ=255)
                                                                                                                      tCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTc                          >  4:181736/1‑149 (MQ=255)
                                                                                                                         aGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAGAGCCAGATCCCCGAGAAGAATGACTTCCt                        >  7:173039/1‑148 (MQ=255)
                                                                                                                           gacgacCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCAc                     >  4:162233/1‑149 (MQ=255)
                                                                                                                              gacCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACt                            >  7:255428/1‑139 (MQ=255)
                                                                                                                              gacCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACt                            <  8:255428/139‑1 (MQ=255)
                                                                                                                                ccGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGa                <  8:299579/149‑1 (MQ=255)
                                                                                                                                 cGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGAc               >  5:171602/1‑149 (MQ=255)
                                                                                                                                      aaaGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATg           <  2:165286/148‑1 (MQ=255)
                                                                                                                                       gtgaacATCTCCCGCTTGCAGCGGTAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGccc                                                          >  6:159474/6‑100 (MQ=37)
                                                                                                                                       gtgaacATCTCCCGCTTGCAGCGGTAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGccc                                                          <  5:159474/95‑1 (MQ=37)
                                                                                                                                        aGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATgggg        <  3:162233/149‑1 (MQ=255)
                                                                                                                                         gggCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCgaga                                     >  8:220306/1‑119 (MQ=255)
                                                                                                                                         gggCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCgaga                                     <  7:220306/119‑1 (MQ=255)
                                                                                                                                         gggAATCTACCGCTGCAAGCGGGATATGTTCACCGACCCACGGCTTTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGAGCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGgc       <  2:292348/149‑1 (MQ=255)
                                                                                                                                          ggCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGACAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGgcg      >  3:217745/1‑149 (MQ=255)
                                                                                                                                           gCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATATACCTGGCCCACGAAATCCAGATCCCCGAGAAGAATGACTTCATCACCCTGACCATGCGgcaa     >  4:59999/1‑147 (MQ=255)
                                                                                                                                              tCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAg  >  3:174574/1‑149 (MQ=255)
                                                                                                                                               cTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGTCCCACGAAAGACAGATCCCCGAGAAGAATGACTTACTCACCCTGACCATGGGGAGcc    >  2:235649/1‑146 (MQ=255)
                                                                                                                                               cTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGcc    <  1:235649/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGTTCTGCACAATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCATGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAG  >  NC_002947/3581847‑3582136

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: