Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 648,378 T→C E139E (GAA→GAG citF ← citrate lyase alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913648,3780TC99.3% 879.6 / ‑3.5 284E139E (GAA→GAGcitFcitrate lyase alpha chain
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (148/133);  total (149/135)
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.63e-01

Only 100 of 296 total aligned reads displayed.

GGACGGGACGCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTG  >  NC_000913/648288‑648467
                                                                                          |                                                                                         
ggacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTc                                                                                         <  2:464093/93‑1 (MQ=255)
ggacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTc                                                                                         >  1:231500/1‑93 (MQ=255)
ggacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTc                                                                                         <  2:3245069/93‑1 (MQ=255)
 gacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGctct                                                                                        >  1:6895907/1‑93 (MQ=255)
 gacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGctct                                                                                        <  1:638528/93‑1 (MQ=255)
 gacgggacgCCGAGGAAAGCCACGTCGATATGCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGctct                                                                                        <  1:4419060/93‑1 (MQ=255)
    gggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGcc                                                                                     >  2:1801929/1‑93 (MQ=255)
       acgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCcagc                                                                                  <  1:326607/93‑1 (MQ=255)
        cgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCcagca                                                                                 <  2:6690844/93‑1 (MQ=255)
         gCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCcagcag                                                                                >  1:3797636/1‑93 (MQ=255)
           cGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGAc                                                                              <  2:4227485/93‑1 (MQ=255)
            gAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGAcc                                                                             >  1:3394988/1‑93 (MQ=255)
            gAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGAcc                                                                             >  2:2169282/1‑93 (MQ=255)
             aGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCa                                                                            >  1:813634/1‑93 (MQ=255)
              ggAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCAc                                                                           >  1:3371416/1‑93 (MQ=255)
              ggAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCAc                                                                           >  1:1483797/1‑93 (MQ=255)
               gAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACg                                                                          >  2:637260/1‑93 (MQ=255)
                aaaGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACgg                                                                         <  1:6492224/93‑1 (MQ=255)
                aaaGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACgg                                                                         <  1:6492227/93‑1 (MQ=255)
                  aGCCACGTCGATATTCAGTTCGCCGCTCTGTCCCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGga                                                                       <  1:3588403/93‑1 (MQ=255)
                  aGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGga                                                                       <  2:1013409/93‑1 (MQ=255)
                    ccACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGgaga                                                                     <  2:7520049/93‑1 (MQ=255)
                     cACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGAt                                                                    <  1:3674682/93‑1 (MQ=255)
                     cACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGAt                                                                    <  2:2717655/93‑1 (MQ=255)
                      aCGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGAtc                                                                   >  2:4000378/1‑93 (MQ=255)
                       cctCGATAATCAGTTCCGCGGTCTGGACCAGATTGACCCCACCCGCGGGAGGAGGGGTCCGCCCCGGGTCTGCCAGCAGACCACGGGAGAtct                                                                  <  2:4985964/91‑1 (MQ=255)
                       cGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGAtct                                                                  >  1:2542553/1‑93 (MQ=255)
                         tCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATctct                                                                >  1:2916155/1‑93 (MQ=255)
                         tCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATctct                                                                >  1:6209793/1‑93 (MQ=255)
                          cGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTTAGAGTGGATCTGCACCGCCTCTGCCAGCAGACCACCGGCGATCTCtt                                                               >  2:2345882/1‑93 (MQ=255)
                            atatTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTcc                                                             <  1:866635/93‑1 (MQ=255)
                             tatTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCg                                                            >  2:3873078/1‑93 (MQ=255)
                             tatTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCg                                                            <  2:6319614/93‑1 (MQ=255)
                             tatTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGAGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCg                                                            >  2:2614799/1‑93 (MQ=255)
                              atTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGc                                                           >  1:1672393/1‑93 (MQ=255)
                               ttCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGCTCTCTTCCGcc                                                          >  2:6121174/1‑93 (MQ=255)
                               ttCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGcc                                                          >  2:5422089/1‑93 (MQ=255)
                                 cAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAg                                                        <  2:4549315/93‑1 (MQ=255)
                                  aGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGt                                                       >  1:3912093/1‑93 (MQ=255)
                                         cgcTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCAc                                                >  2:2488075/1‑93 (MQ=255)
                                          gcTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACg                                               >  1:6659074/1‑93 (MQ=255)
                                              tGTACGAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCGCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAgg                                           <  2:3945741/93‑1 (MQ=255)
                                              tGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAgg                                           >  2:6425888/1‑93 (MQ=255)
                                               gTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGc                                          <  2:4838514/93‑1 (MQ=255)
                                               gTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGc                                          <  2:5974792/93‑1 (MQ=255)
                                                tACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGcg                                         >  1:7374371/1‑92 (MQ=255)
                                                 aCCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCg                                        <  2:3545959/93‑1 (MQ=255)
                                                  ccAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCgg                                       >  2:2117788/1‑93 (MQ=255)
                                                  ccAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCgg                                       >  1:4079954/1‑93 (MQ=255)
                                                   cAGATGCCCACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGa                                      <  2:1630187/93‑1 (MQ=255)
                                                   cAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGa                                      <  1:5600590/93‑1 (MQ=255)
                                                   cAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGa                                      <  1:7662534/93‑1 (MQ=255)
                                                   cAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGa                                      >  1:7743165/1‑93 (MQ=255)
                                                    aGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAg                                     <  1:5888894/93‑1 (MQ=255)
                                                     gATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAgg                                    >  2:1145479/1‑93 (MQ=255)
                                                      aTGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGt                                   >  1:128419/1‑93 (MQ=255)
                                                       tGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGta                                  <  1:5111538/93‑1 (MQ=255)
                                                       tGCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGta                                  >  1:658429/1‑93 (MQ=255)
                                                        gCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGtat                                 >  2:6564195/1‑93 (MQ=255)
                                                        gCACACGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGtat                                 <  2:705275/93‑1 (MQ=255)
                                                          acacGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGCATaa                               >  1:316207/1‑93 (MQ=255)
                                                          acacGACCGCCGTGAGAGTGGATCTGCAACGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCGGGAGGTATaa                               >  1:1411907/1‑93 (MQ=255)
                                                           cacGACCGCCGTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATaaa                              <  2:792791/93‑1 (MQ=255)
                                                           cacGACCGCCGTGAGAGTGGATCCGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATaaa                              <  1:6758148/93‑1 (MQ=255)
                                                               gccgccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCg                          <  2:2421647/92‑1 (MQ=255)
                                                                ccgccgTGAGAGTGTATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCgg                         <  1:4800240/93‑1 (MQ=255)
                                                                ccgccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCgg                         >  1:2136774/1‑93 (MQ=255)
                                                                ccgccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCgg                         <  1:1088470/93‑1 (MQ=255)
                                                                 cgccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCggg                        >  1:4501551/1‑93 (MQ=255)
                                                                  gccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGt                       >  1:6169589/1‑93 (MQ=255)
                                                                  gccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGt                       <  1:4155098/93‑1 (MQ=255)
                                                                  gccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGt                       <  2:2958148/93‑1 (MQ=255)
                                                                   ccgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTa                      <  2:1574767/93‑1 (MQ=255)
                                                                    cgTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTaa                     >  2:3804921/1‑93 (MQ=255)
                                                                    cgTGAGAGGGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTaa                     >  1:3092544/1‑93 (MQ=255)
                                                                     gTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAAc                    >  1:1822004/1‑93 (MQ=255)
                                                                     gTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAAc                    >  2:7672002/1‑93 (MQ=255)
                                                                     gTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAAc                    >  1:7168473/1‑93 (MQ=255)
                                                                     gTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAAc                    <  1:2443711/93‑1 (MQ=255)
                                                                     gTGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGGAAc                    >  2:5649395/1‑93 (MQ=255)
                                                                      tGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAAcc                   >  1:4366214/1‑93 (MQ=255)
                                                                      tGAGAGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAAcc                   >  1:6067577/1‑93 (MQ=255)
                                                                       gagaGTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCa                  >  2:1482096/1‑93 (MQ=255)
                                                                        agagTGGATCTGCACCGGCTCTGCCAGCGGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCAc                 <  2:4388492/93‑1 (MQ=255)
                                                                        agagTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCAc                 >  2:6733087/1‑93 (MQ=255)
                                                                           gTGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGcc              <  1:1984163/93‑1 (MQ=255)
                                                                            tGGATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGccc             >  1:7430133/1‑93 (MQ=255)
                                                                             ggATCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCt            <  1:5441587/93‑1 (MQ=255)
                                                                               aTCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTgg          >  2:7018321/1‑93 (MQ=255)
                                                                                tCTGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGc         >  1:189799/1‑93 (MQ=255)
                                                                                  tGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGa       >  2:5218540/1‑93 (MQ=255)
                                                                                  tGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGa       >  1:1577680/1‑93 (MQ=255)
                                                                                  tGCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGa       >  1:7240056/1‑93 (MQ=255)
                                                                                   gCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGaa      >  1:4634802/1‑93 (MQ=255)
                                                                                   gCACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGaa      <  1:3859931/93‑1 (MQ=255)
                                                                                    cACCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAt     <  2:7596397/93‑1 (MQ=255)
                                                                                    cACCGGCTCTGCAAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAt     <  2:7317228/93‑1 (MQ=255)
                                                                                     aCCGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAtg    >  2:3069549/1‑93 (MQ=255)
                                                                                      ccGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAtgt   <  1:7852155/93‑1 (MQ=255)
                                                                                       cGGCTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAtgtg  <  1:3489258/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGACGGGACGCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTG  >  NC_000913/648288‑648467

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