BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,920,902	T→C	Q31R (CAG→CGG)	b1838 ←	protein phosphatase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,920,902	0	T	C	99.6%	730.5 / ‑5.1	228	Q31R (CAG→CGG)	b1838	protein phosphatase 1
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (117/110); total (117/111)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 232 total aligned reads displayed.

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

 tttGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCt                                                                                        <  1:4313016/93‑1 (MQ=255)

  ttGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTc                                                                                       >  1:4697092/1‑93 (MQ=255)

       gCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGAc                                                                                  <  2:3988329/93‑1 (MQ=255)

        cccACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACa                                                                                 <  1:3296530/93‑1 (MQ=255)

          cACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAAc                                                                               <  2:876904/93‑1 (MQ=255)

          cACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAAc                                                                               >  2:3135779/1‑93 (MQ=255)

          cACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCCCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAAc                                                                               >  1:1769016/1‑93 (MQ=255)

           aCGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAAcc                                                                              >  2:2525116/1‑93 (MQ=255)

            cgatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCg                                                                             <  2:5985845/93‑1 (MQ=255)

             gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgt                                                                            >  1:5416085/1‑93 (MQ=255)

             gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgt                                                                            >  1:3413498/1‑93 (MQ=255)

             gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgg                                                                            >  2:2885171/1‑92 (MQ=255)

              atcgatAACGTCTCCCACTGTGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtg                                                                           >  1:969537/1‑93 (MQ=255)

              atcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtg                                                                           <  2:1808787/93‑1 (MQ=255)

              atcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtg                                                                           <  2:3404741/93‑1 (MQ=255)

              atcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtg                                                                           >  1:1722840/1‑93 (MQ=255)

               tcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtgt                                                                          <  1:1895192/93‑1 (MQ=255)

                cgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGta                                                                         <  2:3962420/93‑1 (MQ=255)

                cgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGta                                                                         >  2:1400462/1‑93 (MQ=255)

                cgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGta                                                                         <  1:3566380/93‑1 (MQ=255)

                cgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGta                                                                         <  2:3872326/93‑1 (MQ=255)

                 gatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGtat                                                                        >  1:3889895/1‑93 (MQ=255)

                  atAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGtata                                                                       >  2:6173634/1‑93 (MQ=255)

                  atAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGtata                                                                       >  2:3762126/1‑93 (MQ=255)

                   tAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTatat                                                                      >  1:2424353/1‑93 (MQ=255)

                    aaCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGAGAACCGTGTATATc                                                                     >  1:821454/1‑93 (MQ=255)

                    aaCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATc                                                                     >  1:1568086/1‑93 (MQ=255)

                     aCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCg                                                                    >  1:3476266/1‑93 (MQ=255)

                     aCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCg                                                                    >  1:5011253/1‑93 (MQ=255)

                        tctcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCa                                                                 <  2:4526295/93‑1 (MQ=255)

                        tctcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCa                                                                 >  2:3442199/1‑93 (MQ=255)

                         ctcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAg                                                                >  2:650997/1‑93 (MQ=255)

                            ccACTGAGATAAGTAAATCTCGCCACGTATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGaaa                                                             <  2:2907627/93‑1 (MQ=255)

                            ccACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGaaa                                                             <  2:3240614/93‑1 (MQ=255)

                               cTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGcc                                                          <  2:1432752/93‑1 (MQ=255)

                               cTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGcc                                                          >  1:1424082/1‑93 (MQ=255)

                                tGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCa                                                         >  2:3273539/1‑93 (MQ=255)

                                 gagaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAg                                                        >  1:117813/1‑93 (MQ=255)

                                  agaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGTAAGCCAGa                                                       <  1:6564848/93‑1 (MQ=255)

                                  agaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGa                                                       <  1:4847074/93‑1 (MQ=255)

                                  agaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGa                                                       >  2:4991144/1‑93 (MQ=255)

                                  agaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGGGTATATCGCCAGAAAGCCAGa                                                       >  1:5066485/1‑93 (MQ=255)

                                   gaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGat                                                      <  1:5259237/93‑1 (MQ=255)

                                   gaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGat                                                      <  2:1268651/93‑1 (MQ=255)

                                    aTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGata                                                     <  1:2078651/93‑1 (MQ=255)

                                    aTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGata                                                     <  1:6241099/93‑1 (MQ=255)

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    <  2:4764020/93‑1 (MQ=255)

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    <  1:5674725/93‑1 (MQ=255)

                                       aGTAAATCTCGCCGCGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                                  <  2:1335252/93‑1 (MQ=255)

                                       aGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                                  >  1:3002110/1‑93 (MQ=255)

                                        gTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTc                                                 >  1:1539560/1‑93 (MQ=255)

                                        gTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTc                                                 <  2:4193195/93‑1 (MQ=255)

                                        gTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTc                                                 >  1:2717214/1‑93 (MQ=255)

                                            aTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCa                                             <  2:724826/93‑1 (MQ=255)

                                            aTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGGGTATATCGCCAGAAAGCCAGATATGTCGcct                                             >  1:4967606/1‑92 (MQ=255)

                                             tctcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCCAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAt                                            <  2:4697456/93‑1 (MQ=255)

                                              ctcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAtt                                           <  1:1337112/93‑1 (MQ=255)

                                              ctcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAtt                                           <  2:2302520/93‑1 (MQ=255)

                                               tcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtg                                          >  2:585112/1‑93 (MQ=255)

                                               tcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATAGCGCCAGAAAGCCAGATATGTCGCCATtg                                          >  2:6450355/1‑93 (MQ=255)

                                                cGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtga                                         <  2:1471283/93‑1 (MQ=255)

                                                 gCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgat                                        >  1:3236394/1‑93 (MQ=255)

                                                   cACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgatga                                      >  2:5391896/1‑93 (MQ=255)

                                                   cACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgatga                                      <  2:2159843/93‑1 (MQ=255)

                                                     cGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAcc                                    <  1:4028810/93‑1 (MQ=255)

                                                      ggATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAccc                                   >  2:5783747/1‑93 (MQ=255)

                                                        aTCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGc                                 >  1:2463900/1‑93 (MQ=255)

                                                         tCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCa                                >  1:2658750/1‑93 (MQ=255)

                                                          cAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCaa                               <  1:4253154/93‑1 (MQ=255)

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              >  1:4149680/1‑93 (MQ=255)

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              >  2:5107862/1‑93 (MQ=255)

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              >  2:1168558/1‑93 (MQ=255)

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              >  1:4593995/1‑93 (MQ=255)

                                                            aaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAtt                             >  2:2742245/1‑93 (MQ=255)

                                                            aaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAtt                             >  1:917578/1‑93 (MQ=255)

                                                             aaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtc                            >  1:685559/1‑93 (MQ=255)

                                                              aCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtct                           >  2:3521853/1‑93 (MQ=255)

                                                              aCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtct                           >  2:112132/1‑93 (MQ=255)

                                                                gACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTctct                         <  1:6222261/93‑1 (MQ=255)

                                                                 aCAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTg                        >  1:5766958/1‑93 (MQ=255)

                                                                  cAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGa                       >  2:1080656/1‑93 (MQ=255)

                                                                   aaTGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGGCCCGCAATTCTCTGAt                      <  1:437920/93‑1 (MQ=255)

                                                                   aaTGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGAt                      <  1:2537165/93‑1 (MQ=255)

                                                                   aaTGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGAt                      >  1:1834965/1‑93 (MQ=255)

                                                                   aaTGCCATAATGTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGAt                      <  1:2226032/93‑1 (MQ=255)

                                                                     tGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaa                    <  1:5974360/93‑1 (MQ=255)

                                                                       ccATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAAc                  <  2:5237955/93‑1 (MQ=255)

                                                                        cATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACg                 >  2:4028026/1‑93 (MQ=255)

                                                                        cATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACg                 >  2:6355085/1‑93 (MQ=255)

                                                                         aTAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACgg                <  1:1610289/93‑1 (MQ=255)

                                                                         aTAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACgg                <  2:5579120/93‑1 (MQ=255)

                                                                          tAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACggg               <  1:52978/93‑1 (MQ=255)

                                                                            aTTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgcg             >  2:1657020/1‑93 (MQ=255)

                                                                             tttGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgcgc            >  1:6118300/1‑93 (MQ=255)

                                                                              ttGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCa           <  2:3901573/93‑1 (MQ=255)

                                                                                 cggcgCAACCGCTCAAGACAGCCGTGTATATCGCCAGAAAGCCAGATATGTCGCCGTTGATGGCCCGCAATTCTCTGATAAACGGGCGCAGGc        <  1:3188067/93‑1 (MQ=255)

                                                                                  ggcgCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCt       <  2:3717553/93‑1 (MQ=255)

                                                                                    cgcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGt     >  1:4011493/1‑93 (MQ=255)

                                                                                       aaCCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTtc  >  2:1699623/1‑93 (MQ=255)

                                                                                       aaCCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTtc  >  1:1079537/1‑93 (MQ=255)

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

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