BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,974,008	C→T	A40V (GCG→GTG)	wzxE →	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,974,008	0	C	T	99.2%	429.3 / ‑5.2	130	A40V (GCG→GTG)	wzxE	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (58/71); total (58/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

Only 100 of 136 total aligned reads displayed.

GGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACC  >  NC_000913/3973922‑3974097

ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaaa                                                                                     >  2:4171434/1‑93 (MQ=255)

ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaaa                                                                                     <  2:1235458/93‑1 (MQ=255)

 gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAAt                                                                                    >  2:49930/1‑93 (MQ=255)

 gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAAt                                                                                    <  2:228499/93‑1 (MQ=255)

  ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAAtt                                                                                   >  2:1721069/1‑93 (MQ=255)

    aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTc                                                                                 >  1:430177/1‑93 (MQ=255)

     gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTcc                                                                                >  2:4047214/1‑93 (MQ=255)

     gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTcc                                                                                <  2:2323568/93‑1 (MQ=255)

       acacTGGTCAAGATTGGTGCCGTGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTTGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGc                                                                              <  1:2948883/93‑1 (MQ=255)

         acTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCa                                                                            <  1:3253055/93‑1 (MQ=255)

           tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGt                                                                          <  2:2316770/93‑1 (MQ=255)

            ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCattg                                                                         >  1:876107/1‑90 (MQ=255)

                aaGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAtt                                                                     >  2:116895/1‑93 (MQ=255)

                 aGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTa                                                                    >  2:737379/1‑93 (MQ=255)

                  gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTAc                                                                   >  2:3067399/1‑93 (MQ=255)

                   aTTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTAcc                                                                  <  1:4190199/93‑1 (MQ=255)

                    ttGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCg                                                                 >  2:2914486/1‑93 (MQ=255)

                     tGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGt                                                                <  1:1502180/93‑1 (MQ=255)

                       gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGc                                                              <  1:2343666/93‑1 (MQ=255)

                       gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGc                                                              <  2:1535550/93‑1 (MQ=255)

                       gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGc                                                              <  1:1885950/93‑1 (MQ=255)

                         gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                            >  1:2696207/1‑93 (MQ=255)

                         gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                            <  1:4341183/93‑1 (MQ=255)

                         gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                            <  2:4648786/93‑1 (MQ=255)

                          ccGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCg                                                           <  2:4387572/93‑1 (MQ=255)

                           cGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                          >  2:2464859/1‑93 (MQ=255)

                           cGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                          <  1:497701/93‑1 (MQ=255)

                            gggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTGCCGCAGTTGGTTACCCTGCTCGGc                                                         >  1:3349946/1‑93 (MQ=255)

                            gggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGc                                                         >  1:1087679/1‑93 (MQ=255)

                             ggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCg                                                        >  1:1462861/1‑93 (MQ=255)

                               ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTg                                                      >  1:1873016/1‑93 (MQ=255)

                                 aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                    <  1:2098265/93‑1 (MQ=255)

                                 aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                    <  1:1424436/93‑1 (MQ=255)

                                  cTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                   <  2:3875176/93‑1 (MQ=255)

                                   tGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                  <  2:2252295/93‑1 (MQ=255)

                                   tGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                  >  1:3482265/1‑93 (MQ=255)

                                     gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                <  2:1720573/93‑1 (MQ=255)

                                     gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                <  2:2944907/93‑1 (MQ=255)

                                      tCGGTAAGTTGCTGGCGGTGTCATTTGGTCTGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                               <  2:3213343/93‑1 (MQ=255)

                                       cGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCTGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                              <  1:2893896/93‑1 (MQ=255)

                                        ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                             >  2:286837/1‑93 (MQ=255)

                                        ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTGGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                             <  1:2452874/93‑1 (MQ=255)

                                         gTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                            <  2:2674997/93‑1 (MQ=255)

                                          tAAGTTGCTGGTGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  2:1275559/93‑1 (MQ=255)

                                          tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           >  2:255162/1‑93 (MQ=255)

                                          tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  1:1318501/93‑1 (MQ=255)

                                          tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  1:1110842/93‑1 (MQ=255)

                                          tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  2:509668/93‑1 (MQ=255)

                                           aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                          >  1:497702/1‑93 (MQ=255)

                                           aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                          <  1:4467963/93‑1 (MQ=255)

                                           aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                          >  1:1162157/1‑93 (MQ=255)

                                           aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                          <  2:2224710/93‑1 (MQ=255)

                                           aaGTGGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                          <  2:3900792/93‑1 (MQ=255)

                                            aGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                         >  2:3363670/1‑93 (MQ=255)

                                            aGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                         <  2:865095/93‑1 (MQ=255)

                                            aGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGGAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                         >  2:2120583/1‑93 (MQ=255)

                                             gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                        >  2:4416493/1‑93 (MQ=255)

                                               tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                      <  1:783799/93‑1 (MQ=255)

                                               tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                      >  2:2504951/1‑93 (MQ=255)

                                                 cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                    <  2:2897737/93‑1 (MQ=255)

                                                   ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                  <  1:2691240/93‑1 (MQ=255)

                                                     cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                >  2:2855388/1‑93 (MQ=255)

                                                     cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                <  1:3683545/93‑1 (MQ=255)

                                                      ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                               >  2:1556868/1‑93 (MQ=255)

                                                      ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                               <  2:802327/93‑1 (MQ=255)

                                                        tgtCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                             >  1:328441/1‑93 (MQ=255)

                                                         gtCATGTGGTCCGGCGGGGCTTGGTCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                            <  1:89854/93‑1 (MQ=255)

                                                          tCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                           <  2:427540/93‑1 (MQ=255)

                                                          tCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACCgt                           >  1:4421658/1‑93 (MQ=255)

                                                           cATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                          >  2:3177537/1‑93 (MQ=255)

                                                            aTTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtgt                         <  2:2168489/93‑1 (MQ=255)

                                                             tttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                        >  1:3934765/1‑93 (MQ=255)

                                                             tttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                        >  2:3036358/1‑93 (MQ=255)

                                                              ttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTaa                       <  1:1447203/93‑1 (MQ=255)

                                                              tcggTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTaa                       <  1:3944061/91‑1 (MQ=255)

                                                               tGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTTATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                      <  2:1198219/93‑1 (MQ=255)

                                                               tGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                      <  1:2507527/93‑1 (MQ=255)

                                                               tGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                      <  1:495414/93‑1 (MQ=255)

                                                                 gTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                    >  1:2064520/1‑93 (MQ=255)

                                                                  tCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                   >  1:794724/1‑93 (MQ=255)

                                                                  tCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                   >  2:1659421/1‑93 (MQ=255)

                                                                   ccGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaaa                  <  1:3826572/93‑1 (MQ=255)

                                                                   ccGGCGGGGCTGGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaaa                  <  1:3826573/93‑1 (MQ=255)

                                                                    cggcggGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                 >  2:4048138/1‑93 (MQ=255)

                                                                    cggcggGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                 >  2:4048139/1‑93 (MQ=255)

                                                                    cggcggGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                 <  1:1075210/93‑1 (MQ=255)

                                                                     ggcggGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                <  2:4636842/93‑1 (MQ=255)

                                                                        ggggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt             >  1:356527/1‑93 (MQ=255)

                                                                        ggggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt             >  1:3948408/1‑93 (MQ=255)

                                                                        ggggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt             >  1:235198/1‑93 (MQ=255)

                                                                         gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt            <  2:2130358/93‑1 (MQ=255)

                                                                         gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt            <  1:2758111/93‑1 (MQ=255)

                                                                         gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAGCTTTAACGGTGTAACCAAATAAGtt            >  2:3222936/1‑93 (MQ=255)

                                                                           gCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGAGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTACCCAGATAAGTTga          >  2:3508891/1‑92 (MQ=255)

                                                                             ttGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGccc        <  2:179483/93‑1 (MQ=255)

                                                                               gggCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg      >  2:3930854/1‑93 (MQ=255)

                                                                               gggCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg      <  2:4155883/93‑1 (MQ=255)

                                                                                  cTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGGTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc   >  1:180793/1‑93 (MQ=255)

                                                                                  cTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc   >  2:761419/1‑93 (MQ=255)

                                                                                   tggtggCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGAGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc  <  2:3473130/93‑1 (MQ=255)

GGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACC  >  NC_000913/3973922‑3974097

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