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.6%	782.2 / NA	245	G470G (GGG→GGC)	tktA	transketolase 1 isozyme
Reads supporting (aligned to +/- strand): ref base C (0/0); major base G (134/110); minor base T (0/1); total (134/111)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 259 total aligned reads displayed.

CGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078155‑3078334

 gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCg                                                                                        >  2:5939802/1‑93 (MQ=255)

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

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

  cgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGt                                                                                       <  2:1027819/93‑1 (MQ=255)

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

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

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

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

     aTTCAAGCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCtt                                                                                    >  1:5555457/1‑93 (MQ=255)

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

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

      ttCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTc                                                                                   >  1:1702245/1‑93 (MQ=255)

       tCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCttag                                                                                  >  2:2322703/1‑91 (MQ=255)

       tCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCg                                                                                  >  1:1018254/1‑93 (MQ=255)

        cAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGc                                                                                 <  2:6277948/93‑1 (MQ=255)

        cAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGc                                                                                 >  1:3062329/1‑93 (MQ=255)

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

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

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

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

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

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

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

                    acGGACGCCATGTAGACATGTTCGGGGGTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATgg                                                                     >  2:4755319/1‑93 (MQ=255)

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

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

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

                        aCGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTc                                                                 <  1:2804969/93‑1 (MQ=255)

                          gCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                                                               <  2:2700922/93‑1 (MQ=255)

                            cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                                             >  1:6534820/1‑93 (MQ=255)

                            cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                                             >  1:4948971/1‑93 (MQ=255)

                             aTGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTggg                                                            >  2:3250108/1‑93 (MQ=255)

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

                              tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgt                                                           >  2:1888947/1‑93 (MQ=255)

                               gTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtg                                                          >  1:4172853/1‑93 (MQ=255)

                                   aCATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTaaa                                                      >  1:163357/1‑93 (MQ=255)

                                   aCATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTaaa                                                      <  2:2014273/93‑1 (MQ=255)

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

                                     aTGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAcc                                                    >  2:5783086/1‑93 (MQ=255)

                                      tGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCa                                                   >  1:3309541/1‑93 (MQ=255)

                                      tGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCa                                                   <  2:1136876/93‑1 (MQ=255)

                                      tGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCa                                                   >  1:2482933/1‑93 (MQ=255)

                                      tGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCa                                                   <  2:5088422/93‑1 (MQ=255)

                                       gTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCAt                                                  <  1:6418083/93‑1 (MQ=255)

                                           ggggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc                                              <  2:4886035/93‑1 (MQ=255)

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

                                             ggTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAACTg                                            >  2:6457906/1‑93 (MQ=255)

                                              gTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGa                                           >  1:5882395/1‑93 (MQ=255)

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

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

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

                                                 aCGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGc                                        >  2:1553813/1‑93 (MQ=255)

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

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

                                                  cgcgCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCCGACCGATGGAGTCGTGGGTGTAAACCATCACCTGActct                                       >  2:55959/1‑93 (MQ=255)

                                                    cgcAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGt                                     <  1:1862585/93‑1 (MQ=255)

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

                                                      cAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGttt                                   <  1:2757628/93‑1 (MQ=255)

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

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

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

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

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

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

                                                          gaAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc                               <  1:1562846/93‑1 (MQ=255)

                                                            aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgt                                                           >  1:6701371/1‑63 (MQ=255)

                                                            aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATGAg                             >  1:2023187/1‑93 (MQ=255)

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

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

                                                            aGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAg                             <  1:5721943/93‑1 (MQ=255)

                                                             gCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc                            >  2:3090490/1‑93 (MQ=255)

                                                             gCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc                            >  2:4851809/1‑93 (MQ=255)

                                                               gACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc                          <  1:749613/93‑1 (MQ=255)

                                                                aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCa                         >  2:711856/1‑93 (MQ=255)

                                                                aCCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCa                         >  2:2605725/1‑93 (MQ=255)

                                                                  cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGc                       >  1:1942771/1‑93 (MQ=255)

                                                                  cTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGc                       >  1:133162/1‑93 (MQ=255)

                                                                    gCTCAACCGGCTGGTGAGTTGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCACCGCAGCCa                     >  1:469124/1‑93 (MQ=255)

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

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

                                                                       cAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtac                  <  2:5243088/93‑1 (MQ=255)

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

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

                                                                         aCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgt                <  2:544102/93‑1 (MQ=255)

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

                                                                          ccGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgta               >  2:282988/1‑93 (MQ=255)

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

                                                                             gCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtacgc            >  2:4681319/1‑92 (MQ=255)

                                                                              cTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCAtacgtacgcg           >  1:5421982/1‑91 (MQ=255)

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

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

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

                                                                                 gTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGtt        <  2:1713682/93‑1 (MQ=255)

                                                                                 gTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGtt        <  2:1332496/93‑1 (MQ=255)

                                                                                  tGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTa       >  1:3868302/1‑93 (MQ=255)

                                                                                  tGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTa       >  1:1186352/1‑93 (MQ=255)

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

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

                                                                                     gTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGt    <  1:6264927/93‑1 (MQ=255)

                                                                                      tCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTg   >  2:5880892/1‑93 (MQ=255)

CGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGCAGCCATACGTACGGCGTTACGTGC  >  NC_000913/3078155‑3078334

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