Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,128 T→A 8.3% M678L (ATG→TTG)  fadJ ← 3‑hydroxyacyl‑CoA dehydrogenase FadJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,1280TA 8.3% 182.9 / 4.4 59M678L (ATG→TTG) fadJ3‑hydroxyacyl‑CoA dehydrogenase FadJ
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base A (2/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCA  >  NC_000913/2456984‑2457276
                                                                                                                                                 |                                                                                                                                                    
ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTg                                                                                                                                                 <  1:436795/150‑1 (MQ=255)
ttCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGAGACCATCTCGACCAAAAGCTCGAAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTg                                                                                                                                                 >  1:1670531/1‑150 (MQ=255)
  cAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCa                                                                                                                                               <  1:90883/150‑1 (MQ=255)
          aTTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCAATAATGCAACGACttg                                                                                                                                       <  1:994707/150‑2 (MQ=255)
            tCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGc                                                                                                                                     >  1:1061884/1‑150 (MQ=255)
             cATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                                                                    >  1:764609/1‑150 (MQ=255)
                 taGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgccacgc                                                                                                                                 <  1:419725/150‑8 (MQ=38)
                      ttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                                                                           <  2:830787/150‑1 (MQ=255)
                        gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGCCCAACCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagag                                                                                                                         >  2:1261264/1‑150 (MQ=255)
                        gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagag                                                                                                                         <  1:543467/150‑1 (MQ=255)
                          ccTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGaa                                                                                                                       >  2:141549/1‑150 (MQ=255)
                           cttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAAt                                                                                                                      >  1:1778381/1‑150 (MQ=255)
                            ttcttATTGCAGGTCAGTTGCACTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGAGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                                                                     <  1:83812/150‑1 (MQ=255)
                             tcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                                                                    <  2:775152/150‑1 (MQ=255)
                                   tGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAg                                                                                                              >  1:1003578/1‑150 (MQ=255)
                                    gCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGACAGAATCGATATAGc                                                                                                             <  1:1748596/150‑1 (MQ=255)
                                      aGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCgg                                                                                                           >  1:220106/1‑150 (MQ=255)
                                           agttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACgg                                                                                                      >  1:843249/1‑150 (MQ=255)
                                              tgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTcc                                                                                                   >  2:1595571/1‑150 (MQ=255)
                                               gcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCa                                                                                                  <  2:1134420/150‑1 (MQ=255)
                                                    tgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                             <  2:1743929/150‑1 (MQ=255)
                                                       tttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaaa                                                                                          >  2:599442/1‑150 (MQ=255)
                                                        ttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAAt                                                                                         <  1:687239/150‑1 (MQ=255)
                                                         tCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATg                                                                                        >  2:1222366/1‑150 (MQ=255)
                                                           cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      <  1:141549/150‑1 (MQ=255)
                                                            aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCg                                                                                     <  2:1176480/150‑1 (MQ=255)
                                                             aaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                                                                    <  1:1974898/150‑1 (MQ=255)
                                                             aaaaCTTTCCCAACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                                                                    <  2:277764/150‑1 (MQ=255)
                                                                 cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaaaa                                                                                <  1:460537/150‑1 (MQ=255)
                                                                 cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaaaa                                                                                <  2:1061884/150‑1 (MQ=255)
                                                                   ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAAcc                                                                              <  2:86218/150‑1 (MQ=255)
                                                                      cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                           <  2:1563585/150‑1 (MQ=255)
                                                                             gcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAAt                                                                    <  1:312626/150‑1 (MQ=255)
                                                                                   aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcgcg                                                              >  1:747346/1‑150 (MQ=255)
                                                                                   aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGA‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCTATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcgcg                                                              <  2:1670531/150‑1 (MQ=255)
                                                                                         aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAata                                                        <  1:93103/150‑1 (MQ=255)
                                                                                            aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc                                                     <  1:599442/150‑1 (MQ=255)
                                                                                                gCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGt                                                 <  2:594258/150‑1 (MQ=255)
                                                                                                   cGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcac                                              <  1:1547807/150‑1 (MQ=255)
                                                                                                       aGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgcac                                          >  1:1421172/1‑150 (MQ=255)
                                                                                                       aGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgcac                                          >  1:1589646/1‑150 (MQ=255)
                                                                                                           gTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTa                                      <  1:743485/150‑1 (MQ=255)
                                                                                                             aaaaCGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACg                                    >  2:1052836/1‑150 (MQ=255)
                                                                                                               aaCGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCTGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                  <  2:778283/150‑1 (MQ=255)
                                                                                                                    gAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                             <  2:1421172/150‑1 (MQ=255)
                                                                                                                    gAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                             <  2:1589646/150‑1 (MQ=255)
                                                                                                                     aaCCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCt                            >  2:917908/1‑150 (MQ=255)
                                                                                                                            cTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtcatc                     >  1:1359325/1‑150 (MQ=255)
                                                                                                                                  ggC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAAcaca               >  2:1942727/1‑150 (MQ=255)
                                                                                                                                   gC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACaa              >  1:1694411/1‑150 (MQ=255)
                                                                                                                                    cgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTca                                               <  1:842004/115‑1 (MQ=255)
                                                                                                                                    cgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTca                                               >  2:842004/3‑117 (MQ=255)
                                                                                                                                    c‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGACGAGAGAATCGATATAGCGGAACGGTCAACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAAc             >  2:228984/1‑150 (MQ=255)
                                                                                                                                        gtggTTGCAATAATGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTa          >  1:1217766/4‑150 (MQ=255)
                                                                                                                                          cGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACt        <  1:1595571/150‑1 (MQ=255)
                                                                                                                                           gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTg       >  2:1480859/1‑150 (MQ=255)
                                                                                                                                              gCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCtt    >  1:12992/1‑150 (MQ=255)
                                                                                                                                              gCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCtt    >  2:441710/1‑150 (MQ=255)
                                                                                                                                              gCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCtt    >  2:429818/1‑150 (MQ=255)
                                                                                                                                                attattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttca  <  1:580526/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCA  >  NC_000913/2456984‑2457276

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: ‑