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%	860.1 / NA	269	Q31R (CAG→CGG)	b1838	protein phosphatase 1
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (131/137); total (132/137)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGc                                                                                         >  1:4300430/1‑93 (MQ=255)

cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGc                                                                                         >  2:2983551/1‑93 (MQ=255)

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

   tGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCa                                                                                      <  1:586583/93‑1 (MQ=255)

       gCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGAc                                                                                  >  1:120736/1‑93 (MQ=255)

        cccACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACa                                                                                 >  2:5200775/1‑93 (MQ=255)

         ccACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACaa                                                                                >  2:2932553/1‑93 (MQ=255)

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

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

             gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgt                                                                            <  2:4777204/93‑1 (MQ=255)

             gatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgt                                                                            <  2:5746166/93‑1 (MQ=255)

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

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

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

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

               tcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtgt                                                                          >  1:3712210/1‑93 (MQ=255)

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

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

                  atAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGtata                                                                       <  2:3327560/93‑1 (MQ=255)

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

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

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

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

                    aaCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCACCCGCTCAAGACAACCGTGTATATc                                                                     <  1:6237808/93‑1 (MQ=255)

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

                          tcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGa                                                               <  2:3105637/93‑1 (MQ=255)

                          tcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGa                                                               >  1:3005020/1‑93 (MQ=255)

                          tcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGa                                                               >  1:6030904/1‑93 (MQ=255)

                            ccACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGaaa                                                             >  2:1167974/1‑93 (MQ=255)

                            ccACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGaaa                                                             >  2:2230438/1‑93 (MQ=255)

                            ccACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGaaa                                                             >  1:4541306/1‑93 (MQ=255)

                             cACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAg                                                            <  1:327076/93‑1 (MQ=255)

                               cTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGACAGca                                                          >  2:5625444/1‑92 (MQ=255)

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

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

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

                                tGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCa                                                         <  1:3667146/93‑1 (MQ=255)

                                 gcgaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAg                                                        >  2:3295886/3‑93 (MQ=255)

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

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

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

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

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    >  2:2862797/1‑93 (MQ=255)

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    >  1:5554076/1‑93 (MQ=255)

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

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    >  1:5435783/1‑93 (MQ=255)

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

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

                                         tAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCg                                                <  1:2215111/93‑1 (MQ=255)

                                          aaaTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGc                                               <  1:1694662/93‑1 (MQ=255)

                                            aTCTCGCCACGTAACAAAACGACAATGCCATTATTCGCGGCGCAACCGCCCAAGACATCCGTGTATACCGCCATAAAGCCAGATATGTCGCCa                                             <  2:6145891/93‑1 (MQ=255)

                                             tctcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAt                                            >  1:3571840/1‑93 (MQ=255)

                                              ctcGCCACGGATCAAAACGACAATGCCATAATCTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAtt                                           <  1:5451072/93‑1 (MQ=255)

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

                                               tcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtg                                          <  1:1059198/93‑1 (MQ=255)

                                                cGCCACGGATCAAAACGACAATGGCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtga                                         <  1:6123609/93‑1 (MQ=255)

                                                cGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTGAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtga                                         <  2:1451952/93‑1 (MQ=255)

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

                                                    aCGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAc                                     >  1:2341245/1‑93 (MQ=255)

                                                    aCGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAc                                     <  2:4604296/93‑1 (MQ=255)

                                                    aCGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAc                                     <  1:4798118/93‑1 (MQ=255)

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

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

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

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

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

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

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              <  1:3803288/93‑1 (MQ=255)

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

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              <  2:1956340/93‑1 (MQ=255)

                                                           aaaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAt                              <  2:5393537/93‑1 (MQ=255)

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

                                                             aaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtc                            <  1:1477928/93‑1 (MQ=255)

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

                                                                gACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTctct                         >  2:2185391/1‑93 (MQ=255)

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

                                                                  cAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGa                       <  2:4989488/93‑1 (MQ=255)

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

                                                                    aTGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATa                     <  2:2505213/93‑1 (MQ=255)

                                                                      gCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaaa                   <  1:642237/93‑1 (MQ=255)

                                                                      gCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATaaa                   >  1:5317056/1‑93 (MQ=255)

                                                                       ccATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAAc                  >  1:2797664/1‑93 (MQ=255)

                                                                          tAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACggg               <  2:6272491/93‑1 (MQ=255)

                                                                          tAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACggg               <  1:2492658/93‑1 (MQ=255)

                                                                          tAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACggg               <  1:3147560/93‑1 (MQ=255)

                                                                            aTTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgcg             <  1:3898389/93‑1 (MQ=255)

                                                                             tttGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgcgc            <  2:1929706/93‑1 (MQ=255)

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

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

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

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

                                                                                 cggcgCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGc        <  2:1003774/93‑1 (MQ=255)

                                                                                 cggcgCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGc        <  1:4779591/93‑1 (MQ=255)

                                                                                   gcgcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTg      <  2:1074035/93‑1 (MQ=255)

                                                                                   gcgcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTg      >  2:4139164/1‑93 (MQ=255)

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

                                                                                    cgcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGt     <  1:3323261/93‑1 (MQ=255)

                                                                                    cgcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGt     <  1:2714972/93‑1 (MQ=255)

                                                                                     gcAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGtt    <  2:5916100/93‑1 (MQ=255)

                                                                                      cAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGttt   >  2:4056181/1‑93 (MQ=255)

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

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