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.6% 916.4 / NA 248D488D (GAC→GATcitFcitrate lyase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (109/138);  minor base C (1/0);  total (110/138)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647241‑647420
                                                                                          |                                                                                         
ttCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATc                                                                                         <  1:5620127/93‑1 (MQ=255)
ttCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATc                                                                                         <  1:2889998/93‑1 (MQ=255)
 tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCt                                                                                        >  1:3530673/1‑93 (MQ=255)
 tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCAAAACGGCGACGACGCGATCt                                                                                        >  1:5735694/1‑93 (MQ=255)
   aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      >  2:4243478/1‑93 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtg                                                                                     >  1:5030488/1‑93 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtg                                                                                     >  1:2704682/1‑93 (MQ=255)
    ggAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtg                                                                                     <  2:6398819/93‑1 (MQ=255)
     gAAGCAGGTGCATGTCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGa                                                                                    <  2:5627911/93‑1 (MQ=255)
       aGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAAt                                                                                  >  1:1034417/1‑93 (MQ=255)
       aGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAAt                                                                                  <  2:687981/93‑1 (MQ=255)
         cAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTc                                                                                >  1:4905737/1‑93 (MQ=255)
         cAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTc                                                                                <  2:5304546/93‑1 (MQ=255)
          aGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCa                                                                               >  1:3303719/1‑93 (MQ=255)
          aGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCa                                                                               >  2:4767642/1‑93 (MQ=255)
           ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              <  2:566648/93‑1 (MQ=255)
            gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             <  2:1560444/93‑1 (MQ=255)
            gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             >  2:4569488/1‑93 (MQ=255)
            gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             <  1:1654791/93‑1 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  1:5762089/93‑1 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  1:540720/93‑1 (MQ=255)
             tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  2:1079265/93‑1 (MQ=255)
              gCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCg                                                                           <  2:6196048/93‑1 (MQ=255)
               cATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCTAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCgg                                                                          >  2:3834545/1‑93 (MQ=255)
                aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                         <  1:3855291/93‑1 (MQ=255)
                   gCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgt                                                                      <  1:2107835/93‑1 (MQ=255)
                    cTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgtg                                                                     >  1:4753815/1‑93 (MQ=255)
                     ttATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTgg                                                                    >  2:3213608/1‑93 (MQ=255)
                     ttATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTgg                                                                    >  1:1255114/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  >  1:3552364/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  >  2:3092970/1‑93 (MQ=255)
                       aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  >  2:1473161/1‑93 (MQ=255)
                          ccTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAc                                                               >  2:2812591/1‑93 (MQ=255)
                          ccTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAc                                                               <  1:1524176/93‑1 (MQ=255)
                               aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                          >  1:1692366/1‑93 (MQ=255)
                                ccTGAGGCACAACATCGATCACCGAGCCATGGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTca                                                         <  1:5093608/93‑1 (MQ=255)
                                  tGATGCACAACATCGATCACCGAGCCAGCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                       <  2:1391188/93‑1 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      >  1:6248574/1‑93 (MQ=255)
                                   gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTCTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      >  2:3284793/1‑93 (MQ=255)
                                     tGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGa                                                    <  2:2350157/93‑1 (MQ=255)
                                     tGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGa                                                    <  2:929917/93‑1 (MQ=255)
                                      gCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAc                                                   <  2:4361351/93‑1 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                 >  1:6447528/1‑93 (MQ=255)
                                        acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                 >  1:3432550/1‑93 (MQ=255)
                                          aacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACgcgc                                               <  2:6632542/93‑1 (MQ=255)
                                          aacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACgcgc                                               >  2:1651853/1‑93 (MQ=255)
                                           acaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCa                                              <  1:4456829/93‑1 (MQ=255)
                                            caTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCAc                                             >  2:482723/1‑93 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           >  1:609375/1‑93 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           <  1:5444378/93‑1 (MQ=255)
                                              tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           >  1:5198341/1‑93 (MQ=255)
                                               cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                          <  2:4146672/93‑1 (MQ=255)
                                               cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                          <  2:4471978/93‑1 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                         >  1:2307009/1‑93 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                         <  2:2295722/93‑1 (MQ=255)
                                                gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGGTCACCGGTCAGCAGACGCGCACGTTc                                         >  2:873137/1‑93 (MQ=255)
                                                gATCACCGAGCCAGCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                         <  2:1138016/93‑1 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        <  2:1845446/93‑1 (MQ=255)
                                                 aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        <  2:624161/93‑1 (MQ=255)
                                                  tCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgc                                       >  2:758332/1‑93 (MQ=255)
                                                  tCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTCCACCGGTCAGCAGACGCGCACGTTcgc                                       <  2:5384129/93‑1 (MQ=255)
                                                     ccGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCa                                    >  1:6575119/1‑93 (MQ=255)
                                                      cGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAg                                   >  1:3052633/1‑93 (MQ=255)
                                                       gAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGc                                  <  1:2160713/93‑1 (MQ=255)
                                                       gAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAACCGCACGTTCGCGCAGc                                  >  2:2357748/1‑93 (MQ=255)
                                                        aGCCATCGCTGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGcc                                 <  1:2472366/93‑1 (MQ=255)
                                                         gCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTTACCGGTCAGCAGACGCGCACGTTCGCGCAGCCa                                >  1:4232949/1‑93 (MQ=255)
                                                         gCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCa                                >  1:2055562/1‑93 (MQ=255)
                                                          ccATCGCGTTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGTTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCAc                               <  1:2149503/93‑1 (MQ=255)
                                                            aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             >  1:1131694/1‑93 (MQ=255)
                                                             tCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCa                            <  2:261975/93‑1 (MQ=255)
                                                               gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                          <  2:6063391/93‑1 (MQ=255)
                                                               gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                          <  1:6681352/93‑1 (MQ=255)
                                                                cgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATg                         >  1:4695860/1‑93 (MQ=255)
                                                                cgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATg                         <  1:737052/93‑1 (MQ=255)
                                                                 ggTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATgg                        <  2:1606910/93‑1 (MQ=255)
                                                                  gTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGa                       >  1:5867929/1‑93 (MQ=255)
                                                                  gTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGa                       <  1:4369557/93‑1 (MQ=255)
                                                                   tAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGaa                      >  1:6489202/1‑93 (MQ=255)
                                                                    aaCGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGaaa                     >  2:3024085/1‑93 (MQ=255)
                                                                      cGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAAcc                   <  1:1739421/93‑1 (MQ=255)
                                                                       gCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCa                  >  2:192542/1‑93 (MQ=255)
                                                                       gCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCa                  >  2:4215300/1‑93 (MQ=255)
                                                                         acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt                >  2:4505306/1‑93 (MQ=255)
                                                                         acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt                >  1:1437825/1‑93 (MQ=255)
                                                                         acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt                >  2:636293/1‑93 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt               >  2:3984648/1‑93 (MQ=255)
                                                                          caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt               >  1:5993577/1‑93 (MQ=255)
                                                                           aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt              <  1:4501934/93‑1 (MQ=255)
                                                                           aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt              >  2:3553520/1‑93 (MQ=255)
                                                                           aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt              >  2:2178443/1‑93 (MQ=255)
                                                                             cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            <  1:6136154/93‑1 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt           <  1:3274964/93‑1 (MQ=255)
                                                                              ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt           <  2:1678920/93‑1 (MQ=255)
                                                                               gCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGAAAACCCCTTTAATg          >  2:2856538/1‑93 (MQ=255)
                                                                                 aaCGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGcc        <  2:4552527/93‑1 (MQ=255)
                                                                                 aaCGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGcc        >  2:5101568/1‑93 (MQ=255)
                                                                                  acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc       <  2:1034417/93‑1 (MQ=255)
                                                                                   cgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCg      >  2:2668088/1‑93 (MQ=255)
                                                                                     acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    <  2:5986711/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647241‑647420

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