BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	800,958	A→G	Q326R (CAG→CGG)	ybhH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	800,958	0	A	G	99.1%	697.3 / ‑5.2	225	Q326R (CAG→CGG)	ybhH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (114/108); total (117/108)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

CCCGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCAGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTGT  >  NC_000913/800868‑801047

cccGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCggg                                                                                         >  2:6236399/1‑93 (MQ=255)

 cctaCAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTACTGAAGGTCGGGa                                                                                        <  2:3068503/90‑1 (MQ=255)

 ccGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGa                                                                                        >  2:6645711/1‑93 (MQ=255)

 ccGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGa                                                                                        >  1:3089609/1‑93 (MQ=255)

 ccGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGa                                                                                        <  1:7065455/93‑1 (MQ=255)

 ccGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGa                                                                                        <  2:2240275/93‑1 (MQ=255)

  cGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGAt                                                                                       >  2:512129/1‑93 (MQ=255)

   gACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATg                                                                                      <  1:4982086/93‑1 (MQ=255)

   gACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATg                                                                                      <  2:3540738/93‑1 (MQ=255)

     cAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGCGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGcc                                                                                    >  1:2390068/1‑93 (MQ=255)

      aaaTCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGccc                                                                                   >  2:6059510/1‑92 (MQ=255)

      aaaTCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCa                                                                                   <  1:4950109/93‑1 (MQ=255)

           gTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCacgacg                                                                              >  1:4392308/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             >  1:104812/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             >  1:4816361/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             >  1:4342890/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             >  1:3895356/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             >  1:133629/1‑93 (MQ=255)

            tCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGt                                                                             <  1:2667916/93‑1 (MQ=255)

             cccTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGCGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGtt                                                                            >  1:4668486/1‑93 (MQ=255)

              ccTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTa                                                                           >  2:4991388/1‑93 (MQ=255)

               cTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAc                                                                          <  1:5605389/93‑1 (MQ=255)

               cTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGCGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAc                                                                          <  2:6698774/93‑1 (MQ=255)

                ttCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAcg                                                                         >  1:1591528/1‑93 (MQ=255)

                ttCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTAcg                                                                         >  2:4533085/1‑93 (MQ=255)

                   tGTAGGATACGGCAATATCAATATTGAAAACCCCAGCGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCCCGACGTTACgcgc                                                                      >  2:6752291/1‑93 (MQ=255)

                    gTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCa                                                                     >  1:4016019/1‑93 (MQ=255)

                    gTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCa                                                                     <  1:2052231/93‑1 (MQ=255)

                      aGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATc                                                                   <  1:791182/93‑1 (MQ=255)

                       ggATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCt                                                                  >  2:5278297/1‑93 (MQ=255)

                       ggATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCt                                                                  <  2:5517845/93‑1 (MQ=255)

                        gATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTg                                                                 <  2:1983558/93‑1 (MQ=255)

                        gATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTg                                                                 >  2:2569590/1‑93 (MQ=255)

                              gCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTc                                                           <  1:110882/93‑1 (MQ=255)

                               caatatcaatatTGAACACCCCAGCGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCg                                                          >  2:1942396/1‑93 (MQ=255)

                                aatatcaatatTGAACACCCCAGTGGTGCTCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCgg                                                         <  1:2448472/93‑1 (MQ=255)

                                aatatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCgg                                                         >  2:2832271/1‑93 (MQ=255)

                                 atatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGga                                                        >  2:6848516/1‑93 (MQ=255)

                                 atatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGga                                                        <  1:1919100/93‑1 (MQ=255)

                                  tatcaatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgac                                                       <  2:7350604/93‑1 (MQ=255)

                                     caatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacgac                                                    <  2:2928201/93‑1 (MQ=255)

                                     caatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGgacgac                                                    <  1:1754005/93‑1 (MQ=255)

                                      aatatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGAcc                                                   <  1:3177510/93‑1 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGGTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  <  2:4283912/93‑1 (MQ=255)

                                       atatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCa                                                  <  1:4415732/93‑1 (MQ=255)

                                        tatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAg                                                 >  2:1445771/1‑93 (MQ=255)

                                        tatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAg                                                 >  1:141924/1‑93 (MQ=255)

                                        tatTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAg                                                 <  2:2154625/93‑1 (MQ=255)

                                          ttGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaa                                               <  2:4193295/93‑1 (MQ=255)

                                          ttGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaa                                               <  1:7799459/93‑1 (MQ=255)

                                           tGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaa                                              <  1:2534546/93‑1 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaaaa                                            >  1:3230962/1‑93 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaaaa                                            <  2:6954389/93‑1 (MQ=255)

                                             aaCACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGaaaaa                                            <  1:2191398/93‑1 (MQ=255)

                                                acCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAatat                                         <  1:5849984/93‑1 (MQ=255)

                                                 ccccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtt                                        >  1:3027580/1‑93 (MQ=255)

                                                 ccccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtt                                        <  2:4776142/93‑1 (MQ=255)

                                                 ccccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtt                                        <  2:6988643/93‑1 (MQ=255)

                                                   ccAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATAtttt                                      <  2:7292136/93‑1 (MQ=255)

                                                    cAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTc                                     >  2:3680082/1‑93 (MQ=255)

                                                     aGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTcc                                    <  1:3417458/93‑1 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   <  1:5460742/93‑1 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   >  2:7671765/1‑93 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   >  2:3488691/1‑93 (MQ=255)

                                                      gtggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCg                                   >  2:3612254/1‑93 (MQ=255)

                                                        ggtgCGCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGt                                 <  2:4057916/93‑1 (MQ=255)

                                                          tgCTCTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGa                               <  2:4631263/93‑1 (MQ=255)

                                                          tgCGCTCGCCGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGa                               <  2:7851867/93‑1 (MQ=255)

                                                              cTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGtt                           <  2:115411/93‑1 (MQ=255)

                                                              cTCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGtt                           <  1:95465/93‑1 (MQ=255)

                                                               tCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGttt                          >  2:1924178/1‑93 (MQ=255)

                                                               tCGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGttt                          >  1:3639007/1‑93 (MQ=255)

                                                                cGACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTa                         <  1:6295026/93‑1 (MQ=255)

                                                                 gACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTAt                        >  1:2513373/1‑93 (MQ=255)

                                                                 gACGTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTAt                        >  1:6756445/1‑93 (MQ=255)

                                                                    gTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCtt                     <  1:1004960/93‑1 (MQ=255)

                                                                    gTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCtt                     >  2:4360924/1‑93 (MQ=255)

                                                                    gTTCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCtt                     >  1:4880699/1‑93 (MQ=255)

                                                                      tCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTcc                   >  1:5176567/1‑93 (MQ=255)

                                                                      tCATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTcc                   >  2:3864552/1‑93 (MQ=255)

                                                                       cATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTccc                  <  1:457778/93‑1 (MQ=255)

                                                                       cATTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTccc                  <  1:2897743/93‑1 (MQ=255)

                                                                        aTTTAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCt                 >  2:4880180/1‑93 (MQ=255)

                                                                         tttAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTg                <  1:2969445/93‑1 (MQ=255)

                                                                          ttAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGa               >  1:6504424/1‑93 (MQ=255)

                                                                          ttAAGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGa               <  1:4013596/93‑1 (MQ=255)

                                                                             aGTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaa            >  2:5622679/1‑93 (MQ=255)

                                                                              gTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaa           >  2:6281400/1‑93 (MQ=255)

                                                                              gTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaa           >  1:694336/1‑93 (MQ=255)

                                                                              gTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaa           <  2:4949888/93‑1 (MQ=255)

                                                                              gTAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaa           <  2:7833943/93‑1 (MQ=255)

                                                                               tAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaaa          >  1:6618249/1‑93 (MQ=255)

                                                                               tAATGAAGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGaaaaaa          >  1:743458/1‑93 (MQ=255)

                                                                                aaTGAAGGTCGGGATGCCATGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAAt         <  2:589436/93‑1 (MQ=255)

                                                                                   gcaGGTCAGGATGCCACGACGTTACGCGCATCTGTTATTCCGACGCCCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGACAACATTct      >  2:3913280/3‑92 (MQ=255)

                                                                                    aaGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGt     <  2:7482327/93‑1 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGATAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    <  2:1854067/93‑1 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    <  1:658387/93‑1 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    >  1:2803192/1‑93 (MQ=255)

                                                                                     aGGTCGGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGtt    <  2:1314568/93‑1 (MQ=255)

CCCGACAAATCGTCCCTTCTGTAGGATACGGCAATATCAATATTGAACACCCCAGTGGTGCGCTCGACGTTCATTTAAGTAATGAAGGTCAGGATGCCACGACGTTACGCGCATCTGTTATTCGGACGACCAGAAAAATATTTTCCGGTGAAGTTTATCTTCCCTGAAAAAATTCGTTGT  >  NC_000913/800868‑801047

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: `‑`