BRESEQ :: Evidence

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_000913	3,078,245	0	C	G	99.7%	1066.8 / NA	321	G470G (GGG→GGC)	tktA	transketolase 1 isozyme
Reads supporting (aligned to +/- strand): ref base C (1/0); new base G (164/156); total (165/156)
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.98e-01

Only 100 of 359 total aligned reads displayed.

CGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078155‑3078334

 gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGGGAGTCGGGCCg                                                                                        >  1:586594/1‑93 (MQ=255)

  cgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGt                                                                                       >  2:3282307/1‑93 (MQ=255)

   ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                                                      <  1:856679/93‑1 (MQ=255)

    gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                                                     <  2:4300954/93‑1 (MQ=255)

    gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGCCt                                                                                     >  1:1248293/1‑93 (MQ=255)

     aTTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCtt                                                                                    <  1:7601116/93‑1 (MQ=255)

     aTTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCtt                                                                                    >  1:289385/1‑93 (MQ=255)

         aaCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGcc                                                                                >  2:8591660/1‑93 (MQ=255)

          aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                                               <  2:1892556/93‑1 (MQ=255)

          aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                                               <  1:7483222/93‑1 (MQ=255)

          aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                                               >  2:459143/1‑93 (MQ=255)

          aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                                               >  1:7462722/1‑93 (MQ=255)

           ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCa                                                                              <  1:5527796/93‑1 (MQ=255)

           ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCa                                                                              >  2:3654391/1‑93 (MQ=255)

              ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                                           <  2:7208286/93‑1 (MQ=255)

              ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                                           >  1:3031359/1‑93 (MQ=255)

              ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                                           <  1:121292/93‑1 (MQ=255)

              ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                                           >  2:5735283/1‑93 (MQ=255)

               gTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAcc                                                                          >  2:8480343/1‑93 (MQ=255)

                tCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGGCGGGCCGTCTTCGCCCAGACCg                                                                         >  2:2297041/1‑93 (MQ=255)

                 cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                                        <  1:147515/93‑1 (MQ=255)

                  acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                                                       >  1:3311681/1‑93 (MQ=255)

                   cacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATg                                                                      <  1:4078960/93‑1 (MQ=255)

                    acGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATgg                                                                     <  1:4810555/93‑1 (MQ=255)

                    acGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATgg                                                                     <  2:7127631/93‑1 (MQ=255)

                     cGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGa                                                                    <  1:437456/93‑1 (MQ=255)

                       gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGt                                                                  >  1:95502/1‑93 (MQ=255)

                             aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCAACCTGCTCATCAGGCTGGTGAGTCGGTCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            <  2:8025262/93‑1 (MQ=255)

                              tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgggg                                                           >  2:3550369/1‑92 (MQ=255)

                                 aGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTa                                                        <  1:6908260/93‑1 (MQ=255)

                                    cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                                                     <  2:6988837/93‑1 (MQ=255)

                                    cATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGCGGGTGTAAAc                                                     >  2:3548592/1‑93 (MQ=255)

                                     aTGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAcc                                                    <  2:339136/93‑1 (MQ=255)

                                     aTGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAcc                                                    <  1:826150/93‑1 (MQ=255)

                                     aTGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCCTCTTCGCACAGACCGATGGAGTCGTGGGTGTAAAcc                                                    >  2:6660194/1‑93 (MQ=255)

                                       gTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCAt                                                  >  1:6141373/1‑93 (MQ=255)

                                        ttCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATc                                                 >  1:4941567/1‑93 (MQ=255)

                                        ttCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATc                                                 <  1:3274768/93‑1 (MQ=255)

                                        ttCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATc                                                 >  1:7895853/1‑93 (MQ=255)

                                         tCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCa                                                <  1:7184152/93‑1 (MQ=255)

                                           ggggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc                                              >  2:3894782/1‑93 (MQ=255)

                                           ggggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGGCGTGGGTGTAAACCATCAcc                                              >  1:6966410/1‑93 (MQ=255)

                                            gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             >  1:2046670/1‑93 (MQ=255)

                                            gggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                                             <  1:4531956/93‑1 (MQ=255)

                                             ggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTg                                            <  1:3721651/93‑1 (MQ=255)

                                             ggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTg                                            <  1:340569/93‑1 (MQ=255)

                                             ggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTg                                            <  2:3831986/93‑1 (MQ=255)

                                             ggTTACACGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACGGATGGAGTCGTGGGTGGAAACCATccccgg                                            >  2:6262249/1‑91 (MQ=255)

                                               ttACGCGCGGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc                                          <  2:2252178/93‑1 (MQ=255)

                                                tACGCGCAGAGAAGCGCCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACg                                         <  2:1206723/93‑1 (MQ=255)

                                                tACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACg                                         >  2:1471991/1‑93 (MQ=255)

                                                tACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACg                                         >  1:353625/1‑93 (MQ=255)

                                                 aCGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGc                                        <  2:1288580/93‑1 (MQ=255)

                                                 aCGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGc                                        <  2:7707409/93‑1 (MQ=255)

                                                     gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGtt                                    >  1:257365/1‑93 (MQ=255)

                                                     gcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGtt                                    <  2:6215988/93‑1 (MQ=255)

                                                       agagAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtc                                  >  2:1595888/1‑93 (MQ=255)

                                                        gagaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca                                 <  2:5047485/93‑1 (MQ=255)

                                                        gagaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca                                 >  1:8798813/1‑93 (MQ=255)

                                                         agaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGGCCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcat                                <  2:6559710/93‑1 (MQ=255)

                                                         agaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcat                                <  1:8763051/93‑1 (MQ=255)

                                                           aaGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtaatca                              <  1:2757974/93‑1 (MQ=255)

                                                            aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAg                             >  2:4391035/1‑93 (MQ=255)

                                                            aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAg                             >  1:8347277/1‑93 (MQ=255)

                                                              cGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcg                           <  1:6240336/93‑1 (MQ=255)

                                                               gACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc                          >  2:293660/1‑93 (MQ=255)

                                                               gACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc                          >  2:2679405/1‑93 (MQ=255)

                                                                aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCCGCGCa                         >  1:1655490/1‑93 (MQ=255)

                                                                 ccTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAg                        >  1:5235268/1‑93 (MQ=255)

                                                                  cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGACGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGc                       >  1:5098127/1‑93 (MQ=255)

                                                                   tGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGcc                      <  2:4855458/93‑1 (MQ=255)

                                                                    gCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCa                     <  2:5513766/93‑1 (MQ=255)

                                                                    gCTCAACCGGCTGGTGAGTCGGCCCTTCTTCGCCCAGACCGATGGAGTCGGGGGTGTAAACCATGACCTGACGCTGGTTGATCAGCGCAGCCa                     >  2:1055247/1‑93 (MQ=255)

                                                                      tCAACCGTCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   <  2:8747229/93‑1 (MQ=255)

                                                                      tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   >  2:7400673/1‑93 (MQ=255)

                                                                      tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   <  2:1749738/93‑1 (MQ=255)

                                                                      tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   <  2:8323498/93‑1 (MQ=255)

                                                                      tCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAta                   >  1:5743102/1‑93 (MQ=255)

                                                                       cAACCGGCTGGTGAGGCGGACCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtac                  <  1:6367471/93‑1 (MQ=255)

                                                                        aaCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacg                 <  2:4518578/93‑1 (MQ=255)

                                                                        aaCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacg                 <  1:8099121/93‑1 (MQ=255)

                                                                        aaCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacg                 >  2:4155944/1‑93 (MQ=255)

                                                                         aCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgt                >  1:2261344/1‑93 (MQ=255)

                                                                          ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               <  1:5642274/93‑1 (MQ=255)

                                                                          ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               <  1:3202638/93‑1 (MQ=255)

                                                                          ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               <  2:7868991/93‑1 (MQ=255)

                                                                           cGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtac              >  2:2400053/1‑93 (MQ=255)

                                                                            ggCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtacg             <  1:1851467/93‑1 (MQ=255)

                                                                            ggCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtacg             <  2:6816602/93‑1 (MQ=255)

                                                                            ggCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtacg             >  2:2093054/1‑93 (MQ=255)

                                                                              cTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGc           >  1:6004330/1‑93 (MQ=255)

                                                                                ggTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGt         >  1:3931274/1‑93 (MQ=255)

                                                                                ggTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGt         <  2:7583584/93‑1 (MQ=255)

                                                                                 gggAGTCGGGCCGGCCTCCGCCAAGCCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGtt        <  1:410061/91‑1 (MQ=255)

                                                                                   gAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTAc      >  2:1403117/1‑93 (MQ=255)

                                                                                   gAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTAc      <  1:6758154/93‑1 (MQ=255)

                                                                                   gAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTAc      <  2:4815235/93‑1 (MQ=255)

                                                                                     gTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGt    >  2:6725623/1‑93 (MQ=255)

                                                                                      tCGGGCCGTCTGCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTg   <  2:5218553/93‑1 (MQ=255)

                                                                                       cGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGc  <  2:474006/93‑1 (MQ=255)

CGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078155‑3078334

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