Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,283 G→C 6.7% intergenic (+32/‑164) 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,2830GC 6.7% 260.8 / 10.8 89intergenic (+32/‑164)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 G (44/39);  new base C (4/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

AGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGC  >  NC_000913/2323137‑2323432
                                                                                                                                                  |                                                                                                                                                     
aGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                                                                                    >  1:721270/1‑150 (MQ=255)
aGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                                                                                    >  2:673233/1‑150 (MQ=255)
 cggCTGGCACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATCCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcg                                                                                                                                                   <  2:183568/148‑1 (MQ=255)
  tGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                                  <  2:265645/150‑1 (MQ=255)
  tGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                                                                  <  1:1194017/150‑1 (MQ=255)
     tGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                                                               >  2:387665/1‑150 (MQ=255)
      ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                                                              >  1:443367/1‑150 (MQ=255)
      ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                                                              >  1:443769/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                                                                                                                       <  1:512705/150‑1 (MQ=255)
                             cccGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                       <  2:969503/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:968698/150‑1 (MQ=255)
                                                           gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                                         <  1:979569/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)
                                                                                                                                                caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTAt    >  2:654200/1‑150 (MQ=255)
                                                                                                                                                caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTAt    <  2:597264/150‑1 (MQ=255)
                                                                                                                                                caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAGTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGGTATTGCCTGACTGTACCCAGAACGGTGTAt    >  2:841959/1‑150 (MQ=255)
                                                                                                                                                  gTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTCATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGAACCCACAACGGTGTATGc  >  2:748243/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGC  >  NC_000913/2323137‑2323432

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: