BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	654,113	A→G	L360P (CTG→CCG)	dcuC ←	DcuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	654,113	0	A	G	83.3%	611.7 / 124.2	282	L360P (CTG→CCG)	dcuC	DcuC
Reads supporting (aligned to +/- strand): ref base A (27/20); new base G (116/119); total (143/139)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 286 total aligned reads displayed.

GGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGCCGATG  >  NC_000913/654023‑654205

gggTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGc                                                                                            >  1:4962446/1‑93 (MQ=255)

 ggTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGGCTGCCGCCGGCa                                                                                           >  1:3653511/1‑93 (MQ=255)

 ggTTAATCCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCa                                                                                           <  2:863113/93‑1 (MQ=255)

   ttAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCAtt                                                                                         >  2:1921867/1‑93 (MQ=255)

     aaTGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGt                                                                                       <  2:1192791/93‑1 (MQ=255)

          cGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGaa                                                                                  <  1:2072119/93‑1 (MQ=255)

          cGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGaa                                                                                  >  1:2234695/1‑93 (MQ=255)

           ggAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAAt                                                                                 <  1:5580872/93‑1 (MQ=255)

            gAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATc                                                                                >  2:6729976/1‑93 (MQ=255)

            gAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATc                                                                                >  1:5029470/1‑93 (MQ=255)

            gAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATc                                                                                >  1:1733860/1‑93 (MQ=255)

             aaGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCa                                                                               <  2:758631/93‑1 (MQ=255)

              agagTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCAc                                                                              >  2:2593220/1‑93 (MQ=255)

               gagTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCAcc                                                                             <  2:2941733/93‑1 (MQ=255)

                 gTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCaa                                                                           <  2:4787355/93‑1 (MQ=255)

                 gTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCaa                                                                           <  1:142437/93‑1 (MQ=255)

                 gTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGACAGCATTGTCAGAATCACCaa                                                                           >  2:4657476/1‑93 (MQ=255)

                  tGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAAt                                                                          <  1:7072396/93‑1 (MQ=255)

                  tGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAAt                                                                          >  1:6849659/1‑93 (MQ=255)

                   gcgcCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATa                                                                         <  1:5300339/93‑1 (MQ=255)

                   gcgcCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATa                                                                         <  1:4122419/93‑1 (MQ=255)

                      ccAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCa                                                                      <  2:2881856/93‑1 (MQ=255)

                      ccAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCa                                                                      >  2:1285831/1‑93 (MQ=255)

                      ccAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCa                                                                      <  1:1586302/93‑1 (MQ=255)

                      ccAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCa                                                                      <  2:3463239/93‑1 (MQ=255)

                        aGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGc                                                                    >  2:4237485/1‑93 (MQ=255)

                         gTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCa                                                                   <  2:4670521/93‑1 (MQ=255)

                         gTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCa                                                                   >  1:2864666/1‑93 (MQ=255)

                         gTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCa                                                                   >  2:6986147/1‑93 (MQ=255)

                          tttCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCAt                                                                  <  1:1898015/93‑1 (MQ=255)

                           ttCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGGCGTGACTGCCGCCGGCATCGTCAGAATCACCAATACCAGCATc                                                                 <  1:5729631/93‑1 (MQ=255)

                             cGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAg                                                               <  2:1568194/93‑1 (MQ=255)

                               ggATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGga                                                             <  2:5334145/93‑1 (MQ=255)

                                gATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGgat                                                            >  2:3850/1‑93 (MQ=255)

                                  tcgtctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGgatga                                                          <  1:2548415/90‑1 (MQ=255)

                                  tcatcTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGgatga                                                          <  2:3110457/93‑1 (MQ=255)

                                  tcatcTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGgatga                                                          >  2:3449776/1‑93 (MQ=255)

                                    atcTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATa                                                        >  2:443552/1‑93 (MQ=255)

                                    atcTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATa                                                        <  1:4196337/93‑1 (MQ=255)

                                    atcTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATa                                                        >  1:4674925/1‑93 (MQ=255)

                                     tctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATAc                                                       >  1:937398/1‑93 (MQ=255)

                                     tctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATAc                                                       <  2:3513535/93‑1 (MQ=255)

                                     tctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATAc                                                       >  2:4005092/1‑93 (MQ=255)

                                     tctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATAc                                                       >  2:6573489/1‑93 (MQ=255)

                                     tctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATAc                                                       <  1:122712/93‑1 (MQ=255)

                                      ctcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACt                                                      <  2:549551/93‑1 (MQ=255)

                                       tcAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTc                                                     <  2:1331082/93‑1 (MQ=255)

                                          acaaACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGcc                                                  >  1:321788/1‑93 (MQ=255)

                                          acaaACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGcc                                                  >  2:64048/1‑93 (MQ=255)

                                           caaACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCg                                                 >  2:3630878/1‑93 (MQ=255)

                                           caaACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCg                                                 >  1:6306149/1‑93 (MQ=255)

                                            aaaCGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGa                                                >  2:7152658/1‑93 (MQ=255)

                                            aaaCGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGa                                                >  1:808038/1‑93 (MQ=255)

                                              aCGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAAc                                              <  1:4047041/93‑1 (MQ=255)

                                              aCGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAAc                                              <  2:5850713/93‑1 (MQ=255)

                                                  aTAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCaaa                                          <  1:5144741/93‑1 (MQ=255)

                                                  aTAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCaaa                                          >  1:6522099/1‑93 (MQ=255)

                                                  aTAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCaaa                                          >  2:1106630/1‑93 (MQ=255)

                                                     aaaCGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGa                                       >  2:5045288/1‑93 (MQ=255)

                                                     aaaCGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGa                                       <  1:3020357/93‑1 (MQ=255)

                                                      aaCGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAg                                      <  2:6150475/93‑1 (MQ=255)

                                                      aaCGGCGCATTGCCTGAACCGGTCGCGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAg                                      >  2:6063206/1‑93 (MQ=255)

                                                       aCGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAgg                                     <  2:3389429/93‑1 (MQ=255)

                                                        cGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGt                                    <  2:5090984/93‑1 (MQ=255)

                                                         ggCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTa                                   <  1:6244085/93‑1 (MQ=255)

                                                          gcgcATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAg                                  >  2:2479292/1‑93 (MQ=255)

                                                          gcgcATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAg                                  >  1:2683060/1‑93 (MQ=255)

                                                             cATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGa                               <  2:7018282/93‑1 (MQ=255)

                                                             cATTGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGa                               >  2:3209212/1‑93 (MQ=255)

                                                               ttGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATa                             >  2:7246607/1‑93 (MQ=255)

                                                               ttGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATa                             >  1:2132721/1‑93 (MQ=255)

                                                               ttGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATa                             >  1:4289716/1‑93 (MQ=255)

                                                               ttGCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATa                             >  1:5620897/1‑93 (MQ=255)

                                                               ttGCCTGAACCGGTCGTGACTGCCGCCGGCATTATCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATa                             >  2:7438730/1‑93 (MQ=255)

                                                                 gCCTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGa                           <  1:1355458/93‑1 (MQ=255)

                                                                   cTGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGaaa                         >  1:3458293/1‑93 (MQ=255)

                                                                    tGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAAt                        >  2:1703526/1‑93 (MQ=255)

                                                                    tGAACCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAAt                        >  1:3360789/1‑93 (MQ=255)

                                                                    tGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAAt                        >  2:599515/1‑93 (MQ=255)

                                                                     gAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATc                       >  1:2478430/1‑93 (MQ=255)

                                                                      aaCCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCa                      <  2:4090878/93‑1 (MQ=255)

                                                                      aaCCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCa                      <  1:5934306/93‑1 (MQ=255)

                                                                      aaCCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCa                      >  2:5177836/1‑93 (MQ=255)

                                                                      aaCCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCa                      <  2:2727777/93‑1 (MQ=255)

                                                                       aCCGGTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAg                     <  1:4478883/93‑1 (MQ=255)

                                                                           gTCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACtt                 <  2:4145857/93‑1 (MQ=255)

                                                                            tCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACttt                <  1:1002889/93‑1 (MQ=255)

                                                                            tCGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACttt                <  1:4943956/93‑1 (MQ=255)

                                                                             cGTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTg               <  1:4871680/93‑1 (MQ=255)

                                                                              gTGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGa              <  2:5938934/93‑1 (MQ=255)

                                                                               tGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGaa             <  1:7166012/93‑1 (MQ=255)

                                                                               tGACTGCCGCCGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGaa             >  1:6819579/1‑93 (MQ=255)

                                                                                    gccgccGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAg        >  1:4545455/1‑93 (MQ=255)

                                                                                    gccgccGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAg        <  2:7322017/93‑1 (MQ=255)

                                                                                    gccgccAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAg        <  2:3789914/93‑1 (MQ=255)

                                                                                     ccgccgGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGc       <  1:3814529/93‑1 (MQ=255)

                                                                                      cgccGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGcc      <  2:48418/93‑1 (MQ=255)

                                                                                      cgccAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGcc      <  1:617630/93‑1 (MQ=255)

                                                                                       gccGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGCCg     >  1:4730615/1‑93 (MQ=255)

                                                                                       gccGGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGCCg     >  2:2250789/1‑93 (MQ=255)

GGGTTAATGCCGGAAGAGTGCGCCAGTTTCGGGATCATCTCAACAAACGCATAAAACGGCGCATTGCCTGAACCGGTCGTGACTGCCGCCAGCATTGTCAGAATCACCAATACCAGCATCAGGATGATACTCGCCGAACCAAACGAGGTAGCGATAGAAATCAGACTTTGAATAAAGCCGATG  >  NC_000913/654023‑654205

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