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.3%	422.0 / ‑4.8	137	Q31R (CAG→CGG)	b1838	protein phosphatase 1
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (56/80); total (56/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

Only 100 of 138 total aligned reads displayed.

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

cTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGc                                                                                         <  1:1353379/93‑1 (MQ=255)

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

 tttGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCt                                                                                        >  2:2632763/1‑93 (MQ=255)

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

     cGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAg                                                                                    <  1:2095303/93‑1 (MQ=255)

         ccACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACaa                                                                                <  2:1620895/93‑1 (MQ=255)

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

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

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

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

           aCGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAAcc                                                                              <  1:1274628/93‑1 (MQ=255)

            cgatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCg                                                                             >  2:472039/1‑93 (MQ=255)

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

            cgatcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCg                                                                             >  2:4627852/1‑93 (MQ=255)

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

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

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

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

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

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

              agcgatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCgtg                                                                           >  2:4616876/3‑93 (MQ=255)

                cgatAACGTCTCTCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGta                                                                         <  1:2387806/93‑1 (MQ=255)

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

                 gatAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGGTCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCATGtat                                                                        >  2:3849679/1‑93 (MQ=255)

                  atAACGTTTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGtata                                                                       >  2:4454477/1‑93 (MQ=255)

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

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

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

                   tAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTatat                                                                      <  2:2520071/93‑1 (MQ=255)

                   tAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTatat                                                                      <  2:3827551/93‑1 (MQ=255)

                    aaCGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATc                                                                     <  1:3990043/93‑1 (MQ=255)

                       gTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGcc                                                                  <  1:2781179/93‑1 (MQ=255)

                       gTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGcc                                                                  <  2:1882990/93‑1 (MQ=255)

                         ctcCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAg                                                                <  2:2770692/93‑1 (MQ=255)

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

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

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

                              aCTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGATAGc                                                           >  2:1280300/1‑93 (MQ=255)

                                tGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCTGAAAGCCa                                                         <  2:4555601/93‑1 (MQ=255)

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

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

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

                                   aaTAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGat                                                      <  1:4406937/92‑1 (MQ=255)

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

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

                                     tAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTAGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGatat                                                    >  1:3294825/1‑93 (MQ=255)

                                       aGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                                  <  1:3965708/93‑1 (MQ=255)

                                       aGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                                  <  2:3310366/93‑1 (MQ=255)

                                       aGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGt                                                  <  1:3999241/93‑1 (MQ=255)

                                       aGTAAATCACGCGACGGATCAAAACGACAATGCCATAATTGGCTGCGCAGCCGCTCAAGACAAGCGTGTATATCGCGAGAAAGCCAGTTATGt                                                  <  2:3980140/93‑1 (MQ=255)

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

                                          aaaTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGc                                               >  1:2613988/1‑93 (MQ=255)

                                           aaTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGcc                                              <  2:4161142/93‑1 (MQ=255)

                                            aTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCa                                             >  2:2408744/1‑93 (MQ=255)

                                            aTCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCa                                             >  2:1502578/1‑93 (MQ=255)

                                             tctcGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCAt                                            <  2:1976811/93‑1 (MQ=255)

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

                                                 gCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgat                                        <  2:3208182/93‑1 (MQ=255)

                                                   cACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACACCCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATtgatga                                      <  2:1825290/93‑1 (MQ=255)

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

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

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

                                                     cGGATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAcc                                    >  1:1881932/1‑93 (MQ=255)

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

                                                      ggATCAAAACGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGAccc                                   <  1:3681357/93‑1 (MQ=255)

                                                        aTCAAAACGACAATGCCATTATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGc                                 >  2:4582602/1‑93 (MQ=255)

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

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

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

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

                                                            aaaCGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAAtt                             <  2:1096087/93‑1 (MQ=255)

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

                                                               cGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc                          <  1:721805/93‑1 (MQ=255)

                                                               cGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc                          <  2:1206475/93‑1 (MQ=255)

                                                               cGACAATGCCATAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATtctc                          >  2:981990/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                          tAATTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACggg               >  1:1629595/1‑93 (MQ=255)

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

                                                                           aaTTTGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgc              <  1:764230/93‑1 (MQ=255)

                                                                            atgtGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGgcg             <  1:1686711/90‑1 (MQ=255)

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

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

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

                                                                               tGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAg          <  2:985916/93‑1 (MQ=255)

                                                                               tGCGGCGCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAg          <  1:42514/93‑1 (MQ=255)

                                                                                gcggcgCAACCGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAgg         >  1:2195246/1‑93 (MQ=255)

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

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

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

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

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

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

CTTTGCGGCCCACGATCGATAACGTCTCCCACTGAGATAAGTAAATCTCGCCACGGATCAAAACGACAATGCCATAATTTGCGGCGCAACTGCTCAAGACAACCGTGTATATCGCCAGAAAGCCAGATATGTCGCCATTGATGACCCGCAATTCTCTGATAAACGGGCGCAGGCTGTTTC  >  NC_000913/1920812‑1920991

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