BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,093,303	+G	coding (187/705 nt)	yggS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,093,299	1	.	G	99.0%	1080.2 / 0.8	315	coding (183/705 nt)	yggS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base G (132/179); total (134/181)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 367 total aligned reads displayed.

CTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑GGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAGC  >  NC_000913/3093207‑3093386

ctgctTGCAGTCAGTAAAACAAACCCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGaa                                                                                          <  1:6218429/93‑1 (MQ=255)

 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑g                                                                                        >  2:3881005/1‑93 (MQ=255)

 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑g                                                                                        <  1:983715/93‑1 (MQ=255)

 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑g                                                                                        <  2:2148549/93‑1 (MQ=255)

 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑g                                                                                        >  1:2276355/1‑93 (MQ=255)

  gctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gg                                                                                       >  2:4999738/1‑93 (MQ=255)

  gctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gg                                                                                       <  2:2406993/93‑1 (MQ=255)

    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggg                                                                                     >  2:5262910/1‑93 (MQ=255)

    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggg                                                                                     >  2:7484738/1‑93 (MQ=255)

    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAACCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggg                                                                                     >  1:3810763/1‑93 (MQ=255)

      gCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggggt                                                                                   <  1:7389444/93‑3 (MQ=255)

       cagtcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggggta                                                                                  >  1:8365636/1‑90 (MQ=255)

       cagtcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggggta                                                                                  <  1:4246438/93‑4 (MQ=255)

        agtcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑gggggtag                                                                                 >  2:1779824/1‑89 (MQ=255)

          tcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGAt                                                                                >  2:3945116/1‑93 (MQ=255)

              tAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATaaaa                                                                            <  1:2590847/93‑1 (MQ=255)

               aaaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAt                                                                           >  2:4335498/1‑93 (MQ=255)

                aaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAtt                                                                          <  2:7250808/93‑1 (MQ=255)

                aaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAtt                                                                          >  2:3727065/1‑93 (MQ=255)

                aaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAtt                                                                          <  1:2093675/93‑1 (MQ=255)

                 aacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTc                                                                         >  2:8168330/1‑93 (MQ=255)

                  acaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCg                                                                        >  1:4401361/1‑93 (MQ=255)

                  acaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCg                                                                        >  1:2760774/1‑93 (MQ=255)

                  acaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCg                                                                        <  1:5159890/93‑1 (MQ=255)

                     aaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCa                                                                     <  1:4750199/93‑1 (MQ=255)

                     aaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCa                                                                     >  1:1329425/1‑93 (MQ=255)

                      aacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCAc                                                                    <  2:1634056/93‑1 (MQ=255)

                       acCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACt                                                                   >  1:5839988/1‑93 (MQ=255)

                        ccTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACtt                                                                  <  2:7999002/93‑1 (MQ=255)

                        ccTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACtt                                                                  <  1:247077/93‑1 (MQ=255)

                          tGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTc                                                                <  1:6015872/93‑1 (MQ=255)

                              aGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGa                                                            <  2:3190937/93‑1 (MQ=255)

                               gcgcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGaa                                                           <  2:7201359/93‑1 (MQ=255)

                                ggcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTGGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAAc                                                          <  1:4977631/92‑1 (MQ=255)

                                 gcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACt                                                         <  2:3435694/93‑1 (MQ=255)

                                  ccATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                        <  2:6472770/93‑1 (MQ=255)

                                  ccATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                        >  1:8005242/1‑93 (MQ=255)

                                  ccATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                        <  1:8810375/93‑1 (MQ=255)

                                  ccATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                        >  1:8731168/1‑93 (MQ=255)

                                   taTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAATACTACGTTCAGGAAGGGGGCAGATAAAATACGCCACTTTCAGGAACTgg                                                       >  2:7806627/2‑93 (MQ=255)

                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                       >  1:7157507/1‑93 (MQ=255)

                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                       <  1:7103031/93‑1 (MQ=255)

                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                       <  1:2699113/93‑1 (MQ=255)

                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                       >  2:134623/1‑93 (MQ=255)

                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                      >  1:8336159/1‑93 (MQ=255)

                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                      <  2:3364536/93‑1 (MQ=255)

                                      cGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCg                                                    <  1:5639774/93‑1 (MQ=255)

                                         aGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTaa                                                 <  1:6860483/93‑1 (MQ=255)

                                          gAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAAc                                                >  2:741292/1‑93 (MQ=255)

                                          gAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAAc                                                <  2:7727879/93‑1 (MQ=255)

                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                               >  1:1665113/1‑93 (MQ=255)

                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                               >  1:784393/1‑93 (MQ=255)

                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                               <  1:1541462/93‑1 (MQ=255)

                                             gCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAgg                                             <  2:6051517/93‑1 (MQ=255)

                                             gCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAgg                                             <  2:2419393/93‑1 (MQ=255)

                                                aTTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAtt                                          >  1:815979/1‑93 (MQ=255)

                                                aTTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAtt                                          >  1:4588335/1‑93 (MQ=255)

                                                aTTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAtt                                          >  2:622082/1‑93 (MQ=255)

                                                   gATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGa                                       >  1:3612000/1‑93 (MQ=255)

                                                   gATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGa                                       <  2:8302474/93‑1 (MQ=255)

                                                   gATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGa                                       <  1:324050/93‑1 (MQ=255)

                                                     tggcGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAAt                                     <  1:3036569/90‑1 (MQ=255)

                                                     tGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAAt                                     <  2:2640479/93‑1 (MQ=255)

                                                         gggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCa                                 >  1:5585551/1‑93 (MQ=255)

                                                         gggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCa                                 <  1:3774912/93‑1 (MQ=255)

                                                         gggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCa                                 >  2:5343179/1‑93 (MQ=255)

                                                            cAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAGCAGGATTAGAATGGCAttt                              <  1:2305992/93‑1 (MQ=255)

                                                              gCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTa                            <  2:8686608/93‑1 (MQ=255)

                                                              gCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTa                            >  2:7032280/1‑93 (MQ=255)

                                                                gTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTAtt                          >  2:6009968/1‑93 (MQ=255)

                                                                gTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTAtt                          <  2:5130583/93‑1 (MQ=255)

                                                                 tCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTg                         >  1:6964632/1‑93 (MQ=255)

                                                                 tCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTg                         <  2:3595276/93‑1 (MQ=255)

                                                                  cAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTgg                        <  1:74943/93‑1 (MQ=255)

                                                                  cAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTgg                        <  2:2791558/93‑1 (MQ=255)

                                                                   aaTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGc                       >  1:3152747/1‑93 (MQ=255)

                                                                    aTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGcc                      >  2:1465245/1‑93 (MQ=255)

                                                                     tttGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGccc                     <  2:8030306/93‑1 (MQ=255)

                                                                      ttGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCg                    <  2:3255139/93‑1 (MQ=255)

                                                                       tGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGt                   >  1:8498555/1‑93 (MQ=255)

                                                                       tGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGt                   >  2:2013670/1‑93 (MQ=255)

                                                                         gTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTg                 <  1:4828228/93‑1 (MQ=255)

                                                                         gTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTg                 <  1:1918754/93‑1 (MQ=255)

                                                                            aaaaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAg              >  2:4948793/1‑93 (MQ=255)

                                                                              aaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTc            <  2:5488336/93‑1 (MQ=255)

                                                                              aaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTc            >  2:6337994/1‑93 (MQ=255)

                                                                               aCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCATTCt           >  1:1218592/1‑93 (MQ=255)

                                                                                cTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCta          <  2:587404/93‑1 (MQ=255)

                                                                                cTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCta          <  1:5431405/93‑1 (MQ=255)

                                                                                atACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCta          <  2:6768913/92‑1 (MQ=255)

                                                                                 tACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaa         <  1:5878857/93‑1 (MQ=255)

                                                                                 tACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaa         <  2:8336759/93‑1 (MQ=255)

                                                                                  aCGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaat        <  2:4249567/93‑1 (MQ=255)

                                                                                  aCGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaat        <  2:4524242/93‑1 (MQ=255)

                                                                                      tCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaaa    >  2:4792705/1‑93 (MQ=255)

                                                                                      tCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaaa    <  2:2712301/93‑1 (MQ=255)

                                                                                      tCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaaa    <  1:6681487/93‑1 (MQ=255)

                                                                                       cAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATTAAAg   <  2:5599853/93‑1 (MQ=255)

                                                                                       cAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAg   <  2:3050779/93‑1 (MQ=255)

                                                                                        aGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAgc  >  2:5770226/1‑93 (MQ=255)

CTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAA‑GGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAGC  >  NC_000913/3093207‑3093386

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