Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,529 G→A 16.9% 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,5290GA16.9% 186.8 / 29.1 83R85R (CGC→CGTcitDcitrate lyase, acyl carrier (gamma) subunit
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (8/6);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

GCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGA  >  NC_000913/650430‑650626
                                                                                                   |                                                                                                 
gCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcg                                                                                                   <  1:580769/100‑1 (MQ=255)
gCACAAACAACATGCTGCGGCGGGTGCGAGATTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcg                                                                                                   >  1:586944/1‑100 (MQ=255)
    aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTg                                                                                               <  1:1692695/100‑1 (MQ=255)
    aaacaaCATGCTGCGGCGGGGGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGcgcg                                                                                                 >  1:1409241/1‑98 (MQ=255)
      acaacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGcc                                                                                             <  1:2444882/100‑1 (MQ=255)
       caacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCca                                                                                            >  2:618691/1‑100 (MQ=255)
        aacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag                                                                                           <  1:926758/100‑1 (MQ=255)
         acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGGGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagc                                                                                          <  1:1820669/100‑1 (MQ=255)
         acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGTAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                          >  1:1843987/1‑100 (MQ=255)
         acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                          >  1:552242/1‑100 (MQ=255)
          caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATGCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                                         >  2:369100/1‑100 (MQ=255)
          caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                                         >  2:1216968/1‑100 (MQ=255)
          caTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                          >  2:1933962/1‑99 (MQ=255)
            tgctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                                                                                       >  1:136263/1‑100 (MQ=255)
            tgctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggt                                                                                       <  2:140984/100‑4 (MQ=255)
            tgctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggt                                                                                       <  2:653733/100‑4 (MQ=255)
             gctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAggg                                                                                      <  2:416746/100‑1 (MQ=255)
             gctgcGGCGGGTGCGAGTTTTACGTTGTGGCGGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAggg                                                                                      <  1:1685967/100‑1 (MQ=255)
             gctgaGGCGGGTACGAGTTTTACGCAGTTGCATCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAggg                                                                                      <  1:1279877/100‑1 (MQ=255)
               tgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggtatc                                                                                    <  2:503391/100‑7 (MQ=255)
                gcgtcgGGTGCGAGTTTTACGTTCTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggtatc                                                                                    <  1:1696591/99‑7 (MQ=255)
                gcggcgGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggtatcc                                                                                   <  2:933332/100‑8 (MQ=255)
                 cggcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                  <  2:1692273/100‑1 (MQ=255)
                 cggcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTc                                                                                  <  1:1332948/100‑1 (MQ=255)
                  ggcggGTGCGAGTTTTACGTTGTTGTAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTcc                                                                                 <  2:1411782/100‑1 (MQ=255)
                   gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCa                                                                                <  1:799922/100‑1 (MQ=255)
                    cggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAg                                                                               <  1:1393292/100‑1 (MQ=255)
                       gTGCGTGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCGCTCCTCCCTTGCCCTACCCGGGATACTGCCGGCGCGTGCCAGCAGGGCTTCCAGTCg                                                                            <  1:1959460/100‑1 (MQ=255)
                       gTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCg                                                                            >  2:997223/1‑100 (MQ=255)
                        tGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCg                                                                            <  2:322831/99‑1 (MQ=255)
                        tGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                                           >  1:1036563/1‑100 (MQ=255)
                          cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         <  2:920931/100‑1 (MQ=255)
                          cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  2:610923/1‑100 (MQ=255)
                          cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  1:118638/1‑100 (MQ=255)
                          cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  2:821779/1‑100 (MQ=255)
                          cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  2:1761013/1‑100 (MQ=255)
                                 tACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaa                                                                  <  1:1017160/100‑1 (MQ=255)
                                  aCGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTaaaa                                                                 >  2:874547/1‑100 (MQ=255)
                                   cGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                                >  2:2218115/1‑100 (MQ=255)
                                     ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                              <  2:239953/100‑1 (MQ=255)
                                     ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                              <  1:854134/100‑1 (MQ=255)
                                     ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                              >  1:997636/1‑100 (MQ=255)
                                       gttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                            <  1:1464694/100‑1 (MQ=255)
                                       gttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                            >  2:2050993/1‑100 (MQ=255)
                                        ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCCTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                                           <  1:2309864/100‑1 (MQ=255)
                                        ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCACCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                                           >  2:1065106/1‑100 (MQ=255)
                                         tgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                                                          >  2:2573730/1‑100 (MQ=255)
                                         tgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                                                          >  2:71111/1‑100 (MQ=255)
                                         tgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                                                          >  1:2468805/1‑100 (MQ=255)
                                          gCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                                         <  1:1051034/100‑1 (MQ=255)
                                             gCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGt                                                      >  1:1117413/1‑100 (MQ=255)
                                               gAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTTCACGTAAAATGCAGTCCAGTgc                                                    <  1:1500056/100‑1 (MQ=255)
                                               gAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc                                                    <  1:1323698/100‑1 (MQ=255)
                                                aaGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgcg                                                   <  2:1605095/100‑1 (MQ=255)
                                                      aaaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTg                                             <  1:1276462/100‑1 (MQ=255)
                                                      aaaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTg                                             <  1:748277/100‑1 (MQ=255)
                                                       aaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGt                                            <  2:1198567/100‑1 (MQ=255)
                                                          cATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                                         <  1:2050980/100‑1 (MQ=255)
                                                           aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                        <  2:708137/100‑1 (MQ=255)
                                                           aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                        >  2:270240/1‑100 (MQ=255)
                                                            ttGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCg                                       >  1:876330/1‑100 (MQ=255)
                                                             tGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGa                                      >  1:1151100/1‑100 (MQ=255)
                                                               gCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACa                                    <  1:115943/100‑1 (MQ=255)
                                                                cAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAt                                   >  1:1036059/1‑100 (MQ=255)
                                                                    cctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg                               >  1:2102721/1‑100 (MQ=255)
                                                                    cctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg                               >  2:2418906/1‑100 (MQ=255)
                                                                       cccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTg                            >  1:770261/1‑100 (MQ=255)
                                                                            ggCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgc                       >  1:736477/1‑100 (MQ=255)
                                                                                  gCCGGGATACCGCTTGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACATAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCa                 <  1:2192997/100‑1 (MQ=255)
                                                                                    cGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACg               >  2:774171/1‑100 (MQ=255)
                                                                                       gATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTg            <  1:1038176/100‑1 (MQ=255)
                                                                                        gccctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGt           >  2:1865014/6‑100 (MQ=255)
                                                                                         tACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa          >  2:1021632/1‑100 (MQ=255)
                                                                                         tACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa          >  2:904732/1‑100 (MQ=255)
                                                                                         ccctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa          >  2:468087/5‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTGCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  2:1572668/4‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACTCCGCGCACGTTGTAg         >  1:1360277/4‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  2:998905/4‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  2:1797891/4‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  1:589880/4‑100 (MQ=255)
                                                                                          cctgCTGGCACGCGCCAGCAGGGCTTACAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACCTTCAGCTGTACGCCGCGCACGTTGTAg         >  2:1577364/4‑100 (MQ=255)
                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGt           <  1:985916/98‑1 (MQ=255)
                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         <  2:772186/100‑1 (MQ=255)
                                                                                          aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         <  2:246330/100‑1 (MQ=255)
                                                                                            cGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCg       <  1:231810/100‑1 (MQ=255)
                                                                                              cTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCggg     <  2:3829/100‑1 (MQ=255)
                                                                                                 gcgcGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATCCAGCTGTACGCCGCGCACGTTGTAGCGGGCga  <  2:1036566/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGA  >  NC_000913/650430‑650626

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