Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 647,331 G→A D488D (GAC→GAT citF ← citrate lyase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913647,3310GA99.3% 475.9 / ‑5.1 139D488D (GAC→GATcitFcitrate lyase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (69/69);  total (70/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 150 total aligned reads displayed.

TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTT  >  NC_000913/647241‑647419
                                                                                          |                                                                                        
ttCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATc                                                                                        <  2:727740/93‑1 (MQ=255)
 tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCt                                                                                       >  2:3749067/1‑93 (MQ=255)
 tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCt                                                                                       >  1:3334077/1‑93 (MQ=255)
   aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                     <  1:1886084/93‑1 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtg                                                                                    <  2:166171/93‑1 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGACATCGCGGTAACGCACAACAGCAACAACGCGATCtgtg                                                                                    >  2:2119838/1‑93 (MQ=255)
     gAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGa                                                                                   <  1:3102175/93‑1 (MQ=255)
        gCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCGATCGCGGTAACGCACAACGGCAACGACGCGCTCTGTGAAtt                                                                                >  1:2525027/1‑93 (MQ=255)
        gCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAAtt                                                                                <  2:385530/93‑1 (MQ=255)
          aGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCa                                                                              >  1:1027200/1‑93 (MQ=255)
          aGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGCGCATTCa                                                                              >  1:168423/1‑93 (MQ=255)
           ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                             <  1:2946554/93‑1 (MQ=255)
            gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                            <  2:2788715/93‑1 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                           <  2:604921/93‑1 (MQ=255)
              gCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCg                                                                          >  2:1395517/1‑93 (MQ=255)
               cATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCCACGACGCGATCTGTGAATTCACTCgg                                                                         >  2:1522277/1‑93 (MQ=255)
               cATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCgg                                                                         <  1:3741457/93‑1 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCggt                                                                        >  2:3182114/1‑92 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                        <  1:2288258/93‑1 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                        >  2:1661561/1‑93 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                        <  2:1788558/93‑1 (MQ=255)
                  ggCTTATTCCTTCACCTGATGCACAACATCGATCCCCGAGCCATCGCGGTAACGCACAACGGCAACGCCGCGATCTGTGAATTCAATCGGCtg                                                                      >  1:4470550/1‑93 (MQ=255)
                     ttATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTgggg                                                                   >  1:695471/1‑93 (MQ=255)
                     ttATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTgg                                                                   >  2:2799736/1‑93 (MQ=255)
                       aTTCCTTCTCCTGATGGACAACAGCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                 <  1:330856/93‑1 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                 >  2:1886469/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                 >  1:1664610/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                 >  2:2471740/1‑93 (MQ=255)
                        ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                >  2:253646/1‑93 (MQ=255)
                        ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                >  2:1989988/1‑93 (MQ=255)
                          ccTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAc                                                              <  1:2013279/93‑1 (MQ=255)
                           cTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAcc                                                             >  1:4322170/1‑93 (MQ=255)
                           cTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAcc                                                             >  1:848782/1‑93 (MQ=255)
                            ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                            >  2:1287429/1‑93 (MQ=255)
                             tCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCgg                                                           >  1:1098184/1‑93 (MQ=255)
                             tCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATAGGCTGTGGTTCACCgg                                                           >  2:4070266/1‑93 (MQ=255)
                              cACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGt                                                          >  2:2024954/1‑93 (MQ=255)
                               cccTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                         <  1:3503618/92‑1 (MQ=255)
                               aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                         <  2:4414998/93‑1 (MQ=255)
                               aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                         <  2:4414997/93‑1 (MQ=255)
                               aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                         >  1:3660403/1‑93 (MQ=255)
                                  tGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                      <  1:3524141/93‑1 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                     >  1:4060719/1‑93 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                     <  1:2501359/93‑1 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                     <  2:519322/93‑1 (MQ=255)
                                     tGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGa                                                   <  1:2748496/93‑1 (MQ=255)
                                      gCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAc                                                  <  1:4012491/93‑1 (MQ=255)
                                       cacaACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcg                                                 <  1:1009254/93‑1 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                <  1:2523205/93‑1 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                >  1:1352922/1‑93 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                <  1:324270/93‑1 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAACGGTCAGCAGAcgc                                                >  1:2536158/1‑93 (MQ=255)
                                          aacaTCGACCACCGAGCCATCGGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCa                                             <  1:2634066/93‑1 (MQ=255)
                                           acaTCTGTCACCGACCCATCGCGGTAACGCACAACGGCATCGACGCGGTCTGTGAATTCAATCGGCTGTGCCTCACCGGTCAGCAGACGCGCa                                             <  1:3864980/93‑1 (MQ=255)
                                           acaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCa                                             <  2:3206879/93‑1 (MQ=255)
                                            caTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCAc                                            <  1:2561120/93‑1 (MQ=255)
                                             aTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACg                                           >  2:4676288/1‑93 (MQ=255)
                                              tCGATCCCCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                          <  1:4133162/93‑1 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                          <  1:294967/93‑1 (MQ=255)
                                               cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                         <  2:2391992/93‑1 (MQ=255)
                                               cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                         <  1:595866/93‑1 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                        >  1:796533/1‑93 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                        <  1:1767774/93‑1 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATAGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                        >  2:1446187/1‑93 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAGCGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                       <  2:4149668/93‑1 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                       <  2:3568680/93‑1 (MQ=255)
                                                   cACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgcg                                     >  2:4691026/1‑93 (MQ=255)
                                                   cACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgcg                                     >  1:100402/1‑93 (MQ=255)
                                                    aCCGAGCCATCGCGGTAACGCACAACGGCAACGAGGCGATATGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCgcgc                                    <  2:2864950/93‑1 (MQ=255)
                                                    aCCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCgcgc                                    >  1:1992824/1‑93 (MQ=255)
                                                       gAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGc                                 <  1:4608596/93‑1 (MQ=255)
                                                        aGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGcc                                >  2:309602/1‑93 (MQ=255)
                                                           cATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACt                             >  1:2348107/1‑93 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCCCTc                            >  1:3380013/1‑93 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                            <  2:2963939/93‑1 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                            <  2:3736994/93‑1 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                            <  1:1115442/93‑1 (MQ=255)
                                                               gcgGTAACGCACAACGGCCACGACGCGATCTGTGAATTCAATCGGCTGTGGTTTACCGGTCAGCAGACGCGCCCGTTCGCGCAGCCACTCAAt                         >  2:2738861/1‑93 (MQ=255)
                                                               gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                         >  1:4685566/1‑93 (MQ=255)
                                                                 ggTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATgg                       >  1:4501300/1‑93 (MQ=255)
                                                                    aaCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGaaa                    >  2:2399865/1‑93 (MQ=255)
                                                                    aaCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGaaa                    <  1:3623007/93‑1 (MQ=255)
                                                                    aaCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGaaa                    >  2:3697510/1‑93 (MQ=255)
                                                                     aCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAAc                   <  2:404849/93‑1 (MQ=255)
                                                                        cacaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCAc                >  1:674852/1‑93 (MQ=255)
                                                                        cacaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCAc                <  2:769529/93‑1 (MQ=255)
                                                                         acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt               <  1:4412574/93‑1 (MQ=255)
                                                                         acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt               >  1:2201833/1‑93 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt              >  1:3288057/1‑93 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt              <  1:1441936/93‑1 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTAAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt              <  2:1752864/93‑1 (MQ=255)
                                                                           aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt             >  1:1052815/1‑93 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt          >  2:351414/1‑93 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt          <  1:4589950/93‑1 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt          <  1:3876150/93‑1 (MQ=255)
                                                                                cAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGc        <  1:370834/93‑1 (MQ=255)
                                                                                 aaCGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGcc       <  2:2191983/93‑1 (MQ=255)
                                                                                  acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc      >  2:1915713/1‑93 (MQ=255)
                                                                                      cgcgATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCtt  <  1:4137030/93‑1 (MQ=255)
                                                                                      agcgATCGGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCtt  <  2:2496763/92‑1 (MQ=255)
                                                                                          |                                                                                        
TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTT  >  NC_000913/647241‑647419

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