Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,029,933 T→A 6.4% K347* (AAG→TAG) ‡ fadB ← dodecenoyl‑CoA delta‑isomerase, enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,029,9330TA 6.4% 400.1 / 12.2 124K347* (AAG→TAG) ‡fadBdodecenoyl‑CoA delta‑isomerase, enoyl‑CoA hydratase, 3‑hydroxybutyryl‑CoA epimerase, 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (51/65);  new base A (3/5);  total (54/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 128 total aligned reads displayed.

CACAATATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTACATATT  >  NC_000913/4029784‑4030082
                                                                                                                                                     |                                                                                                                                                     
cacaATATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACtt                                                                                                                                                       >  1:691052/1‑150 (MQ=255)
 acaATATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTg                                                                                                                                                      >  2:1169763/1‑150 (MQ=255)
 acaATATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTg                                                                                                                                                      <  2:870788/150‑1 (MQ=255)
   aaTATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTc                                                                                                                                                    <  1:443428/150‑1 (MQ=255)
       tCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTg                                                                                                                                                <  2:691052/150‑1 (MQ=255)
       tCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTg                                                                                                                                                >  1:960855/1‑150 (MQ=255)
        ccACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACAAGTCGttaa                                                                                                                                               >  1:894724/1‑148 (MQ=255)
         cACGCGGTCAAATCCGGCGTAGTCGAGCGTTTGGTTGATTGTGGAGATCACGCCAGCCAGTTTCAGAGCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGat                                                                                                                                              <  1:737503/150‑1 (MQ=255)
         cACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACAAGTCGttaac                                                                                                                                              <  2:442787/150‑4 (MQ=255)
         cACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACAAGTCGttaac                                                                                                                                              <  2:894724/150‑4 (MQ=255)
         cACGCGGTCAAATCCGGCATAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACAAGTCGttaac                                                                                                                                              <  2:1189518/150‑4 (MQ=255)
          aCGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACAAGTCGttaacc                                                                                                                                             <  2:887813/150‑5 (MQ=255)
                   aaTCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATg                                                                                                                                    <  1:863157/150‑1 (MQ=255)
                     tCCTGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGAc                                                                                                                                  <  2:1166279/150‑1 (MQ=255)
                     tCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGAc                                                                                                                                  <  2:933168/150‑1 (MQ=255)
                     tCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGAc                                                                                                                                  <  2:265493/150‑1 (MQ=255)
                      ccGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATg                                                                                                                                    <  1:920310/147‑1 (MQ=255)
                      ccGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATg                                                                                                                                    >  2:920310/1‑147 (MQ=255)
                         gCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAAcc                                                                                                                              >  1:761447/1‑150 (MQ=255)
                         gCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAAcc                                                                                                                              >  2:1269904/1‑150 (MQ=255)
                         gCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAAcc                                                                                                                              >  2:739320/1‑150 (MQ=255)
                          cGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCg                                                                                                                             <  1:691637/150‑1 (MQ=255)
                              gTCGAGCGTTGGGTGGATTGTGAACATCACGCAAGACAGATAAAGACAAAATAAATAGCCGAGCACAACCTTCTTGTAAAGAAGTTTAGCGGCTTCGGTCATGACGAGCTTTAACCAGATGACCTTGAAATCTTTCATGACAACCCGCAc                                                                                                                         >  2:789668/1‑150 (MQ=255)
                                cGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGc                                                                                                                       <  1:1169763/150‑1 (MQ=255)
                                 gAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGcc                                                                                                                      <  1:990123/150‑1 (MQ=255)
                                      ttGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTcc                                                                                                                 >  1:500875/1‑150 (MQ=255)
                                       tGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCa                                                                                                                <  2:320893/150‑1 (MQ=255)
                                         ggTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACg                                                                                                              <  1:545839/150‑1 (MQ=255)
                                         ggTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACg                                                                                                              <  1:788423/150‑1 (MQ=255)
                                         ggTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACg                                                                                                              >  2:877032/1‑150 (MQ=255)
                                          gTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGc                                                                                                             <  1:540680/150‑1 (MQ=255)
                                           tGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCa                                                                                                            <  1:1348040/150‑1 (MQ=255)
                                            ggATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACg                                                                                                              >  2:223588/1‑147 (MQ=255)
                                            ggATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACg                                                                                                              <  1:223588/147‑1 (MQ=255)
                                              aTTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCCTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGAc                                                                                                         <  1:1012486/150‑1 (MQ=255)
                                                 gtgGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTgg                                                                                                      >  1:614493/1‑150 (MQ=255)
                                                 gtgGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTgg                                                                                                      >  1:932222/1‑150 (MQ=255)
                                                   ggAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTa                                                                                                    <  1:877032/150‑1 (MQ=255)
                                                    gagaTCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTaa                                                                                                   >  2:714287/1‑150 (MQ=255)
                                                     agaTCACGCCAGCCAGTTTCAGACCATCGATCCTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAg                                                                                                  <  2:358511/150‑1 (MQ=255)
                                                      gaTCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGc                                                                                                 >  2:501430/1‑150 (MQ=255)
                                                      gaTCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGGTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGc                                                                                                 <  2:116354/150‑1 (MQ=255)
                                                           cGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATgc                                                                                            <  1:225534/150‑1 (MQ=255)
                                                             ccagccagTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATgccg                                                                                          <  1:834129/150‑1 (MQ=255)
                                                              cagccagTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAAcc                                                                                                                              >  2:960739/1‑113 (MQ=255)
                                                               agccagTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATgccgcc                                                                                        >  1:1372831/1‑150 (MQ=255)
                                                                gccagTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATgccgccg                                                                                       >  2:240982/1‑150 (MQ=255)
                                                                  cagTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATgccgccgcc                                                                                     <  1:354014/150‑1 (MQ=255)
                                                                   agTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGccc                                                                                    <  2:1296792/150‑1 (MQ=255)
                                                                    gTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCa                                                                                   <  1:980016/150‑1 (MQ=255)
                                                                         aGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATg                                                                              >  1:337480/1‑150 (MQ=255)
                                                                            ccATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCt                                                                           >  2:427369/1‑150 (MQ=255)
                                                                            ccATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCt                                                                           <  2:741496/150‑1 (MQ=255)
                                                                            ccATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCt                                                                           >  1:676080/1‑150 (MQ=255)
                                                                             cATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTg                                                                          <  1:1066354/150‑1 (MQ=255)
                                                                             cATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTg                                                                          >  2:431272/1‑150 (MQ=255)
                                                                              aTCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGc                                                                         >  2:1297045/1‑150 (MQ=255)
                                                                              aTCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCa                                                                        >  2:1172106/1‑150 (MQ=255)
                                                                               tCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCa                                                                        <  1:1048463/150‑1 (MQ=255)
                                                                               tCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCa                                                                        >  1:786876/1‑150 (MQ=255)
                                                                                       gCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCAc                                                                >  2:459986/1‑150 (MQ=255)
                                                                                        ccGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCAcc                                                               >  1:629215/1‑150 (MQ=255)
                                                                                         cgcgCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACcg                                                              >  2:69881/1‑150 (MQ=255)
                                                                                         cgcgCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACcg                                                              <  1:605925/150‑1 (MQ=255)
                                                                                            gcTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCgg                                                           <  2:75226/150‑1 (MQ=255)
                                                                                             cTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGc                                                          >  2:105811/1‑150 (MQ=255)
                                                                                             cTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGc                                                          >  1:893678/1‑150 (MQ=255)
                                                                                              tCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGcc                                                         <  1:711547/150‑1 (MQ=255)
                                                                                                    tgctTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTc                                                   >  1:729282/1‑150 (MQ=255)
                                                                                                     gctTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCg                                                  >  2:634140/1‑150 (MQ=255)
                                                                                                     gctTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCg                                                  >  1:884317/1‑150 (MQ=255)
                                                                                                      ctTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCgg                                                 >  2:108229/1‑150 (MQ=255)
                                                                                                      ctTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCgg                                                 <  1:1269904/150‑1 (MQ=255)
                                                                                                            cagcagTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTc                                           >  2:1051223/1‑150 (MQ=255)
                                                                                                            cagcagTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTc                                           >  2:1415701/1‑150 (MQ=255)
                                                                                                             agcagTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCa                                          <  1:1417612/150‑1 (MQ=255)
                                                                                                             agcagTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCa                                          <  2:337480/150‑1 (MQ=255)
                                                                                                                agTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACg                                       <  2:786876/150‑1 (MQ=255)
                                                                                                                       cgGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTgg                                <  1:229437/150‑1 (MQ=255)
                                                                                                                         gCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTg                              <  1:105811/150‑1 (MQ=255)
                                                                                                                         gCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTg                              <  2:832654/150‑1 (MQ=255)
                                                                                                                            tCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGt                           >  2:295864/1‑150 (MQ=255)
                                                                                                                                tCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTttct                       >  1:26228/1‑150 (MQ=255)
                                                                                                                                   tGCCGATGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCg                    <  1:1272833/150‑1 (MQ=255)
                                                                                                                                   tGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCg                    >  2:648379/1‑150 (MQ=255)
                                                                                                                                    gCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGcc                                                                            <  1:550297/93‑1 (MQ=255)
                                                                                                                                    gCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGcc                                                                            >  2:550297/1‑93 (MQ=255)
                                                                                                                                    gCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACACAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGc                   >  1:372646/1‑150 (MQ=255)
                                                                                                                                     ccGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCt                  >  2:514902/1‑150 (MQ=255)
                                                                                                                                       gAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCttt                <  2:1057876/150‑1 (MQ=255)
                                                                                                                                        aGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTg               <  2:426531/150‑1 (MQ=255)
                                                                                                                                           gTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCt            >  1:58012/1‑150 (MQ=255)
                                                                                                                                           atcaaCGACAAGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCt            >  1:437750/4‑150 (MQ=255)
                                                                                                                                             caaCGACAAGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCttt          <  2:437750/149‑1 (MQ=255)
                                                                                                                                              aaCGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCtttt         <  2:394730/150‑1 (MQ=255)
                                                                                                                                               aCGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTa        <  2:1372831/150‑1 (MQ=255)
                                                                                                                                                cGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTAc       <  1:69881/150‑1 (MQ=255)
                                                                                                                                                 gACTTGGCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTCTCTTCGCTTTGCCTTTTACa      <  1:1112943/150‑1 (MQ=255)
                                                                                                                                                     tGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTACATAtt  <  1:427369/150‑1 (MQ=255)
                                                                                                                                                     tGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTACATAtt  >  2:1229714/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACAATATCCACGCGGTCAAATCCGGCGTAGTCGAGCGTTGGGTGGATTGTGGAGATCACGCCAGCCAGTTTCAGACCATCGATCTTGCCGCGCTCAAGCTGCTTGTTCAGCAGTTTCGCGGCTTCGGTCATGCCGAGGGTTAACGACTTGTCGTTGATATCTTTCATGACAACCGGCACGCCTTTCCACGCAGACTGGTAAGCGATGCCGCCGCCCATAATGCCTGCACCCAGCACCGCGGCCTGTTTCGGGGTTTCAACGTCTTTGGTGAGTTTCTTCGCTTTGCCTTTTACATATT  >  NC_000913/4029784‑4030082

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