Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,529 G→A 6.2% R85R (CGC→CGT citD ← citrate lyase, acyl carrier (gamma) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,5290GA 6.2% 261.4 / 6.0 82R85R (CGC→CGTcitDcitrate lyase, acyl carrier (gamma) subunit
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base A (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGC  >  NC_000913/650382‑650677
                                                                                                                                                   |                                                                                                                                                    
cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcgc                                                                                                                                                     >  1:319650/1‑149 (MQ=255)
           aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                                                                                                                                          >  3:707917/1‑149 (MQ=255)
              aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCt                                                                                                                                       <  3:790421/149‑1 (MQ=255)
                ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                                                                     >  4:400066/1‑149 (MQ=255)
                ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                                                                     >  1:417767/1‑149 (MQ=255)
                  aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                                                                   <  3:123836/149‑1 (MQ=255)
                     tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                                                                                                <  1:390993/149‑1 (MQ=255)
                       cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCTGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                                                                                              <  1:167139/149‑1 (MQ=255)
                          aCCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCa                                                                                                                           >  1:45848/1‑149 (MQ=255)
                          aCCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCa                                                                                                                           >  2:813643/1‑149 (MQ=255)
                             aTCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCCTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCCGGGCTTCCCGTCGTGCCCGt                                                                                                                        >  2:488902/1‑149 (MQ=255)
                             aTCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGt                                                                                                                        <  3:241282/149‑1 (MQ=255)
                             aTCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGt                                                                                                                        <  1:437385/149‑1 (MQ=255)
                                 cgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaaa                                                                                                                    <  4:658430/149‑1 (MQ=255)
                                 cgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaaa                                                                                                                    <  3:400066/149‑1 (MQ=255)
                                   gCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGGAAAATg                                                                                                                  >  4:587526/1‑149 (MQ=255)
                                      ttGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                                                                               >  4:484855/1‑149 (MQ=255)
                                      atGGCACCAGGCCCAAACAACATGCTGCGGCGGGTCCGATATTTACGGTGTTGGAGCGAAGCAGAAATCTTTGGGTATCCTTCCATTTCAGGCCCGGGTTCCCGCTGGCCCCGGTCAGCAGGGCTTCCCGTCGTGAACGTAAAAttgtg                                                                                                               >  4:599364/2‑145 (MQ=255)
                                           aCCAGGCACAAACAACATGCTGCAGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAg                                                                                                          >  4:49164/1‑149 (MQ=255)
                                              aGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCCGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc                                                                                                       >  2:68607/1‑149 (MQ=255)
                                               ggCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgcg                                                                                                      >  2:575628/1‑149 (MQ=255)
                                                  acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTTCAGTGCGCCt                                                                                                   <  3:49164/149‑1 (MQ=255)
                                                    aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                                                                                 >  2:535742/1‑149 (MQ=255)
                                                    aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                                                                                 <  2:417767/149‑1 (MQ=255)
                                                    aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                                                                                 >  4:415470/1‑149 (MQ=255)
                                                    aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                                                                                                 <  1:488902/149‑1 (MQ=255)
                                                       caacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTc                                                                                              <  1:407224/149‑1 (MQ=255)
                                                          caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                                                                           <  2:319650/149‑1 (MQ=255)
                                                            tgctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTAGTCATCGa                                                                                         >  2:700786/1‑149 (MQ=255)
                                                             gctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGAc                                                                                        <  3:415470/149‑1 (MQ=255)
                                                                   gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATAGCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg                                                                                  <  1:68607/149‑1 (MQ=255)
                                                                     gggTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                                                                                >  2:356133/1‑149 (MQ=255)
                                                                     gggTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                                                                                <  2:45848/149‑1 (MQ=255)
                                                                           gAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgc                                                                          >  3:97639/1‑149 (MQ=255)
                                                                                    gttgttGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                                 >  1:338547/1‑149 (MQ=255)
                                                                                     ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGtt                                                                >  3:382597/1‑149 (MQ=255)
                                                                                          gCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGc                                                           >  1:403768/1‑149 (MQ=255)
                                                                                            aGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCGCGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCgg                                                         <  4:382597/149‑1 (MQ=255)
                                                                                            aGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCgg                                                         <  4:6863/149‑1 (MQ=255)
                                                                                                    ggAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAAc                                                 <  1:356133/149‑1 (MQ=255)
                                                                                                     gAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcg                                                                           <  4:314357/122‑1 (MQ=255)
                                                                                                     gAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcg                                                                           >  3:314357/1‑122 (MQ=255)
                                                                                                       aaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTc                                              >  1:547283/1‑149 (MQ=255)
                                                                                                       aaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTc                                              <  1:440457/149‑1 (MQ=255)
                                                                                                        aTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTcc                                             <  2:547283/149‑1 (MQ=255)
                                                                                                        aTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTcc                                             >  1:288566/1‑149 (MQ=255)
                                                                                                           aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGa                                          >  4:537096/1‑149 (MQ=255)
                                                                                                              ggCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtg                                       >  4:419534/1‑149 (MQ=255)
                                                                                                              ggCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtg                                       >  3:521745/1‑149 (MQ=255)
                                                                                                               gCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCgg                                                         <  1:542952/130‑1 (MQ=255)
                                                                                                               gCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCgg                                                         >  2:542952/1‑130 (MQ=255)
                                                                                                                 aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                                 >  3:319204/1‑120 (MQ=255)
                                                                                                                 aaTCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt                                                                 <  4:319204/120‑1 (MQ=255)
                                                                                                                   tcaacgcATGGTGGAGCGGGGATACCGCTGGCGCGTGCTAGCAGGCCTTCCAGCCGTACACGTCAAATGCAGTCCAGTGTGACTTTTGCATCAACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCACAATggtggt                                  <  3:599364/143‑1 (MQ=255)
                                                                                                                     ctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCt                                                                                                   >  4:750417/1‑82 (MQ=255)
                                                                                                                     ctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCt                                                                                                   <  3:750417/82‑1 (MQ=255)
                                                                                                                      tccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCg                               >  4:756054/1‑149 (MQ=255)
                                                                                                                      tccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCg                               >  2:459392/1‑149 (MQ=255)
                                                                                                                        ccctGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc                                                                                 >  3:561233/4‑97 (MQ=255)
                                                                                                                        ccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc                                                                                 <  4:561233/97‑1 (MQ=255)
                                                                                                                        ccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGaa                             >  4:161521/1‑149 (MQ=255)
                                                                                                                          aTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAAtt                           >  1:720303/1‑149 (MQ=255)
                                                                                                                               agagCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATc                      <  2:668219/149‑1 (MQ=255)
                                                                                                                                    cGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCaa                 >  3:101193/1‑149 (MQ=255)
                                                                                                                                    cGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCaa                 <  2:288566/149‑1 (MQ=255)
                                                                                                                                     gggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                                                                                            >  2:45778/1‑73 (MQ=255)
                                                                                                                                     gggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                                                                                            <  1:45778/73‑1 (MQ=255)
                                                                                                                                      ggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAAc               >  3:732474/1‑149 (MQ=255)
                                                                                                                                      ggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAAc               >  3:628536/1‑149 (MQ=255)
                                                                                                                                        aTACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTg             >  4:382361/1‑149 (MQ=255)
                                                                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt           <  1:34540/149‑1 (MQ=255)
                                                                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt           <  3:161521/149‑1 (MQ=255)
                                                                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGtt           >  1:574018/1‑149 (MQ=255)
                                                                                                                                           ctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCa                                            <  1:345690/113‑1 (MQ=255)
                                                                                                                                           ctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCa                                            >  2:345690/3‑115 (MQ=255)
                                                                                                                                           ctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt          >  2:174297/3‑149 (MQ=255)
                                                                                                                                           ccGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt          <  3:561116/149‑1 (MQ=255)
                                                                                                                                           ccGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt          <  4:101193/149‑1 (MQ=255)
                                                                                                                                            tgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc                                                                     <  2:269272/88‑1 (MQ=255)
                                                                                                                                            tgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc                                                                     >  1:269272/2‑89 (MQ=255)
                                                                                                                                            cGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTtc         <  1:1929/149‑1 (MQ=255)
                                                                                                                                             gCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTtct        <  1:535742/149‑1 (MQ=255)
                                                                                                                                              cTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTtctc       <  2:574018/149‑1 (MQ=255)
                                                                                                                                                 gcgcGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAAc    <  2:389258/149‑1 (MQ=255)
                                                                                                                                                   gcgTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACgc  >  2:404088/1‑149 (MQ=255)
                                                                                                                                                   gcgTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACgc  <  4:628536/149‑1 (MQ=255)
                                                                                                                                                   gcgTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACgc  <  4:732474/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGC  >  NC_000913/650382‑650677

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