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.5% 697.0 / ‑5.5 194D488D (GAC→GATcitFcitrate lyase alpha chain
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (95/97);  total (97/97)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 200 total aligned reads displayed.

TCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647242‑647420
                                                                                         |                                                                                         
tCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCt                                                                                        <  2:2147687/93‑1 (MQ=255)
  aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      >  1:960633/1‑93 (MQ=255)
  aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      <  1:1020399/93‑1 (MQ=255)
  aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      >  2:2144180/1‑93 (MQ=255)
  aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      <  2:561580/93‑1 (MQ=255)
  aGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgt                                                                                      >  2:3277477/1‑93 (MQ=255)
     aaGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCtgtgca                                                                                   >  2:420805/1‑91 (MQ=255)
      aGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAAt                                                                                  <  1:2418992/93‑1 (MQ=255)
          ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              <  2:3705942/93‑1 (MQ=255)
          ggTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCaa                                                                              >  2:221287/1‑93 (MQ=255)
           gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             <  1:1635560/93‑1 (MQ=255)
           gTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAAt                                                                             <  2:2517873/93‑1 (MQ=255)
            tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            >  1:4648427/1‑93 (MQ=255)
            tGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATc                                                                            <  2:4428373/93‑1 (MQ=255)
             gCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCg                                                                           >  2:4352170/1‑93 (MQ=255)
               aTGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGc                                                                         <  1:3858627/93‑1 (MQ=255)
                tgtcTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCt                                                                        <  2:3264172/90‑1 (MQ=255)
                 ggCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtg                                                                       >  1:1773946/1‑93 (MQ=255)
                   cTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgtg                                                                     >  1:2702044/1‑93 (MQ=255)
                   cTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCtgtg                                                                     <  1:4663125/93‑1 (MQ=255)
                      aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGtt                                                                  <  2:2343335/93‑1 (MQ=255)
                      aTTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTAAATCGGCTGTGGtt                                                                  >  2:3457969/1‑93 (MQ=255)
                       ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                 <  1:2829278/93‑1 (MQ=255)
                       ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                 >  2:1170505/1‑93 (MQ=255)
                       ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                 >  2:4573943/1‑93 (MQ=255)
                       ttCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTc                                                                 >  1:2192149/1‑93 (MQ=255)
                        tCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCa                                                                >  1:1698897/1‑93 (MQ=255)
                        tCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCa                                                                >  1:3519415/1‑93 (MQ=255)
                         ccTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAc                                                               <  2:2750172/93‑1 (MQ=255)
                         ccTTCACCTGACGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCAc                                                               <  2:3807069/93‑1 (MQ=255)
                           ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             <  2:3499773/93‑1 (MQ=255)
                           ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             >  2:2277751/1‑93 (MQ=255)
                           ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             >  2:2180764/1‑93 (MQ=255)
                           ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             >  2:3423112/1‑93 (MQ=255)
                           ttCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCg                                                             <  2:531013/93‑1 (MQ=255)
                            tCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTTTGGTTCACCgg                                                            >  2:4680147/1‑93 (MQ=255)
                            tCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCgg                                                            <  1:700475/93‑1 (MQ=255)
                              aCCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTc                                                          >  2:140297/1‑93 (MQ=255)
                               ccTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTca                                                         <  1:3382485/93‑1 (MQ=255)
                                cTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcag                                                        >  1:1308257/1‑93 (MQ=255)
                                cTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcag                                                        <  2:3587258/93‑1 (MQ=255)
                                cTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcag                                                        >  2:3738627/1‑93 (MQ=255)
                                cTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGGTCACCGGTcag                                                        >  1:3263996/1‑93 (MQ=255)
                                 tGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                       <  1:772086/93‑1 (MQ=255)
                                 tGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                       <  1:848441/93‑1 (MQ=255)
                                 tGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagc                                                       <  2:1044835/93‑1 (MQ=255)
                                  gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      >  2:4236224/1‑93 (MQ=255)
                                  gATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagca                                                      <  2:2974910/93‑1 (MQ=255)
                                   aTGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTcagcag                                                     >  2:4120799/1‑93 (MQ=255)
                                       acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgc                                                 >  1:1843806/1‑93 (MQ=255)
                                       acaacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCACCAGAcgc                                                 >  1:3696700/1‑93 (MQ=255)
                                        caagaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgcg                                                >  2:358355/1‑93 (MQ=255)
                                        caacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGAcgcg                                                <  2:1741925/93‑1 (MQ=255)
                                         aacaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACgcgc                                               >  1:664634/1‑93 (MQ=255)
                                          acaTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCa                                              <  2:1933530/93‑1 (MQ=255)
                                            aTCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACg                                            >  1:1040095/1‑93 (MQ=255)
                                             tccaTCACCGGCCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGAGACCAGACGCACCCGt                                           >  2:644918/4‑93 (MQ=255)
                                             tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           <  1:3088113/93‑1 (MQ=255)
                                             tCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGt                                           <  2:2926998/93‑1 (MQ=255)
                                              cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                          <  1:2436908/93‑1 (MQ=255)
                                              cGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGtt                                          <  2:2792878/93‑1 (MQ=255)
                                               gATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCAGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTc                                         >  1:3040745/1‑93 (MQ=255)
                                                aTCACCGCGCCATCGCGGTAACGCACAACGGCAACGACGCGCTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGGACGTTcg                                        >  2:895734/1‑93 (MQ=255)
                                                aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGCTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        >  1:2011874/1‑93 (MQ=255)
                                                aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        >  1:2546471/1‑93 (MQ=255)
                                                aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        <  2:670714/93‑1 (MQ=255)
                                                aTCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcg                                        <  1:3747326/93‑1 (MQ=255)
                                                 tCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgc                                       <  2:61027/93‑1 (MQ=255)
                                                 tCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTcgc                                       >  2:3955411/1‑93 (MQ=255)
                                                     cGAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAg                                   <  1:4502254/93‑1 (MQ=255)
                                                      gAGCCATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGGATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGc                                  >  2:2322661/1‑93 (MQ=255)
                                                         ccATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCAc                               <  1:4381535/93‑1 (MQ=255)
                                                         ccATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCAc                               <  2:352551/93‑1 (MQ=255)
                                                         ccATCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCAc                               >  2:1694373/1‑93 (MQ=255)
                                                           agcgcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             <  2:1231871/91‑1 (MQ=255)
                                                           aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             >  1:1109823/1‑93 (MQ=255)
                                                           aTCGCGGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTc                             >  2:3655383/1‑93 (MQ=255)
                                                           aTCGCGGGAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACTGGTCAGCAGACGCGCACGTTCTCGCAGACtctc                             >  2:3024625/1‑93 (MQ=255)
                                                             cgcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCaa                           <  1:523785/93‑1 (MQ=255)
                                                              gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                          <  2:96916/93‑1 (MQ=255)
                                                              gcgGTAACGCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAAt                          <  2:2875381/93‑1 (MQ=255)
                                                                      gCACAACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCa                  >  1:1426018/1‑93 (MQ=255)
                                                                       cacaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCAc                 <  2:1793716/93‑1 (MQ=255)
                                                                        acaACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACt                <  1:1216345/93‑1 (MQ=255)
                                                                         caACGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACtt               <  1:2974894/93‑1 (MQ=255)
                                                                          aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt              >  2:1280100/1‑93 (MQ=255)
                                                                          aaCGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACttt              >  1:593751/1‑93 (MQ=255)
                                                                            cGGCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTaa            <  1:2912/93‑1 (MQ=255)
                                                                             ggCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAAt           <  1:3663285/93‑1 (MQ=255)
                                                                              gCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATg          >  2:89067/1‑93 (MQ=255)
                                                                              gCAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATg          >  1:2281448/1‑93 (MQ=255)
                                                                               cAACGACGCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGc         <  2:4153612/93‑1 (MQ=255)
                                                                                 acgacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGccc       <  2:4601977/93‑1 (MQ=255)
                                                                                   gacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGc     >  1:651269/1‑93 (MQ=255)
                                                                                   gacgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGc     <  1:4329180/93‑1 (MQ=255)
                                                                                    acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    <  1:1136407/93‑1 (MQ=255)
                                                                                    acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    >  1:785049/1‑93 (MQ=255)
                                                                                    acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    <  2:1681017/93‑1 (MQ=255)
                                                                                    acgCGATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCt    >  2:2297988/1‑93 (MQ=255)
                                                                                     cgcgATCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCtt   >  2:1555215/1‑93 (MQ=255)
                                                                                         |                                                                                         
TCAGGAAGCAGGTGCATGGCTTATTCCTTCACCTGATGCACAACATCGATCACCGAGCCATCGCGGTAACGCACAACGGCAACGACGCGGTCTGTGAATTCAATCGGCTGTGGTTCACCGGTCAGCAGACGCGCACGTTCGCGCAGCCACTCAATGGAAACCACTTTAATGCCCGCTTC  >  NC_000913/647242‑647420

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