Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,028,758 G→A 5.5% V5V (GTC→GTT fadA ← 3‑ketoacyl‑CoA thiolase (thiolase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,028,7580GA 5.5% 233.7 / 3.2 73V5V (GTC→GTTfadA3‑ketoacyl‑CoA thiolase (thiolase I)
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (2/2);  total (35/38)
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.47e-02

ATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATAT  >  NC_000913/4028615‑4028905
                                                                                                                                               |                                                                                                                                                   
atatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAAc                                                                                                                                                <  2:573467/149‑1 (MQ=255)
atatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAAc                                                                                                                                                <  3:27194/149‑1 (MQ=255)
  atCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCt                                                                                                                                              >  1:240489/1‑149 (MQ=255)
     gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGtt                                                                                                                                           >  3:329262/1‑149 (MQ=255)
      tcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTc                                                                                                                                          <  4:99697/149‑1 (MQ=255)
              cgccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGa                                                                                                                                  >  2:103748/1‑149 (MQ=255)
                ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACt                                                                                                                                >  2:477343/1‑149 (MQ=255)
                ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACt                                                                                                                                >  1:443809/1‑149 (MQ=255)
                ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACt                                                                                                                                >  1:382668/1‑149 (MQ=255)
                 cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTc                                                                                                                               <  4:11553/149‑1 (MQ=255)
                 cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTc                                                                                                                               <  3:56993/149‑1 (MQ=255)
                   cTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCt                                                                                                                             <  2:265479/149‑1 (MQ=255)
                    ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATtgt                                                                                                                                    >  3:440030/1‑141 (MQ=255)
                    ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATtgt                                                                                                                                    <  4:440030/141‑1 (MQ=255)
                      ccAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTaa                                                                                                                          <  4:544989/149‑1 (MQ=255)
                       cAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAAAGGCACATCGGGCTGCGAATTGTCTCGACATTGACAACCTGTTCCATTGTGACACCTTAAg                                                                                                                         <  4:233516/149‑1 (MQ=255)
                        aGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGc                                                                                                                        >  4:133453/1‑149 (MQ=255)
                          cgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCg                                                                                                                      >  2:570705/1‑149 (MQ=255)
                           gcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGGTCTGCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGt                                                                                                                     >  4:90942/1‑149 (MQ=255)
                            ccGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGtt                                                                                                                    >  3:416397/1‑149 (MQ=255)
                              gggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGtttt                                                                                                                  >  1:205238/1‑149 (MQ=255)
                               ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTc                                                                                                                 <  4:500268/149‑1 (MQ=255)
                                 ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAg                                                                                                               <  1:372007/149‑1 (MQ=255)
                                 ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGAGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAg                                                                                                               <  4:477511/149‑1 (MQ=255)
                                    cgcgCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTc                                                                                                            <  4:329262/149‑1 (MQ=255)
                                    cgcgCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTc                                                                                                            <  1:69255/149‑1 (MQ=255)
                                       gcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGcc                                                                                                         >  4:422434/1‑149 (MQ=255)
                                               gCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGAc                                                                                                 >  4:217667/1‑149 (MQ=255)
                                                 tACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCATGTCGCCAACCGGACgg                                                                                               >  1:324695/1‑149 (MQ=255)
                                                 tACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATAACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCa                                                                                                                >  1:537180/1‑132 (MQ=255)
                                                 tACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATAACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCa                                                                                                                <  2:537180/132‑1 (MQ=255)
                                                 tACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATAACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCa                                                                                                                >  3:114453/1‑132 (MQ=255)
                                                 tACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATAACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCa                                                                                                                <  4:114453/132‑1 (MQ=255)
                                                      aTTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGggctgg                                                                                          >  2:452836/1‑149 (MQ=255)
                                                       ttAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGggctggc                                                                                         >  4:14891/1‑149 (MQ=255)
                                                                cGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCgga                                                                                <  2:497595/149‑1 (MQ=255)
                                                                   agagaTCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACg                                                                                                <  4:138773/130‑1 (MQ=255)
                                                                   agagaTCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACg                                                                                                >  3:138773/1‑130 (MQ=255)
                                                                       aTCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTc                                                                                                            >  1:515317/1‑114 (MQ=255)
                                                                       aTCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTc                                                                                                            <  2:515317/114‑1 (MQ=255)
                                                                       aTCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGt                                                                         >  4:55590/1‑149 (MQ=255)
                                                                        tcttctGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTa                                                                        <  1:328364/149‑1 (MQ=255)
                                                                         cttctGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTAc                                                                       >  1:281651/1‑149 (MQ=255)
                                                                              gcacgcacGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTc                                                                  <  3:133453/149‑1 (MQ=255)
                                                                                  gcacgTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTa                                                              >  4:552258/1‑149 (MQ=255)
                                                                                          cGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCCTTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACgcgc                                                      <  4:197004/149‑1 (MQ=255)
                                                                                           gAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCt                                                     >  3:65332/1‑149 (MQ=255)
                                                                                             aaaGCACCGCCATTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCttt                                                   <  2:205238/149‑1 (MQ=255)
                                                                                                     gCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCa                                           <  3:332161/149‑1 (MQ=255)
                                                                                                     gCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCa                                           <  3:162383/149‑1 (MQ=255)
                                                                                                         ttCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATa                                                                           >  3:555564/1‑113 (MQ=255)
                                                                                                         ttCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATa                                                                           <  4:555564/113‑1 (MQ=255)
                                                                                                            gAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTc                                    >  2:361407/1‑149 (MQ=255)
                                                                                                             aaCGGCCCATCGGGGTGCGAATTGCATCGACGATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTcc                                   >  4:286354/1‑149 (MQ=255)
                                                                                                             aaCGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCa                                           <  3:527452/141‑1 (MQ=255)
                                                                                                             aaCGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCa                                           >  4:527452/1‑141 (MQ=255)
                                                                                                             aaCGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTcc                                   <  3:55590/149‑1 (MQ=255)
                                                                                                               cGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCgg                                 <  3:90942/149‑1 (MQ=255)
                                                                                                                  cccATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCAc                              <  3:297906/149‑1 (MQ=255)
                                                                                                                   ccATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCAGGCACt                             >  2:42672/1‑149 (MQ=255)
                                                                                                                       cGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCAt                         <  1:570705/149‑1 (MQ=255)
                                                                                                                           gTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAg                     >  2:359921/1‑149 (MQ=255)
                                                                                                                                    gCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAg            <  4:100790/149‑1 (MQ=255)
                                                                                                                                    gCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAg            <  3:217667/149‑1 (MQ=255)
                                                                                                                                      aTCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGGTAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGt          >  4:129974/1‑149 (MQ=255)
                                                                                                                                        cGACAATGACAACCTGTTCCATTGTGACTCTTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGc        <  1:359921/149‑1 (MQ=255)
                                                                                                                                          aCAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTg      >  1:164612/1‑149 (MQ=255)
                                                                                                                                            aaTGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGat    >  4:247558/1‑149 (MQ=255)
                                                                                                                                             aTGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGata   <  4:79488/149‑1 (MQ=255)
                                                                                                                                              tGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGatat  <  3:155953/149‑1 (MQ=255)
                                                                                                                                              tGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGatat  <  3:14891/149‑1 (MQ=255)
                                                                                                                                              tGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGata   >  4:120026/1‑148 (MQ=255)
                                                                                                                                              tGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGata   <  3:120026/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCATCGACAATGACAACCTGTTCCATTGTGACTCCTTAAGCCGTTTTCAGGTCGCCAACCGGACGGGCTGGCTCAACCGGAGGATAGTACGGTTCGTTATGACGCGCTTTATTACGCAGACCTTCCGGCACTTCATACAGCGGGCCGAGGTGCTGATAT  >  NC_000913/4028615‑4028905

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