BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,072	T→C	D73G (GAC→GGC)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,072	0	T	C	98.7%	710.2 / ‑1.9	223	D73G (GAC→GGC)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (124/95); total (127/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 232 total aligned reads displayed.

ATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596980‑3597161

atgaAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGt                                                                                           <  1:4106030/93‑1 (MQ=255)

atgaAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGt                                                                                           >  1:1556074/1‑93 (MQ=255)

  gaagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCa                                                                                         <  1:4694574/93‑1 (MQ=255)

     gaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCatat                                                                                      <  2:3584941/93‑1 (MQ=255)

       aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  2:2203923/1‑93 (MQ=255)

         agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttttc                                                                                  >  1:890221/1‑92 (MQ=255)

         agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:3831892/1‑93 (MQ=255)

         agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  <  1:3837115/93‑1 (MQ=255)

          gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 <  1:122828/93‑1 (MQ=255)

          gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:1305750/1‑93 (MQ=255)

            ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCTCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  2:1265156/93‑1 (MQ=255)

            ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  2:3643120/93‑1 (MQ=255)

            ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  2:3666999/93‑1 (MQ=255)

             gTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAt                                                                              <  2:1097345/93‑1 (MQ=255)

              tGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAtt                                                                             <  2:2896220/93‑1 (MQ=255)

               gACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAttt                                                                            <  2:566339/93‑1 (MQ=255)

                 ccAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGc                                                                          >  2:3532635/1‑93 (MQ=255)

                 ccAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGc                                                                          <  2:1952919/93‑1 (MQ=255)

                  cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         <  1:2016549/93‑1 (MQ=255)

                  cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  1:2908433/1‑93 (MQ=255)

                   aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  1:1325627/1‑93 (MQ=255)

                   aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:3805668/1‑93 (MQ=255)

                    atgacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGt                                                                       <  2:4654202/90‑1 (MQ=255)

                     tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      <  2:2678260/93‑1 (MQ=255)

                     tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      <  2:3364793/93‑1 (MQ=255)

                      cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     >  1:3057278/1‑93 (MQ=255)

                      cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     <  1:3187575/93‑1 (MQ=255)

                        caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   <  2:3907624/93‑1 (MQ=255)

                        caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   >  2:402394/1‑93 (MQ=255)

                        caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   >  2:4727104/1‑93 (MQ=255)

                         atatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgtt                                                                  <  1:2443190/93‑1 (MQ=255)

                            tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGct                                                               >  1:3469132/1‑92 (MQ=255)

                            tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGcc                                                               >  1:1784456/1‑93 (MQ=255)

                            tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGcc                                                               >  2:4078372/1‑93 (MQ=255)

                              tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             >  1:2939726/1‑93 (MQ=255)

                               aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            <  2:2000579/93‑1 (MQ=255)

                               aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            >  2:725915/1‑93 (MQ=255)

                                aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           <  1:1701236/93‑1 (MQ=255)

                                     gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      >  1:4230014/1‑93 (MQ=255)

                                     gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      <  2:52840/93‑1 (MQ=255)

                                     gtcgtTAACGACTTTGTTCGACACCGCAACCGCCTGTTTCGGGTCACAGGCACCGCCATATTTTACGATTCGCAGTTTTTTATCTTTAGTgcc                                                      >  1:893393/1‑93 (MQ=255)

                                        gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   >  1:1437755/1‑93 (MQ=255)

                                        gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   >  1:3248975/1‑93 (MQ=255)

                                        gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   >  2:1567031/1‑93 (MQ=255)

                                         ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  >  2:2702584/1‑93 (MQ=255)

                                          tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                 >  2:2334400/1‑93 (MQ=255)

                                          tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                 >  1:64480/1‑93 (MQ=255)

                                          tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                 <  2:794944/93‑1 (MQ=255)

                                            aCGACTTTGTTCGCCACCGCAACCTCCTGTTGCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  2:3275575/93‑1 (MQ=255)

                                            aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  2:4358732/93‑1 (MQ=255)

                                             cGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAg                                              <  1:4002593/93‑1 (MQ=255)

                                               aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTATGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            <  1:1689595/93‑1 (MQ=255)

                                                cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           <  2:413243/93‑1 (MQ=255)

                                                 tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  1:2657779/1‑93 (MQ=255)

                                                 tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  2:2830179/1‑93 (MQ=255)

                                                   tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        <  1:2624314/93‑1 (MQ=255)

                                                   tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        >  1:1087268/1‑93 (MQ=255)

                                                    gTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGGTGat                                       >  2:4548420/1‑93 (MQ=255)

                                                     ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      <  1:1385716/93‑1 (MQ=255)

                                                        gCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   <  2:2775959/93‑1 (MQ=255)

                                                        cccACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   <  1:1964667/92‑1 (MQ=255)

                                                         ccACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCg                                  >  1:4502385/1‑93 (MQ=255)

                                                         ccACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCg                                  >  1:1642543/1‑93 (MQ=255)

                                                         ccACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCg                                  <  1:1880739/93‑1 (MQ=255)

                                                          cACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGc                                 >  2:4524791/1‑93 (MQ=255)

                                                           aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  2:3372970/1‑93 (MQ=255)

                                                            ccGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCaa                               <  2:4666279/93‑1 (MQ=255)

                                                              gCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAcc                             <  2:1449753/93‑1 (MQ=255)

                                                                aaCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGc                           >  2:3017304/1‑93 (MQ=255)

                                                                 aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  2:1672589/1‑93 (MQ=255)

                                                                 aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTGTAGCGTTGATATCCGCAACCGcc                          >  2:2591887/1‑93 (MQ=255)

                                                                  ccGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCt                         <  2:652737/93‑1 (MQ=255)

                                                                    tccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       <  2:3165042/92‑1 (MQ=255)

                                                                    gCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       >  1:116834/1‑93 (MQ=255)

                                                                     ccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGct                      >  1:2399795/1‑93 (MQ=255)

                                                                     ccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGct                      >  1:3430722/1‑93 (MQ=255)

                                                                     ccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGct                      <  2:1877988/93‑1 (MQ=255)

                                                                     ccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGct                      >  1:4046205/1‑93 (MQ=255)

                                                                     ccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGct                      >  1:4179366/1‑93 (MQ=255)

                                                                      cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  2:807698/1‑93 (MQ=255)

                                                                      cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  2:1776179/93‑1 (MQ=255)

                                                                       tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    <  2:903587/93‑1 (MQ=255)

                                                                       tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    <  1:3764878/93‑1 (MQ=255)

                                                                         tttCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCACCCGCCTGCTCTgc                  <  2:3207847/93‑1 (MQ=255)

                                                                            cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGcc               >  2:2271019/1‑93 (MQ=255)

                                                                            cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGcc               >  1:1209160/1‑93 (MQ=255)

                                                                            cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGcc               >  2:4448854/1‑93 (MQ=255)

                                                                              ggTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCgg             >  1:3567363/1‑93 (MQ=255)

                                                                              ggTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCgg             >  1:2104217/1‑93 (MQ=255)

                                                                                tCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTa           >  1:961668/1‑93 (MQ=255)

                                                                                 cacaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaa          >  1:478448/1‑93 (MQ=255)

                                                                                  acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  2:757624/1‑93 (MQ=255)

                                                                                     ggCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTc      <  2:2488880/93‑1 (MQ=255)

                                                                                     gacATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTc      >  2:2864015/3‑93 (MQ=255)

                                                                                      gCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     >  1:4408505/1‑93 (MQ=255)

                                                                                      gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  1:2188534/93‑1 (MQ=255)

                                                                                       cATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCt    <  1:4210703/93‑1 (MQ=255)

                                                                                        aTCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTg   <  2:2262050/93‑1 (MQ=255)

                                                                                        aTCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTg   >  1:1290562/1‑93 (MQ=255)

                                                                                         tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:4438839/1‑93 (MQ=255)

ATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596980‑3597161

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