Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,128 T→A 5.4% M678L (ATG→TTG)  fadJ ← enoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,1280TA 5.4% 249.1 / 5.1 74M678L (ATG→TTG) fadJenoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base A (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

GCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCA  >  NC_000913/2456980‑2457276
                                                                                                                                                    |                                                                                                                                                    
gCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCat                                                                                                                                                      <  1:269880/149‑1 (MQ=255)
                 cATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGc                                                                                                                                     <  4:636218/149‑1 (MQ=255)
                  atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                                                                    >  4:441902/1‑149 (MQ=255)
                    ataGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAACCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCaa                                                                                                                                              >  1:127754/1‑137 (MQ=255)
                    ataGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCaa                                                                                                                                              <  2:127754/137‑1 (MQ=255)
                     taGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcgc                                                                                                                                 >  4:37194/1‑149 (MQ=255)
                         tttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCACGCGCCg                                                                                                                             >  1:647950/1‑149 (MQ=255)
                          ttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCga                                                                                                                            <  3:854689/149‑1 (MQ=255)
                            gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgaga                                                                                                                          >  4:85024/1‑149 (MQ=255)
                                 tcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                                                                     <  2:651467/149‑1 (MQ=255)
                                   ttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                   >  2:246259/1‑149 (MQ=255)
                                     aTTGCAGGTCAGTTGCTGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGGCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTTGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGata                                                                                                                 <  2:807010/149‑1 (MQ=255)
                                       tGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGAtata                                                                                                               <  1:439915/149‑1 (MQ=255)
                                          aGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCg                                                                                                            >  4:232724/1‑149 (MQ=255)
                                            gTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGa                                                                                                          >  3:179484/1‑149 (MQ=255)
                                               agttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACg                                                                                                       >  3:376861/1‑149 (MQ=255)
                                                     agttgTTTTCCAAAAACTTACCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAc                                                                                                 >  2:143801/1‑149 (MQ=255)
                                                      gttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAcc                                                                                                >  4:219272/1‑149 (MQ=255)
                                                       ttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCg                                                                                               >  3:730927/1‑149 (MQ=255)
                                                       ttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCg                                                                                               >  4:505047/1‑149 (MQ=255)
                                                         gttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCTACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                             >  2:328947/1‑149 (MQ=255)
                                                            ttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaaa                                                                                          <  3:441902/149‑1 (MQ=255)
                                                              ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                                                                                  >  3:151769/1‑123 (MQ=255)
                                                              ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                                                                                  <  4:151769/123‑1 (MQ=255)
                                                               cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATgg                                                                                       <  4:568838/149‑1 (MQ=255)
                                                                agaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGTAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      >  3:819137/3‑149 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGTCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAc                                                                                                                                             <  2:509978/94‑1 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAc                                                                                                                                             >  1:509978/1‑94 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTc                                                                                                    <  4:807303/135‑1 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTc                                                                                                    >  3:807303/1‑135 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      <  3:185349/149‑1 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      >  2:629846/1‑149 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      >  2:71811/1‑149 (MQ=255)
                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGAGGAACGGTCCACCGAGAAATGGc                                                                                      <  2:667209/149‑1 (MQ=255)
                                                                      tttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaaaa                                                                                <  1:246259/149‑1 (MQ=255)
                                                                       ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAAc                                                                               <  1:828398/149‑1 (MQ=255)
                                                                           ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                           <  2:762598/149‑1 (MQ=255)
                                                                           ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                           <  3:793027/149‑1 (MQ=255)
                                                                           ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                           <  3:505047/149‑1 (MQ=255)
                                                                            cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                                                                          <  4:179484/149‑1 (MQ=255)
                                                                            cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                                                                          >  1:535092/1‑149 (MQ=255)
                                                                                 gcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGc                                                                                                             >  1:284328/1‑109 (MQ=255)
                                                                                 gcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGc                                                                                                             <  2:284328/109‑1 (MQ=255)
                                                                                 gcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaaa                                                                     <  1:337104/149‑1 (MQ=255)
                                                                                     ccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAcc                                                                 >  4:431980/1‑149 (MQ=255)
                                                                                       aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcgc                                                               <  4:730927/149‑1 (MQ=255)
                                                                                               cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAatat                                                       >  4:797035/1‑149 (MQ=255)
                                                                                                aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATc                                                      >  1:502457/1‑149 (MQ=255)
                                                                                                        gCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcac                                              <  2:535092/149‑1 (MQ=255)
                                                                                                               gTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      <  4:390883/102‑1 (MQ=255)
                                                                                                               gTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                      >  3:390883/1‑102 (MQ=255)
                                                                                                               gTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCt                                       >  2:818521/1‑149 (MQ=255)
                                                                                                                 aaaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTAc                                     <  3:232724/149‑1 (MQ=255)
                                                                                                                   aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                   <  1:315958/149‑1 (MQ=255)
                                                                                                                   aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                   <  3:30199/149‑1 (MQ=255)
                                                                                                                      gggAACCATACTGCGTGTCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGACAATCGATATAGCGGAACGGt                                                                                                     <  4:654128/80‑1 (MQ=255)
                                                                                                                      gggAACCATACTGCGTGGCGAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATa                                <  3:85024/149‑1 (MQ=255)
                                                                                                                      gggAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                                                                     >  3:654128/1‑80 (MQ=255)
                                                                                                                        gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAc                              <  2:775981/149‑1 (MQ=255)
                                                                                                                        gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAc                              <  3:652897/149‑1 (MQ=255)
                                                                                                                           ccATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTg                           >  3:404895/1‑149 (MQ=255)
                                                                                                                             aTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCt                         <  1:731330/149‑1 (MQ=255)
                                                                                                                              tACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtc                        <  3:727013/149‑1 (MQ=255)
                                                                                                                                 tGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtcatc                     <  1:13274/149‑1 (MQ=255)
                                                                                                                                 tGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtcatc                     >  2:13274/1‑149 (MQ=255)
                                                                                                                                  gCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtcatca                    <  1:818521/149‑1 (MQ=255)
                                                                                                                                    gggcgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAac                  >  4:646916/6‑149 (MQ=255)
                                                                                                                                           gtggTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGAt                                 <  4:139162/124‑1 (MQ=255)
                                                                                                                                           gtggTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGAt                                 >  3:139162/4‑127 (MQ=255)
                                                                                                                                            tCGTTGCATTATTGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTa          <  2:161915/149‑1 (MQ=255)
                                                                                                                                             cGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTAc         <  3:4951/149‑1 (MQ=255)
                                                                                                                                             cGTTGCATTATTGCAACCACTTCGCCCGCGCCAAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTAc         <  2:647950/149‑1 (MQ=255)
                                                                                                                                               ttGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCaa                                                          >  2:112674/1‑98 (MQ=255)
                                                                                                                                               ttGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCaa                                                          <  1:112674/98‑1 (MQ=255)
                                                                                                                                                 gCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCt     <  3:646916/149‑1 (MQ=255)
                                                                                                                                                   attattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttc   >  2:33241/1‑149 (MQ=255)
                                                                                                                                                   attattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttc   <  3:447302/149‑1 (MQ=255)
                                                                                                                                                    ttattGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttca  <  1:145354/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCA  >  NC_000913/2456980‑2457276

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