Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,532 T→C 16.6% A84A (GCA→GCG 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,5320TC16.6% 184.7 / 28.3 84A84A (GCA→GCGcitDcitrate lyase, acyl carrier (gamma) subunit
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (8/6);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

AAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGA  >  NC_000913/650434‑650626
                                                                                                  |                                                                                              
aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTg                                                                                               <  1:1692695/100‑1 (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:653733/100‑4 (MQ=255)
        tgctgcGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggt                                                                                       <  2:140984/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:1761013/1‑100 (MQ=255)
                      cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  2:610923/1‑100 (MQ=255)
                      cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  1:118638/1‑100 (MQ=255)
                      cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         <  2:920931/100‑1 (MQ=255)
                      cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc                                                                         >  2:821779/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                                                              <  1:854134/100‑1 (MQ=255)
                                 ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                              <  2:239953/100‑1 (MQ=255)
                                 ttgttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                                              >  1:997636/1‑100 (MQ=255)
                                   gttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                            >  2:2050993/1‑100 (MQ=255)
                                   gttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAg                                                            <  1:1464694/100‑1 (MQ=255)
                                    ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCCTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                                           <  1:2309864/100‑1 (MQ=255)
                                    ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCACCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                                           >  2:1065106/1‑100 (MQ=255)
                                     tgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                                                          >  2:71111/1‑100 (MQ=255)
                                     tgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                                                          >  2:2573730/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:748277/100‑1 (MQ=255)
                                                  aaaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTg                                             <  1:1276462/100‑1 (MQ=255)
                                                   aaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGt                                            <  2:1198567/100‑1 (MQ=255)
                                                      cATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                                         <  1:2050980/100‑1 (MQ=255)
                                                       aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                        >  2:270240/1‑100 (MQ=255)
                                                       aTTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                        <  2:708137/100‑1 (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                               >  2:2418906/1‑100 (MQ=255)
                                                                cctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg                               >  1:2102721/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:904732/1‑100 (MQ=255)
                                                                                     tACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTa          >  2:1021632/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:1797891/4‑100 (MQ=255)
                                                                                      cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  1:589880/4‑100 (MQ=255)
                                                                                      cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         >  2:998905/4‑100 (MQ=255)
                                                                                      cctgCTGGCACGCGCCAGCAGGGCTTACAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACCTTCAGCTGTACGCCGCGCACGTTGTAg         >  2:1577364/4‑100 (MQ=255)
                                                                                      aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGt           <  1:985916/98‑1 (MQ=255)
                                                                                      aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         <  2:246330/100‑1 (MQ=255)
                                                                                      aCCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg         <  2:772186/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)
                                                                                                  |                                                                                              
AAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGA  >  NC_000913/650434‑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: