Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,280 A→G 7.1% intergenic (+29/‑167) atoC → / → atoD fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,323,2800AG 7.1% 248.3 / 7.8 86intergenic (+29/‑167)atoC/atoDfused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base G (4/2);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTG  >  NC_000913/2323132‑2323427
                                                                                                                                                    |                                                                                                                                                   
aTGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcac                                                                                                                                                    <  1:386416/150‑1 (MQ=255)
aTGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcac                                                                                                                                                    >  2:583081/1‑150 (MQ=255)
     aGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                                                                               >  2:673233/1‑150 (MQ=255)
     aGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                                                                               >  1:721270/1‑150 (MQ=255)
      cggCTGGCACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATCCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcg                                                                                                                                              <  2:183568/148‑1 (MQ=255)
       tGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                             <  1:1194017/150‑1 (MQ=255)
       tGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                             <  2:265645/150‑1 (MQ=255)
          tGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                                                          >  2:387665/1‑150 (MQ=255)
           ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                                                         >  1:443769/1‑150 (MQ=255)
           ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                                                         >  1:443367/1‑150 (MQ=255)
           ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATACGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCCCAAtt                                                                                                                                         >  2:501181/1‑150 (MQ=255)
                 acaaGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGc                                                                                                                                   <  1:1216513/150‑1 (MQ=255)
                 acaaGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                                                                   <  1:1286512/150‑1 (MQ=39)
                   tggagtttatacatcaatgcacggcgactgatgcccagcattaaagcagtgcATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACGTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                                                                 >  1:96018‑M2/53‑150 (MQ=255)
                     gaagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                                                               <  1:820416/150‑1 (MQ=39)
                      aagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                                                                              <  1:1199491/150‑1 (MQ=255)
                      aagGAAACCGAACCCGCACTGCATTAATGCTGGGAATCAGTCGCCGTTCATTGATGTATAAACTCCAGGAATACGTTATCGATACGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                                                                              <  2:46553/150‑1 (MQ=255)
                      aagGAAACCGAACCCGCACAGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                                                                              <  2:1105945/150‑1 (MQ=255)
                       agGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                             >  2:516786/1‑150 (MQ=255)
                          aaaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGATATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGca                                                                                                                          >  1:482455/1‑149 (MQ=255)
                           aaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCAGTGCATAGATGTATAAACTCCAGGAATACTGTATCGAGCCGGCGGATGTATAACACAAAAACTTTCAATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCGGCTc                                                                                                                         >  2:253418/1‑150 (MQ=255)
                              catcaatgcacggcgactgatgcccagcattaaagcagtgcATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAtt                                                                                                                      >  2:1063742‑M2/42‑150 (MQ=255)
                                tcaatgcacggcgactgatgcccagcattaaagcagtgcATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                                                                    >  1:597264‑M2/40‑150 (MQ=255)
                                  cccGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                  <  2:969503/150‑1 (MQ=255)
                                  cccGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                  <  1:512705/150‑1 (MQ=255)
                                    cGCACTGCTTTAATTCTGGGCATCAGTCGCCGTTCATTGATGTATAAAATCCAGGAATACGGTATAGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                                <  2:626036/150‑1 (MQ=255)
                                    cGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                                >  1:773602/1‑150 (MQ=255)
                                       ACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTAGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                                                                             >  1:1298728‑M1/1‑150 (MQ=255)
                                       ACTGCTTTAATGCGGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                                                                             >  1:1125634‑M1/1‑150 (MQ=255)
                                                   cTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATACCCATACa                                                                                                 >  1:148948/1‑150 (MQ=255)
                                                     gggCATGATGCGCCGCGCCTTGATGTATAAACTCCAGGAATACTGTCTCGATCCGGCGGATGTATAACACCAAACCTTGCTATGCAGATATTTGCTGCGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCGATCCAAt                                                                                               <  2:135337/150‑1 (MQ=255)
                                                          tCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATc                                                                                          <  1:392894/150‑1 (MQ=255)
                                                           cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                         >  2:332660/1‑150 (MQ=255)
                                                                gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                                    <  1:979569/150‑1 (MQ=255)
                                                                gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                                    <  1:968698/150‑1 (MQ=255)
                                                                gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACAACAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAAACCCCATCCAATTTATCCCTTCa                                                                                    >  2:1172016/1‑150 (MQ=255)
                                                                  cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGCTCCGGCGGCTGGATCACACCAAAACTTGCTATGCAGACATTTGCACAGTGCGCCATTTTCTGCATAGCCGCTCATTCGCCTTATAAATCGCCTTACAATGTCTCCCTTCata                                                                                  >  1:206444/1‑150 (MQ=255)
                                                                  cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATGTTCTGCATAGCCGCTCATTCTCCTTATAAATACCCATCCAATTTCTCCCTTCata                                                                                  >  1:69923/1‑150 (MQ=255)
                                                                   gtgcATTGATGTATAAACTCCAGGAAGACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCACAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                                                                 >  1:966090‑M2/5‑150 (MQ=255)
                                                                    tgcATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAACTCCCCATCCAATTTATCCCTTCATAtt                                                                                >  2:1158370‑M2/4‑150 (MQ=255)
                                                                     gcATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                                                                               >  1:914709‑M2/3‑150 (MQ=255)
                                                                      cagtGATGTATAACCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGAAGAGCCGCGCATTCTCCTTTAAAATCCCGATACATTTTATCCCTTCatatgaa                                                                              >  2:75760/4‑147 (MQ=255)
                                                                        ttGATGTATAAACTCCAGGAATACTGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGAGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                            <  1:173042/150‑1 (MQ=255)
                                                                        ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATTCAGAAATTTGCACAGTGCGCAATATTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                            <  1:1171454/150‑1 (MQ=255)
                                                                        ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                            <  1:1176761/150‑1 (MQ=255)
                                                                         tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                                                           <  2:400792/150‑1 (MQ=255)
                                                                         tGATGTATAAACTACAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                                                           <  2:1038622/150‑1 (MQ=255)
                                                                          gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTa                                                                          >  1:470547/1‑150 (MQ=255)
                                                                                taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATCCAATTAGTTaaa                                                                    <  2:773602/150‑1 (MQ=255)
                                                                                  aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                                                                  >  2:888394/1‑150 (MQ=255)
                                                                                  aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAACCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                                                                  >  1:730046/1‑150 (MQ=255)
                                                                                    cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc                                                                >  2:943052/1‑150 (MQ=255)
                                                                                      ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                                                              <  1:270260/150‑1 (MQ=255)
                                                                                         ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTTTGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                                                           <  1:875070/150‑1 (MQ=255)
                                                                                           aaaactgTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCATCATATTCAATTAGTTAAATAACTAAATcc                                                         <  1:934171/144‑1 (MQ=255)
                                                                                            aTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa                                                        >  1:1301853/1‑150 (MQ=255)
                                                                                              aCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGAGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                                      <  2:951641/150‑1 (MQ=255)
                                                                                               cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCAATCCAATTTATCACTTCATATTCAATTAGTTAAATAACTAAATCCaata                                                     >  1:613677/1‑150 (MQ=255)
                                                                                                        tCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtta                                            >  2:971342/1‑149 (MQ=255)
                                                                                                         ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTTCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCt                                           >  2:884279/1‑150 (MQ=255)
                                                                                                         ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCt                                           <  2:148661/150‑1 (MQ=255)
                                                                                                         ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTACTGCATATCCGCTCATTCTACTAAAAAATCACCAACCAATTTATCCATACATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTca                                           >  1:797680/1‑149 (MQ=255)
                                                                                                          cggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCAAAGCCGCTCATTATCATTATAAAACCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATACAATAATCTCATTATg                                          >  2:984295/1‑150 (MQ=255)
                                                                                                              ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCAc                                      >  1:665170/1‑150 (MQ=255)
                                                                                                               gATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACt                                     >  2:1019653/1‑150 (MQ=255)
                                                                                                                     taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCtt                               <  1:1280256/150‑1 (MQ=255)
                                                                                                                           aaaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTAtt                         <  1:790685/149‑1 (MQ=255)
                                                                                                                            aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTg                        >  1:1089049/1‑150 (MQ=255)
                                                                                                                            aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTg                        <  1:406496/150‑1 (MQ=255)
                                                                                                                               aCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCAAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCt                     >  1:470194/1‑150 (MQ=255)
                                                                                                                                cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTACTTAAATAACTAAATCCAAAAATCTCATTCAGGCACTCCCCATGCTATTGCCTg                    >  2:1073192/1‑150 (MQ=255)
                                                                                                                                 ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGa                   >  2:1054039/1‑150 (MQ=255)
                                                                                                                                  tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGAc                  <  1:302847/150‑1 (MQ=255)
                                                                                                                                    cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTg                <  1:423988/150‑1 (MQ=255)
                                                                                                                                    cTATGCAGAAAATTGCACAGTGCGCTATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTATATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTg                <  2:705050/150‑1 (MQ=255)
                                                                                                                                     tATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGt               <  1:669808/150‑1 (MQ=255)
                                                                                                                                         cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCATTGCTATTGCCTGACTGTAccc           >  2:477858/1‑150 (MQ=255)
                                                                                                                                         cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCCCTCCCCTTGCTATTGCATGACTGTAcca           >  1:251162/1‑149 (MQ=255)
                                                                                                                                         cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTAccc           >  2:545168/1‑150 (MQ=255)
                                                                                                                                          aGAAATTTGCACAGTGCGCAATTTTCTTCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCca          >  2:506270/1‑150 (MQ=255)
                                                                                                                                          aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACAAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCca          <  2:550897/150‑1 (MQ=255)
                                                                                                                                          aGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCca          >  1:385237/1‑150 (MQ=255)
                                                                                                                                           gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCcac         <  2:584758/150‑1 (MQ=255)
                                                                                                                                            aaaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTTCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCcaca        <  2:267705/150‑1 (MQ=255)
                                                                                                                                              aTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATTCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAAc      <  2:938787/150‑1 (MQ=255)
                                                                                                                                              aTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAAc      <  2:614017/150‑1 (MQ=255)
                                                                                                                                               tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACg     >  1:918203/1‑150 (MQ=255)
                                                                                                                                                  gCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCAAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgtg  <  2:914709/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTG  >  NC_000913/2323132‑2323427

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: