Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,335,881 T→G D17E (GAT→GAG acnA → aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,8810TG95.1% 393.4 / 7.3 122D17E (GAT→GAGacnAaconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (63/53);  total (66/56)
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.92e-01

Only 100 of 122 total aligned reads displayed.

CCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCT  >  NC_000913/1335785‑1335976
                                                                                                |                                                                                               
ccTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGaaa                                                                                              >  1:3351462/1‑100 (MQ=255)
ccTGAAGAGAATCAGGGCTGCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGaaa                                                                                              <  1:3166294/100‑1 (MQ=255)
  tGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAAc                                                                                            <  2:3293720/100‑1 (MQ=255)
  tGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAAc                                                                                            <  2:3846434/100‑1 (MQ=255)
     agagAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCACAGAGAAAACTTa                                                                                         >  2:2076279/1‑100 (MQ=255)
       agaATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTAcc                                                                                       <  2:744504/100‑1 (MQ=255)
        gaATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCa                                                                                      >  1:536851/1‑100 (MQ=255)
         aaTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCac                                                                                     <  2:3261888/100‑1 (MQ=255)
          aTCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCGAAGAGAAAACTCACCact                                                                                    >  2:2855840/1‑100 (MQ=255)
            cAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCactac                                                                                  <  2:3173600/100‑1 (MQ=255)
             aGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCactact                                                                                 <  2:560453/100‑1 (MQ=255)
              gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCActacta                                                                                <  1:2492592/100‑1 (MQ=255)
              gggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCActacta                                                                                >  2:1739279/1‑100 (MQ=255)
               ggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACtactac                                                                               <  1:2259827/100‑1 (MQ=255)
               ggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACtactac                                                                               >  1:23907/1‑100 (MQ=255)
               ggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACtactac                                                                               <  2:2329972/100‑1 (MQ=255)
               ggCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACtactac                                                                               >  1:1311009/1‑100 (MQ=255)
                 cTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAg                                                                             <  2:1948080/100‑1 (MQ=255)
                    cGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAACGACACGTTGCAGGCCAAAGAGAAAACGTACCACTACTACAGCCt                                                                          >  1:3744782/1‑100 (MQ=255)
                     gCAACCCTGTCATTAAGGCGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAACGAGAACACTTACCACTACTACAGCCTg                                                                         >  2:971104/1‑100 (MQ=255)
                      cAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGc                                                                        <  1:402300/100‑1 (MQ=255)
                       aaCCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGcc                                                                       <  1:3715449/100‑1 (MQ=255)
                         cccTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGc                                                                     <  2:1038401/100‑1 (MQ=255)
                         cccTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGc                                                                     >  1:1272153/1‑100 (MQ=255)
                            tGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtg                                                                  <  1:1540049/100‑1 (MQ=255)
                             gTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtgc                                                                 <  1:2845940/100‑1 (MQ=255)
                             gTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtgc                                                                 >  1:978975/1‑100 (MQ=255)
                              tCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgct                                                                <  1:2154984/100‑1 (MQ=255)
                               cATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACATGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgctg                                                               >  2:1890501/1‑100 (MQ=255)
                               cATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtgctg                                                               >  2:97999/1‑100 (MQ=255)
                                aTTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtgctgc                                                              <  1:860162/100‑1 (MQ=255)
                                aTTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTtgctgc                                                              >  1:1066240/1‑100 (MQ=255)
                                  tAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTa                                                            <  1:1090694/100‑1 (MQ=255)
                                  tAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTa                                                            >  1:2261003/1‑100 (MQ=255)
                                   aaGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTaa                                                           >  1:1319679/1‑100 (MQ=255)
                                   aaGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTaa                                                           >  1:2541978/1‑100 (MQ=255)
                                   aaGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTaa                                                           <  2:1111584/100‑1 (MQ=255)
                                      gaggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATc                                                        >  1:1271323/1‑100 (MQ=255)
                                        ggagCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCAc                                                      <  1:799078/100‑1 (MQ=255)
                                          agCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTg                                                    <  1:906329/100‑1 (MQ=255)
                                            cTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTggg                                                  <  2:3267283/100‑1 (MQ=255)
                                              aTGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCg                                                <  1:1019497/100‑1 (MQ=255)
                                              aTGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCg                                                >  2:3061975/1‑100 (MQ=255)
                                                gtcgtcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGat                                              <  1:14595/100‑1 (MQ=255)
                                                  gatctaCCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGatat                                            <  1:460034/95‑1 (MQ=255)
                                                  cgtcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGatat                                            >  2:2346678/1‑100 (MQ=255)
                                                   gtcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATc                                           >  1:2064136/1‑100 (MQ=255)
                                                   gtcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATc                                           <  2:3286444/100‑1 (MQ=255)
                                                   gtcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTCGCTGCTAAATCACTGGGCGATATc                                           >  1:1738384/1‑100 (MQ=255)
                                                    tcAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCa                                          >  2:3858706/1‑100 (MQ=255)
                                                      aaCCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAcc                                        <  1:1311026/100‑1 (MQ=255)
                                                       aCCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGatatagatc                                       >  2:451295/1‑95 (MQ=255)
                                                       aCCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCAccc                                       <  2:3475653/100‑1 (MQ=255)
                                                        cccTACGAGAAGCCAGTAAGGACACGTTGCAGGCAAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCg                                      <  2:2916704/100‑1 (MQ=255)
                                                         ccTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGt                                     >  2:3630192/1‑100 (MQ=255)
                                                         ccTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGt                                     >  1:2600660/1‑100 (MQ=255)
                                                          cTACGAGAAGGCAGTAAGGACACGTTGCAGGCCAGAGAGAAAACTTACCACTAGTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTc                                    >  2:3689044/1‑100 (MQ=255)
                                                          cTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAGAGAGAAAACTTACCACTAATACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTc                                    >  2:3687975/1‑100 (MQ=255)
                                                          cTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTc                                    >  2:885314/1‑100 (MQ=255)
                                                          cTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTc                                    <  1:1959145/100‑1 (MQ=255)
                                                          cTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTc                                    >  2:1603347/1‑100 (MQ=255)
                                                           tACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCt                                   >  1:3225123/1‑100 (MQ=255)
                                                            aCGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTa                                  >  1:464187/1‑100 (MQ=255)
                                                             caagaAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACCACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTAc                                 <  2:1293646/98‑1 (MQ=255)
                                                             cGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTAc                                 >  2:201882/1‑100 (MQ=255)
                                                             cGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTAc                                 >  1:671149/1‑100 (MQ=255)
                                                              gagaAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTAcc                                <  2:1295651/100‑1 (MQ=255)
                                                              gagaAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACGTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTAcc                                <  1:3602814/100‑1 (MQ=255)
                                                                gaAGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCa                              >  1:3780202/1‑100 (MQ=255)
                                                                  aGCCAGTAAGGACACGTTGCATGCCAAATAGAAAACGTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACGCAAg                            >  2:1385473/1‑100 (MQ=255)
                                                                  aGCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAg                            >  1:1092010/1‑100 (MQ=255)
                                                                   nccAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGt                           >  1:1589207/2‑100 (MQ=255)
                                                                   gCCAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGt                           <  1:3877312/100‑1 (MQ=255)
                                                                    ccAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTc                          >  1:3444243/1‑100 (MQ=255)
                                                                    ccAGTAAGGACACGTGGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCGGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTc                          <  2:303542/100‑1 (MQ=255)
                                                                    ccAGCAAGGACACGTTGCAGGCCAAAGAGAAAACTTCCCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTTGGCGATATCACCCGTCTACCCAAGTc                          <  1:3265158/100‑1 (MQ=255)
                                                                     cAGTAAGGACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCa                         >  2:894553/1‑100 (MQ=255)
                                                                           ggACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCaaa                   <  1:29371/100‑1 (MQ=255)
                                                                            gACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTCCAGCCTGCCGCTTGCTGCTAAATCCCTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAg                  >  2:2865400/1‑100 (MQ=255)
                                                                            gACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAg                  >  1:2992479/1‑100 (MQ=255)
                                                                            gACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAg                  >  1:2365683/1‑100 (MQ=255)
                                                                            gACACGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAg                  >  1:1452262/1‑100 (MQ=255)
                                                                              cacGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtt                >  2:1356805/1‑100 (MQ=255)
                                                                               acGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt               <  2:3175369/100‑1 (MQ=255)
                                                                               acGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGttt               <  1:1814064/100‑1 (MQ=255)
                                                                                cGTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGtttt              >  2:2563870/1‑100 (MQ=255)
                                                                                 gTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTg             >  2:2519770/1‑100 (MQ=255)
                                                                                 gTTGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTg             <  1:1768062/100‑1 (MQ=255)
                                                                                  ttGCAGGCCAAAGAGAAAACTTACCACTGCTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc            <  2:493179/100‑1 (MQ=255)
                                                                                  ttGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGc            <  1:1066818/100‑1 (MQ=255)
                                                                                   tGCAGGCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCt           <  1:3282908/100‑1 (MQ=255)
                                                                                       ggCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGaa       >  1:2758232/1‑100 (MQ=255)
                                                                                       ggCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGaa       >  2:613178/1‑100 (MQ=255)
                                                                                       ggCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGaa       >  1:789878/1‑100 (MQ=255)
                                                                                       ggCCAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACACGTCTACCCAAGTCACTCAAAGTTTTGCTCGaa       >  1:1669888/1‑100 (MQ=255)
                                                                                         ccAAAGAGAAAACTTACCAGTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGTTCGaaaa     >  2:3659092/1‑100 (MQ=255)
                                                                                          cAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAAc    >  2:1856947/1‑100 (MQ=255)
                                                                                          cAAAGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAAc    >  1:3432988/1‑100 (MQ=255)
                                                                                           aaaGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAAcc   <  1:478955/100‑1 (MQ=255)
                                                                                            aaGAGAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACct  >  1:2927004/1‑100 (MQ=255)
                                                                                                |                                                                                               
CCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGTTTTGCTCGAAAACCT  >  NC_000913/1335785‑1335976

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