Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,291 T→A 7.3% intergenic (+40/‑156) 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,2910TA 7.3% 230.1 / 6.0 82intergenic (+40/‑156)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 T (39/37);  new base A (4/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGG  >  NC_000913/2323142‑2323439
                                                                                                                                                     |                                                                                                                                                    
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: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:406496/150‑1 (MQ=255)
                                                                                                                  aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTg                                    >  1:1089049/1‑150 (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:597264/150‑1 (MQ=255)
                                                                                                                                           caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTAt           >  2:654200/1‑150 (MQ=255)
                                                                                                                                           caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAGTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGGTATTGCCTGACTGTACCCAGAACGGTGTAt           >  2:841959/1‑150 (MQ=255)
                                                                                                                                             gTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTCATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGAACCCACAACGGTGTATGc         >  2:748243/1‑150 (MQ=255)
                                                                                                                                                    aTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAggg  <  1:673233/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGG  >  NC_000913/2323142‑2323439

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: