Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 648,555 A→G R339R (CGT→CGC citF ← citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913648,5550AG100.0% 230.5 / NA 65R339R (CGT→CGCcitFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/32);  total (33/32)

CGCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGT  >  NC_000913/648407‑648700
                                                                                                                                                    |                                                                                                                                                 
cGCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCg                                                                                                                                                 <  2:312181/151‑1 (MQ=255)
 gCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCgg                                                                                                                                                >  2:576269/1‑151 (MQ=255)
       agagCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGc                                                                                                                                          <  2:54558/151‑1 (MQ=255)
       agagCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGc                                                                                                                                          >  2:264968/1‑151 (MQ=255)
           ccGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGAtt                                                                                                                                      <  1:464179/151‑1 (MQ=255)
             ggTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTgc                                                                                                                                    <  1:479413/151‑1 (MQ=255)
               tCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTgcgc                                                                                                                                  >  1:84627/1‑151 (MQ=255)
                cAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCgca                                                                                                                                 <  2:354177/151‑1 (MQ=255)
                 aGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCgcag                                                                                                                                <  2:548847/151‑1 (MQ=255)
                     cGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGCTGCGGGCCAGCGATTGCGCAGCCTg                                                                                                                            >  2:175647/1‑151 (MQ=255)
                         aaCGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCt                                                                                                                         <  2:732105/150‑1 (MQ=255)
                         aaCGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTg                                                                                                                        >  2:634971/1‑151 (MQ=255)
                                   tGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGctct                                                                                                              >  1:151748/1‑151 (MQ=255)
                                   tGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGctc                                                                                                               >  2:309099/1‑150 (MQ=255)
                                    gggTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTg                                                                                                             >  2:524682/1‑151 (MQ=255)
                                          cAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACAt                                                                                                        <  1:125318/150‑1 (MQ=255)
                                             tttCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATcc                                                                                                      >  1:414484/1‑149 (MQ=255)
                                             tttCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATcc                                                                                                      >  1:755296/1‑149 (MQ=255)
                                                ccAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCcagcag                                                                                                 <  2:719282/151‑1 (MQ=255)
                                                   gcgcgcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGttt                                                                                              <  1:696875/151‑1 (MQ=255)
                                                       gcTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCgg                                                                                          >  1:218847/1‑151 (MQ=255)
                                                         tCAGTACCACCACGTCGCGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGa                                                                                         >  2:423309/1‑150 (MQ=255)
                                                         tCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGAt                                                                                        >  1:311332/1‑151 (MQ=255)
                                                          cAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGAt                                                                                        <  1:103945/150‑1 (MQ=255)
                                                          cAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGAt                                                                                        <  1:446024/150‑1 (MQ=255)
                                                          cAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       <  2:12164/151‑1 (MQ=255)
                                                          cAGTACCACCACGGCGAGACGAGCAACCGAGGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       <  2:173339/151‑1 (MQ=255)
                                                             tACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATc                                                                                       >  2:547591/1‑148 (MQ=255)
                                                              accaccACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGAc                                                                                   <  2:554676/151‑1 (MQ=255)
                                                               ccaccaCGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGAcc                                                                                  <  2:189946/151‑1 (MQ=255)
                                                                caccacGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCt                                                                                 <  2:522633/151‑1 (MQ=255)
                                                                    acGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCttttt                                                                             <  2:673160/151‑1 (MQ=255)
                                                                     cGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTc                                                                            <  2:702041/151‑1 (MQ=255)
                                                                         tagaCGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGCTGGCGCTGATTTCGATGTGATTGGGGTGGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTGTGCGGATCAGACATTTTTCGt                                                                          >  2:544782/2‑149 (MQ=255)
                                                                            aCGATCAACCGATGCGCCTTTCGATCCCCAGTTAGCGTGCTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGCCCTTTTTCGTGCAg                                                                      <  2:139508/150‑1 (MQ=255)
                                                                               aTCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCa                                                                  <  1:692188/151‑1 (MQ=255)
                                                                                  aaCCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAAcc                                                               >  1:348736/1‑151 (MQ=255)
                                                                                     cGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                            >  2:638752/1‑151 (MQ=255)
                                                                                         gcgcCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtc                                                         <  2:726988/150‑1 (MQ=255)
                                                                                          cgcCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                       >  1:731266/1‑151 (MQ=255)
                                                                                          cgcCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                       >  2:461348/1‑151 (MQ=255)
                                                                                           gcCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgcg                                                      >  2:197125/1‑151 (MQ=255)
                                                                                             cTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGt                                                    >  2:628889/1‑151 (MQ=255)
                                                                                                  gAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAc                                               >  2:68330/1‑151 (MQ=255)
                                                                                                   aaCCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAcc                                              <  1:747002/151‑1 (MQ=255)
                                                                                                   aaCCCCAGTTAGCGTACTGGTTGGCGCTGATTTAGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAcc                                              >  2:333739/1‑151 (MQ=255)
                                                                                                    aCCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCg                                             <  1:614544/151‑1 (MQ=255)
                                                                                                     ccccAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGc                                            <  2:600555/151‑1 (MQ=255)
                                                                                                          gTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAGCCATCGTCGCGGTAATACCGCCAAgg                                       <  2:359158/151‑1 (MQ=255)
                                                                                                            tAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGc                                     >  2:68364/1‑151 (MQ=255)
                                                                                                             aGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCg                                    >  1:594271/1‑151 (MQ=255)
                                                                                                                 tACTGGTTGGGGCTGATTTCGTTGTGATTGGGGTTGCGGGCGAGCGATTGCGCAGCATGGCTGTTAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACGTTTTTCGTGCAGGTCAACCATCGTCGCGGTACTTCCGCCAAGGGCGa                                   <  2:407240/148‑1 (MQ=255)
                                                                                                                    tGGTTGGCGCTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCgg                             >  1:235235/1‑151 (MQ=255)
                                                                                                                          gcgcTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGaa                       <  1:189859/151‑1 (MQ=255)
                                                                                                                            gcTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAat                      >  1:185966/1‑150 (MQ=255)
                                                                                                                            gcTGATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAat                      >  1:417677/1‑150 (MQ=255)
                                                                                                                               gATTTCGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcg                  >  1:420390/1‑151 (MQ=255)
                                                                                                                                    cGATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCt             <  1:305720/151‑1 (MQ=255)
                                                                                                                                     gATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCt             >  1:610040/1‑150 (MQ=255)
                                                                                                                                     gATGTGATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACAACCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGCCACCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCTAATAGCGCGGc              <  2:219643/149‑1 (MQ=255)
                                                                                                                                       tgtgAGTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGTTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCGAGGGCGAAGTCGGCGCGAATATCGCGGCTAc           <  2:54847/150‑1 (MQ=255)
                                                                                                                                          gATTGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCa        <  1:461369/150‑1 (MQ=255)
                                                                                                                                            ttGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAttt     >  2:232523/1‑151 (MQ=255)
                                                                                                                                            ttGGGGTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAttt     >  1:520381/1‑151 (MQ=255)
                                                                                                                                               gggTTGCGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGt  >  1:482229/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGT  >  NC_000913/648407‑648700

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