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.7% 1028.8 / ‑5.7 285D488D (GAC→GATcitFcitrate lyase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (166/118);  total (167/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 298 total aligned reads displayed.

TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647241‑647420
                                                                                          |                                                                                         
ttCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATc                                                                                         >  1:8415650/1‑93 (MQ=255)
ttCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATc                                                                                         >  1:7078079/1‑93 (MQ=255)
 tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCt                                                                                        >  2:9505373/1‑93 (MQ=255)
  cAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtg                                                                                       <  1:7767018/93‑1 (MQ=255)
  cAGGAAGCAGCTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtg                                                                                       <  1:5626264/93‑1 (MQ=255)
   aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      >  2:5994528/1‑93 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtg                                                                                     >  1:1407994/1‑93 (MQ=255)
     gAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGa                                                                                    >  1:6441680/1‑93 (MQ=255)
     gAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGa                                                                                    >  1:4565453/1‑93 (MQ=255)
     gAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGa                                                                                    >  2:2396699/1‑93 (MQ=255)
           ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              >  2:3165386/1‑93 (MQ=255)
           ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              >  1:6781406/1‑93 (MQ=255)
           ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGGAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              >  2:3457002/1‑93 (MQ=255)
            gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             >  2:9683647/1‑93 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  1:5290682/93‑1 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  2:4364347/93‑1 (MQ=255)
              gCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCg                                                                           >  1:8101024/1‑93 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                         >  1:3542322/1‑93 (MQ=255)
                 tGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCt                                                                        >  1:5686732/1‑93 (MQ=255)
                 tGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCt                                                                        >  1:441249/1‑93 (MQ=255)
                   gCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgt                                                                      <  2:1276040/93‑1 (MQ=255)
                    cTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgtg                                                                     <  2:8742052/93‑1 (MQ=255)
                    cTTATTCCTCCCCCTGATGCACAACATCGATCGCCGAGGCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgtg                                                                     <  1:3087436/93‑1 (MQ=255)
                     ttATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTgg                                                                    >  1:3854726/1‑93 (MQ=255)
                      tATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGt                                                                   >  2:2542220/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  <  2:452989/93‑1 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  >  2:6572781/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  <  1:1230939/93‑1 (MQ=255)
                        ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGCTc                                                                 >  1:6049603/1‑93 (MQ=255)
                            ttCACCTGATGCACAACATCGATCACCGAGCCATCTCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             >  1:232018/1‑93 (MQ=255)
                            ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             >  1:2515152/1‑93 (MQ=255)
                            tacACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             <  2:5686731/91‑1 (MQ=255)
                            ctCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACACCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             <  1:2632470/92‑1 (MQ=255)
                             tCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCgg                                                            <  1:6724235/93‑1 (MQ=255)
                              cACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGt                                                           >  2:8896538/1‑93 (MQ=255)
                              cACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGt                                                           <  2:9795333/93‑1 (MQ=255)
                               cccTGATGCACAACTTCGGTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                          <  1:8924031/92‑1 (MQ=255)
                               aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGCATTCAATCGGCTGTGGTTCACCGGTc                                                          >  2:9911225/1‑93 (MQ=255)
                                ccTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTca                                                         <  1:4629439/93‑1 (MQ=255)
                                 cTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcag                                                        >  1:4166634/1‑93 (MQ=255)
                                  tGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                       >  2:3434170/1‑93 (MQ=255)
                                   ggtGCACAACAGTGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      <  1:8560265/91‑1 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      <  2:1803643/93‑1 (MQ=255)
                                    aTGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagcag                                                     <  2:9360385/93‑1 (MQ=255)
                                     tGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGa                                                    >  1:9883083/1‑93 (MQ=255)
                                     agCACAAAATCGAGCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTTAATTCAATCGGCTGTGGTTCACCGGTCAGCAGa                                                    <  1:8316763/92‑1 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                 >  2:626291/1‑93 (MQ=255)
                                         caacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGACGGAGGGTAACCGCGCAACAGAcgcc                                                >  1:360700/1‑92 (MQ=255)
                                           acaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGTTCAGCAGACgcgcc                                              >  2:7251386/1‑92 (MQ=255)
                                           acaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCa                                              <  2:8549500/93‑1 (MQ=255)
                                             aTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACg                                            >  1:1696791/1‑93 (MQ=255)
                                             aTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACg                                            >  2:9089630/1‑93 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTGAGCAGACGCGCACGt                                           >  1:5077106/1‑93 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           >  1:4775162/1‑93 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           >  1:5772483/1‑93 (MQ=255)
                                               cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                          <  2:7164267/93‑1 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                         <  2:4670137/93‑1 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        >  1:8278768/1‑93 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        >  1:5476471/1‑93 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGGTCACCGGTCAGCAGACGCGCACGTTcg                                        >  1:4028051/1‑93 (MQ=255)
                                                  tCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGAGCACGTTcgc                                       >  1:2749069/1‑93 (MQ=255)
                                                   cACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgcg                                      >  2:8019934/1‑93 (MQ=255)
                                                     ccGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCa                                    >  2:1487137/1‑93 (MQ=255)
                                                     ccGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCa                                    <  2:8368455/93‑1 (MQ=255)
                                                      cGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAg                                   <  1:1511602/93‑1 (MQ=255)
                                                      cGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAg                                   <  1:1919007/93‑1 (MQ=255)
                                                       gAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGc                                  >  1:8740393/1‑93 (MQ=255)
                                                          ccATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGcccc                               >  2:3094827/1‑91 (MQ=255)
                                                          ccATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCAc                               >  2:1085298/1‑93 (MQ=255)
                                                           cATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACt                              <  2:1981686/93‑1 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             >  2:6781220/1‑93 (MQ=255)
                                                            aTCGCGGGAAAGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             >  2:8518755/1‑93 (MQ=255)
                                                             tCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCa                            >  1:2673735/1‑93 (MQ=255)
                                                             tCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCa                            >  2:6043012/1‑93 (MQ=255)
                                                               gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                          <  2:5587157/93‑1 (MQ=255)
                                                               gcgGTAACGAACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGGTCGCGCAGCCACTCAAt                          >  2:7511518/1‑93 (MQ=255)
                                                                 ggTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATgg                        <  2:8834568/93‑1 (MQ=255)
                                                                  gTAACGCACACCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGa                       <  1:7918464/93‑1 (MQ=255)
                                                                  gTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGa                       >  1:130163/1‑93 (MQ=255)
                                                                  gTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGa                       <  1:6173821/93‑1 (MQ=255)
                                                                  gTAACGCACAACGGCAACGACGCGAACTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCaaccga                       >  2:421859/1‑89 (MQ=255)
                                                                     aCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAAc                    >  1:2443710/1‑93 (MQ=255)
                                                                      cGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTTGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAAcc                   <  1:9500826/93‑1 (MQ=255)
                                                                      cGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAAcc                   <  1:1123722/93‑1 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCTGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCAcct               >  1:1417194/1‑91 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTTGCGCAGCCACTCAATGGAAACCACtt               >  2:9945358/1‑93 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt               >  2:3992616/1‑93 (MQ=255)
                                                                             cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            >  1:8217862/1‑93 (MQ=255)
                                                                             cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            <  1:5994529/93‑1 (MQ=255)
                                                                             cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            <  1:4248608/93‑1 (MQ=255)
                                                                             cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            <  1:1014638/93‑1 (MQ=255)
                                                                              gtcAACGACGCGAGCTGTGAATTCAATCGGCGGTGGTTCACCGGTCAGCAGACGCGCACGGTCGCGCAGCCACTCAATGGAAACCACTTTAAt           <  1:6862486/91‑1 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt           <  2:2819032/93‑1 (MQ=255)
                                                                               gCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGCTCGCGCAGCCACTCAATGGAAACCACTTTAATg          >  2:6132293/1‑93 (MQ=255)
                                                                                 aaCGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGcc        <  1:7411344/93‑1 (MQ=255)
                                                                                  acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc       <  2:698104/93‑1 (MQ=255)
                                                                                  acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc       >  2:4901546/1‑93 (MQ=255)
                                                                                  acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc       >  1:3755506/1‑93 (MQ=255)
                                                                                     acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    >  1:6960724/1‑93 (MQ=255)
                                                                                     acgCGATCTGTGAATTCAACCGCCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    <  2:7616923/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647241‑647420

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