Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,078,245 C→G G470G (GGG→GGC tktA ← transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,2450CG98.5% 426.2 / NA 140G470G (GGG→GGCtktAtransketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (62/73);  total (62/78)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

GCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078156‑3078334
                                                                                          |                                                                                         
gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCg                                                                                        >  2:2518575/1‑93 (MQ=255)
gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGTCAACCCGGCTGGGGAGCCGGCCCg                                                                                        >  2:3757860/1‑93 (MQ=255)
  ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                                                      >  2:264884/1‑93 (MQ=255)
       cAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGc                                                                                 <  2:243701/93‑1 (MQ=255)
       cAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGc                                                                                 >  1:1849149/1‑93 (MQ=255)
        caCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                <  1:889088/92‑1 (MQ=255)
        caCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                <  1:2286878/92‑1 (MQ=255)
        aaCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                <  1:889207/93‑1 (MQ=255)
         aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCcccc                                                                               >  2:1968508/1‑93 (MQ=255)
         aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCACAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGACCGTCTTCGccc                                                                               <  2:3822575/93‑1 (MQ=255)
          ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCa                                                                              <  1:409463/93‑1 (MQ=255)
          ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCa                                                                              <  2:3320703/93‑1 (MQ=255)
            tGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGa                                                                            <  1:1158246/93‑1 (MQ=255)
            tGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGa                                                                            >  2:3061956/1‑93 (MQ=255)
             ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                                           >  1:4203333/1‑93 (MQ=255)
              gTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                                                                          <  1:4390778/93‑1 (MQ=255)
                cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGTTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                                        <  2:3019452/93‑1 (MQ=255)
                cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                                        <  2:4715119/93‑1 (MQ=255)
                 acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       >  2:3327148/1‑93 (MQ=255)
                 acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       >  2:4023044/1‑93 (MQ=255)
                 acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       <  1:2737995/93‑1 (MQ=255)
                  cacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATg                                                                      <  1:1461284/93‑1 (MQ=255)
                    cGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACGGGCTGGATAGTCGGGCCGTCTTCGCCCAGACCGATGGa                                                                    <  2:3693912/93‑1 (MQ=255)
                    cGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCATCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGa                                                                    <  1:2391761/93‑1 (MQ=255)
                    cGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGa                                                                    <  2:216593/93‑1 (MQ=255)
                     ggACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGCGTCGGGCCGTCTTCGCCCAGCCCGGTggcg                                                                   >  2:855803/1‑91 (MQ=255)
                     ggACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAg                                                                   <  1:2198170/93‑1 (MQ=255)
                      gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGt                                                                  <  2:4614635/93‑1 (MQ=255)
                        cGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCg                                                                >  2:4692570/1‑93 (MQ=255)
                         tccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGGCGGGCCGTCTTCGCCCAGACCGGTGGAGTCGt                                                               <  2:2386414/92‑1 (MQ=255)
                         gCCATGTAGTCATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               <  2:955880/93‑1 (MQ=255)
                         gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCgg                                                               >  2:4119334/1‑92 (MQ=255)
                         gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               >  2:3443269/1‑93 (MQ=255)
                         gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               <  1:186879/93‑1 (MQ=255)
                         gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               <  1:172843/93‑1 (MQ=255)
                          ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                                              >  2:1982127/1‑93 (MQ=255)
                          ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                                              >  1:3114835/1‑93 (MQ=255)
                          ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                                              >  1:2557411/1‑93 (MQ=255)
                           cATGTAGACATGTTCGGGGTAACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                                             <  2:1835746/93‑1 (MQ=255)
                            aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            <  1:1083608/93‑1 (MQ=255)
                             tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgggg                                                           >  2:1641713/1‑92 (MQ=255)
                             tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgt                                                           >  1:3323173/1‑93 (MQ=255)
                             tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgt                                                           <  1:4379677/93‑1 (MQ=255)
                               tAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                                         >  2:956167/1‑93 (MQ=255)
                               tAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                                         <  1:3618267/93‑1 (MQ=255)
                               gaGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                                         <  1:3618266/92‑1 (MQ=255)
                                aGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGGCGTCTTCGCCCAGACCGATGGAGCCGGGGGTGTa                                                        >  2:1000115/1‑93 (MQ=255)
                                   cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                                                     <  1:702381/93‑1 (MQ=255)
                                   cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                                                     >  2:410157/1‑93 (MQ=255)
                                   cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                                                     <  2:1494991/93‑1 (MQ=255)
                                       ttCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATc                                                 >  1:908298/1‑93 (MQ=255)
                                        tCGTGGTTACGGGCAGAGAAGCGACCTGCTCCACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCa                                                <  1:1144762/93‑1 (MQ=255)
                                         cGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAc                                               <  1:3416246/93‑1 (MQ=255)
                                           gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  1:4554234/93‑1 (MQ=255)
                                           gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  2:3231510/93‑1 (MQ=255)
                                           gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  2:3231387/93‑1 (MQ=255)
                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  2:688489/93‑1 (MQ=255)
                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  1:736930/93‑1 (MQ=255)
                                              ttACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGGGTAAACCATTACCTGAc                                          >  1:2929202/1‑93 (MQ=255)
                                              gtACGCGCAGAGCAGCGCCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  2:1034123/92‑1 (MQ=255)
                                               tACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAcc                                         >  2:1225367/1‑92 (MQ=255)
                                                 cgcgCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAGGCt                                       >  2:1296693/1‑93 (MQ=255)
                                                   cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGt                                     <  2:4345745/93‑1 (MQ=255)
                                                   cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGt                                     <  2:3547800/93‑1 (MQ=255)
                                                   cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGt                                     <  1:4084975/93‑1 (MQ=255)
                                                    gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGtt                                    <  1:1334086/93‑1 (MQ=255)
                                                    gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAACTGACGCTGtt                                    >  1:3862276/1‑93 (MQ=255)
                                                     cAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGttt                                   <  1:3376525/93‑1 (MQ=255)
                                                      agagAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtc                                  <  1:384097/93‑1 (MQ=255)
                                                      agagAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtc                                  <  1:2405559/93‑1 (MQ=255)
                                                      agagAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtc                                  >  2:4596657/1‑93 (MQ=255)
                                                       gagaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca                                 >  1:3334919/1‑93 (MQ=255)
                                                         gaAGCGACCTGCTCATCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcata                               <  2:4217138/93‑2 (MQ=255)
                                                           aGCGGCCTGCTCATCCGGCTGGTGAGCCGGTCCGTCTTCGCCCAGACCGATGGAGTCGGGAGTGTAATCCATCTCCTGACGCTGTTTCATCAg                             <  2:4503318/93‑1 (MQ=255)
                                                           aGCGACCTGCTCAACCGGCTGGTGAGTCGGTCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAg                             <  1:1796323/93‑1 (MQ=255)
                                                            gCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc                            <  2:3845024/93‑1 (MQ=255)
                                                              gACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc                          >  2:3973104/1‑93 (MQ=255)
                                                               aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCCCCTGACGCTGTTTCATCAGCGCa                         <  1:3489072/93‑1 (MQ=255)
                                                                 cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGc                       <  1:4285927/93‑1 (MQ=255)
                                                                   gCTCACCCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACTGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCa                     <  2:23895/93‑1 (MQ=255)
                                                                   gCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGGGGGTGTAAACCATCACCTGACGCTGGTTGAGCAGCGCAGGCa                     >  2:385243/1‑93 (MQ=255)
                                                                    cTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAt                    >  2:3538027/1‑93 (MQ=255)
                                                                    cTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAt                    <  2:766211/93‑1 (MQ=255)
                                                                     tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   >  1:1336576/1‑93 (MQ=255)
                                                                     tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   <  1:4672493/93‑1 (MQ=255)
                                                                      cAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtac                  <  2:1512818/93‑1 (MQ=255)
                                                                      cAACCGGCTGGGGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtac                  <  1:1473095/93‑1 (MQ=255)
                                                                          ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               >  2:2738332/1‑93 (MQ=255)
                                                                             gCTTGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCGGTTTCATCAGCGCAGCCATACGTACgg            <  1:1511580/93‑1 (MQ=255)
                                                                             gCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACgg            <  1:2760099/93‑1 (MQ=255)
                                                                             gCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACgg            >  2:2196071/1‑93 (MQ=255)
                                                                               tGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCg          >  2:692339/1‑93 (MQ=255)
                                                                               tGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCg          >  1:3489571/1‑93 (MQ=255)
                                                                               tGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGGCGTGGGTGTACACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCg          >  1:52386/1‑93 (MQ=255)
                                                                                ggTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGt         >  2:3461800/1‑93 (MQ=255)
                                                                                  tGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTa       >  1:458758/1‑93 (MQ=255)
                                                                                       cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  >  1:270541/1‑93 (MQ=255)
                                                                                       cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  1:1507063/93‑1 (MQ=255)
                                                                                       cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  1:921293/93‑1 (MQ=255)
                                                                                       cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  1:1015599/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078156‑3078334

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