Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,138 C→G 6.2% V674V (GTG→GTC)  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,1380CG 6.2% 269.9 / 6.2 78V674V (GTG→GTC) fadJenoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base G (2/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCA  >  NC_000913/2456997‑2457286
                                                                                                                                             |                                                                                                                                                    
cATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGc                                                                                                                                               <  4:636218/149‑1 (MQ=255)
 atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                                                                              >  4:441902/1‑149 (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)
                                               aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTc                                                                                                              >  3:807303/1‑135 (MQ=255)
                                               aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTc                                                                                                              <  4:807303/135‑1 (MQ=255)
                                               aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                                <  3:185349/149‑1 (MQ=255)
                                               aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                                >  2:71811/1‑149 (MQ=255)
                                               aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                                >  2:629846/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                                                                                     <  3:505047/149‑1 (MQ=255)
                                                          ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                                     <  2:762598/149‑1 (MQ=255)
                                                          ccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                                     <  3:793027/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                                                                                                >  3:390883/1‑102 (MQ=255)
                                                                                              gTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGc                                                                                                <  4:390883/102‑1 (MQ=255)
                                                                                              gTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCt                                                 >  2:818521/1‑149 (MQ=255)
                                                                                                aaaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTAc                                               <  3:232724/149‑1 (MQ=255)
                                                                                                  aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                             <  3:30199/149‑1 (MQ=255)
                                                                                                  aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                             <  1:315958/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                                        <  3:652897/149‑1 (MQ=255)
                                                                                                       gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAc                                        <  2:775981/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                                                                    <  1:112674/98‑1 (MQ=255)
                                                                                                                              ttGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCaa                                                                    >  2:112674/1‑98 (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)
                                                                                                                                     attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACTGATAACCTGCTCATCAACACAACGTACTGCttcatt          <  4:819137/149‑1 (MQ=255)
                                                                                                                                        gCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAg       >  4:426010/1‑149 (MQ=255)
                                                                                                                                         cAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGc      >  3:388246/1‑149 (MQ=255)
                                                                                                                                          aaCCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCa     >  2:462190/1‑149 (MQ=255)
                                                                                                                                           aCCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCAt    >  4:325443/1‑149 (MQ=255)
                                                                                                                                             cACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCa  >  2:649694/1‑149 (MQ=255)
                                                                                                                                             cACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCa  >  2:765246/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCA  >  NC_000913/2456997‑2457286

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