Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,994 C→G 7.3% S112R (AGC→AGG nemA → chromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,9940CG 7.3% 206.3 / 8.0 69S112R (AGC→AGGnemAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme
Reads supporting (aligned to +/- strand):  ref base C (37/27);  new base G (2/3);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

CTGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAA  >  NC_000913/1726846‑1727143
                                                                                                                                                    |                                                                                                                                                     
cTGCCCAGGCCAAAGGATATGCATGTGCGCCTTGCATCCATAGGCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCCGGCGTTCATGCTGAAAAAGGTCATATGGCCGTTCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGcc                                                                                                                                                      <  1:1203876/150‑1 (MQ=255)
 tGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGATAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCt                                                                                                                                                     >  1:1154072/1‑150 (MQ=255)
     cAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCaa                                                                                                                                                 >  1:22390/1‑150 (MQ=255)
      aGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAAc                                                                                                                                                >  1:1009118/1‑150 (MQ=255)
         cAAAAGGATACGCAGGTGCGCCCGGCATCCATAGTCCGGAGCAAATAGCCGCATGGAACAAAATCACCGCGGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTg                                                                                                                                             <  2:841168/150‑1 (MQ=255)
           aaaGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGc                                                                                                                                           <  2:257547/150‑1 (MQ=255)
              ggATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtg                                                                                                                                        <  1:437773/150‑3 (MQ=255)
              ggAGATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtg                                                                                                                                        <  2:1223997/150‑3 (MQ=255)
               gATATGCAGGTGCGCCTGGCATGCATAGTCCGGCGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGGTTGCAGGCTGGCgtga                                                                                                                                       <  2:499244/150‑4 (MQ=255)
               gATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAATATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTc                                                                                                                                       >  1:916097/1‑150 (MQ=255)
                         tGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGt                                                                                                                             <  1:503529/150‑1 (MQ=255)
                            gcCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCGGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAgc                                                                                                                          <  1:463083/150‑1 (MQ=255)
                            gcCGGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCGGGCGATCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAgc                                                                                                                          <  1:458498/150‑1 (MQ=255)
                             ccTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCCGGCTCCGGTAgcg                                                                                                                         >  2:1302886/1‑150 (MQ=255)
                                   agccATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCa                                                                                                                   <  1:814392/148‑1 (MQ=255)
                                          gTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTa                                                                                                            >  1:26583/1‑150 (MQ=255)
                                            ccGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgc                                                                                                          <  1:394168/150‑1 (MQ=255)
                                              ggAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAgcgc                                                                                                        <  1:793598/150‑1 (MQ=255)
                                               gAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTcgcgcg                                                                                                       <  2:546205/150‑1 (MQ=255)
                                               gAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGcgcg                                                                                                       <  2:132964/150‑1 (MQ=255)
                                                aGCAACTTGCCGCATCGAAAAAAATCACCGCTGGCCTTCATGCTGAACATGCTCATCGGGCCGCCCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTACCGCgg                                                                                                      <  1:768876/150‑1 (MQ=255)
                                                  cAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGa                                                                                                    >  2:1116681/1‑150 (MQ=255)
                                                        gCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCCCCGGTAGCGCCGGCAGCACTTAGCGCGGGAACACGt                                                                                              >  1:49068/1‑150 (MQ=255)
                                                           gCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACt                                                                                           >  1:453247/1‑150 (MQ=255)
                                                           gCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACt                                                                                           >  2:700117/1‑150 (MQ=255)
                                                             aTGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTtc                                                                                         <  1:899771/150‑1 (MQ=255)
                                                             aTGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGGCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTtc                                                                                         <  1:1155785/150‑1 (MQ=255)
                                                                    aaaaTCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTgcgc                                                                                  >  1:462330/1‑150 (MQ=255)
                                                                      aaTCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGa                                                                                >  2:944741/1‑150 (MQ=255)
                                                                            cGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACCCGTACTTCTCTGCGCGATGaaaa                                                                          >  1:948498/1‑150 (MQ=255)
                                                                              cTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATg                                                                        >  2:812029/1‑150 (MQ=255)
                                                                              cTGGCGGTCATGCTGAAACTGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATg                                                                        <  1:987877/150‑1 (MQ=255)
                                                                                ggCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGt                                                                      <  1:904956/150‑1 (MQ=255)
                                                                                ggCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGt                                                                      <  1:909808/150‑1 (MQ=255)
                                                                                       aTGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCCGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATTGTCAGGCGa                                                               >  2:828949/1‑150 (MQ=255)
                                                                                       aTGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTTGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGa                                                               >  1:1140866/1‑150 (MQ=255)
                                                                                        tGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCACCCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGAt                                                              >  2:793246/1‑150 (MQ=255)
                                                                                        tGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGAt                                                              >  2:652377/1‑150 (MQ=255)
                                                                                          cTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATcc                                                            >  1:319091/1‑150 (MQ=255)
                                                                                             aaaaTGGGCATATGGGCGTGAAGATGTGGCACACCGGACGCATTTCTCATGCCAGCCTGCAACCTGTCGGTCTGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTTCTTCTCTGCGCGATGAAAATGGTCAGGCGATCCgtg                                                         <  1:973080/150‑1 (MQ=255)
                                                                                              aaaTGGTCATATGGCCGTTCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCgtgt                                                        <  2:1100441/150‑1 (MQ=255)
                                                                                                aTGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAATATGGTCAGGCGATCCGTGTTg                                                      >  1:1089068/1‑150 (MQ=255)
                                                                                                  ggTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTGTGCCCGATGAAAATGGTCAGGCGATCCGTGTTGaa                                                    >  1:400832/1‑150 (MQ=255)
                                                                                                     cATAAGGCCGTCCAGCTGGGGAACACCGGACGCATTTCTCACGCCAGCCTTCAACCTGTCGGTCAGGCACCGGTAGCGCCTTCCGCACTTAGCGCGGGAACACGTACTGCTCTGCGCGATGAAAATGGTCAGTCGATCCGTGTTGAcaca                                                 <  2:939388/150‑1 (MQ=255)
                                                                                                           gCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATg                                           <  1:461860/150‑1 (MQ=255)
                                                                                                            ccGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCCGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCGTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGCTCCGTGTTGAAACATCCATGc                                          >  1:27144/1‑150 (MQ=255)
                                                                                                             cGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGcc                                         <  2:677260/150‑1 (MQ=255)
                                                                                                                gCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCcgcg                                      >  1:900012/1‑150 (MQ=255)
                                                                                                                        ggCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTg                              <  1:733004/150‑1 (MQ=255)
                                                                                                                        ggCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTg                              >  2:1037826/1‑150 (MQ=255)
                                                                                                                          cacaCCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGGGCGCTAGaa                            >  1:853739/1‑150 (MQ=255)
                                                                                                                          cacaCCGGACGCATATCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTGCAGCACTTAGCGCGCGAACACGTACTTATCTGCGCGATGAAAATGGTCAGGCGATCCGTGAGGAAACATCCAAGCCGCGTGCGCTTGaa                            >  2:1284484/1‑150 (MQ=255)
                                                                                                                              ccGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTGCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGGGTGGAAACATCCATGCCGCGTGCGCTTGACCTgg                        >  1:621266/1‑150 (MQ=255)
                                                                                                                                   cGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAagag                   <  1:1116681/150‑1 (MQ=255)
                                                                                                                                     cATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAt                 <  1:63516/150‑1 (MQ=255)
                                                                                                                                     cATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAt                 <  1:548851/150‑1 (MQ=255)
                                                                                                                                      aTTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCCGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAtt                >  2:35137/1‑150 (MQ=255)
                                                                                                                                      aTTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGCTCCGTGTTGAAACATCCCTGCCGCGTGCGCTTGAACTGGAAGAGAtt                >  1:787508/1‑150 (MQ=255)
                                                                                                                                      aTTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAtt                >  2:35136/1‑150 (MQ=255)
                                                                                                                                       tttCTCACGGCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCACGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGTACTGGAAGAGATTc               <  2:722877/150‑1 (MQ=255)
                                                                                                                                        ttCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCCGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAGTcc              >  2:834258/1‑150 (MQ=255)
                                                                                                                                         ctgaccGCCAGGTTGCAGGCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGTGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCa             >  1:1273008/6‑150 (MQ=255)
                                                                                                                                          tgaccGCCAGGTTGCAGGCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAg            >  1:902028/5‑150 (MQ=255)
                                                                                                                                           tcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGTACACGTTCTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCTTGCCGCGTGCGCTTGACCTGGAAGAGATTCCAgg           >  2:646416/1‑150 (MQ=255)
                                                                                                                                           tcACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCAATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAgg           >  2:642543/1‑150 (MQ=255)
                                                                                                                                              cGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCGGCGCGATGAAAATGGTCAGGCGATCCGGGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCAAGGTAt        >  1:665181/1‑150 (MQ=255)
                                                                                                                                                 ccgCCTGCAACCTGGCGGTCAGGCACCGGTAGAGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGt     >  2:261562/3‑150 (MQ=255)
                                                                                                                                                  aGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTc    >  1:1260577/1‑150 (MQ=255)
                                                                                                                                                  aGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTc    <  2:1017819/150‑1 (MQ=255)
                                                                                                                                                   gCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCa   >  1:716402/1‑150 (MQ=255)
                                                                                                                                                    ccTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCaa  >  1:1065461/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCCCAGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAA  >  NC_000913/1726846‑1727143

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