Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,727,001 C→G 7.5% P115A (CCT→GCT)  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,727,0010CG 7.5% 204.0 / 8.1 67P115A (CCT→GCT) nemAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base G (2/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

AGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGATTTC  >  NC_000913/1726852‑1727150
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                  >  2:700117/1‑150 (MQ=255)
                                                     gCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACt                                                                                                  >  1:453247/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:909808/150‑1 (MQ=255)
                                                                          ggCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGt                                                                             <  1:904956/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                                     >  2:1037826/1‑150 (MQ=255)
                                                                                                                  ggCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTg                                     <  1:733004/150‑1 (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:548851/150‑1 (MQ=255)
                                                                                                                               cATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGAt                        <  1:63516/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)
                                                                                                                                                     ccTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGATTTc  <  2:22390/150‑1 (MQ=255)
                                                                                                                                                     ccTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTAAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAAACATTTc  >  2:588756/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCAAAAGGATATGCAGGTGCGCCTGGCATCCATAGTCCGGAGCAAATTGCCGCATGGAAAAAAATCACCGCTGGCGTTCATGCTGAAAATGGTCATATGGCCGTGCAGCTGTGGCACACCGGACGCATTTCTCACGCCAGCCTGCAACCTGGCGGTCAGGCACCGGTAGCGCCTTCAGCACTTAGCGCGGGAACACGTACTTCTCTGCGCGATGAAAATGGTCAGGCGATCCGTGTTGAAACATCCATGCCGCGTGCGCTTGAACTGGAAGAGATTCCAGGTATCGTCAATGATTTC  >  NC_000913/1726852‑1727150

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: