BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	260,240	G→A	G210E (GGA→GAA)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,240	0	G	A	99.2%	956.4 / ‑3.6	251	G210E (GGA→GAA)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (127/121); total (130/121)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 255 total aligned reads displayed.

GCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGG  >  NC_000913/260150‑260329

  agcaATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACt                                                                                       >  1:4795961/1‑93 (MQ=255)

     aaTCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTggc                                                                                    <  1:5559225/93‑1 (MQ=255)

      aTCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTggcg                                                                                   <  2:1551785/93‑1 (MQ=255)

      aTCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTggcg                                                                                   >  1:6221506/1‑93 (MQ=255)

      aTCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTggcg                                                                                   >  1:6012157/1‑93 (MQ=255)

        ccGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTggcggc                                                                                 <  2:3026704/93‑1 (MQ=255)

          gcaggcagAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCAt                                                                               >  1:110245/1‑93 (MQ=255)

             ggcagAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAg                                                                            >  1:4325531/1‑93 (MQ=255)

                agAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTAc                                                                         <  2:1158288/93‑1 (MQ=255)

                agAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTAc                                                                         >  1:1628810/1‑93 (MQ=255)

                agAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTAc                                                                         <  1:1666415/93‑1 (MQ=255)

                agAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTAc                                                                         >  2:5338874/1‑93 (MQ=255)

                 gAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTAcc                                                                        >  1:5496310/1‑93 (MQ=255)

                 gAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTAcc                                                                        >  1:2139418/1‑93 (MQ=255)

                   aCTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCaa                                                                      >  2:777643/1‑93 (MQ=255)

                   aCTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCaa                                                                      >  2:209299/1‑93 (MQ=255)

                      gATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAAtt                                                                   >  1:5862301/1‑93 (MQ=255)

                      gATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAAtt                                                                   >  2:2278449/1‑93 (MQ=255)

                       aTTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTg                                                                  >  2:1810278/1‑93 (MQ=255)

                        ttAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGc                                                                 <  2:5929569/93‑1 (MQ=255)

                          aaaGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAg                                                               >  2:158929/1‑93 (MQ=255)

                           aaGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAgg                                                              >  1:3552253/1‑93 (MQ=255)

                            aGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGCAACTGGCGGCATGAGTACCAAATTGCAGGc                                                             >  1:3405520/1‑93 (MQ=255)

                            aGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGc                                                             <  1:5765190/93‑1 (MQ=255)

                             gATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGcc                                                            <  1:1302257/93‑1 (MQ=255)

                              aTGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCg                                                           >  2:3371412/1‑93 (MQ=255)

                              aTGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCg                                                           >  1:926143/1‑93 (MQ=255)

                               tGTTGAGGGGATTGATGCCGCGCGGGGGGCGTGTGGCGGGGACGGGTTTTGAGCCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGc                                                          <  2:4862553/93‑1 (MQ=255)

                                gTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCt                                                         >  2:4108995/1‑93 (MQ=255)

                                gTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCt                                                         <  1:3233575/93‑1 (MQ=255)

                                   tACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGAc                                                      >  2:269112/1‑93 (MQ=255)

                                   tACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGAc                                                      <  1:6420359/93‑1 (MQ=255)

                                   tACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAACTGCAGGCCGCTGAc                                                      >  1:3618073/1‑93 (MQ=255)

                                    aCGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACg                                                     >  2:6446012/1‑93 (MQ=255)

                                     cGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGt                                                    <  2:5393774/93‑1 (MQ=255)

                                     cGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGt                                                    <  1:5681776/93‑1 (MQ=255)

                                     cGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCCCTGACGt                                                    <  1:706476/93‑1 (MQ=255)

                                      ggCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTg                                                   >  1:5070588/1‑93 (MQ=255)

                                      ggCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTg                                                   >  1:3899076/1‑93 (MQ=255)

                                       gCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTgg                                                  >  2:2458990/1‑93 (MQ=255)

                                        cATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGc                                                 <  2:3746022/93‑1 (MQ=255)

                                        cATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGc                                                 <  2:1709754/93‑1 (MQ=255)

                                         aTTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCt                                                <  2:3213639/93‑1 (MQ=255)

                                         aTTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCt                                                <  1:546199/93‑1 (MQ=255)

                                          ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt                                               <  2:1435223/93‑1 (MQ=255)

                                          ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt                                               <  1:4002065/93‑1 (MQ=255)

                                          ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt                                               <  2:4566151/93‑1 (MQ=255)

                                          ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt                                               <  2:165293/93‑1 (MQ=255)

                                          ttGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt                                               <  2:6008859/93‑1 (MQ=255)

                                            gatgaCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGc                                             <  1:3365251/93‑1 (MQ=255)

                                             atgaCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGcc                                            <  2:2141382/93‑1 (MQ=255)

                                             atgaCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGcc                                            >  2:2658530/1‑93 (MQ=255)

                                                aCGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTg                                         >  2:5042562/1‑93 (MQ=255)

                                                 cGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGc                                        <  2:6422119/93‑1 (MQ=255)

                                                  gCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCg                                       <  1:1469381/93‑1 (MQ=255)

                                                    aCTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCggg                                     >  2:2547564/1‑93 (MQ=255)

                                                    aCTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCggg                                     <  2:1480912/93‑1 (MQ=255)

                                                     cTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGTCGTGGCTTGCTGTGCGGGt                                    <  1:6319817/93‑1 (MQ=255)

                                                     cTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGt                                    >  1:1573005/1‑93 (MQ=255)

                                                      tGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTa                                   <  2:2136660/93‑1 (MQ=255)

                                                        cgcgcgATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATc                                 <  1:2639982/93‑1 (MQ=255)

                                                         gcgcgATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCg                                >  2:1988564/1‑93 (MQ=255)

                                                         gcgcgATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCg                                <  2:5798789/93‑1 (MQ=255)

                                                            cgATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGaca                             <  2:5261587/93‑1 (MQ=255)

                                                             gATTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGacac                            >  1:2845979/1‑93 (MQ=255)

                                                              aTTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACAcc                           >  1:1790625/1‑93 (MQ=255)

                                                              aTTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACAcc                           >  1:4492232/1‑93 (MQ=255)

                                                              aTTGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACAcc                           >  1:3848645/1‑93 (MQ=255)

                                                               ttGCCGGTGACAGCGTTTCAGGCCTCGAAGCTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCa                          <  2:6141775/93‑1 (MQ=255)

                                                               ttGCCGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCa                          >  1:2490975/1‑93 (MQ=255)

                                                                  ccGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCatta                       <  1:4268405/93‑1 (MQ=255)

                                                                   cGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCattat                      <  1:1490670/93‑1 (MQ=255)

                                                                   cGGTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCattat                      <  1:491687/93‑1 (MQ=255)

                                                                    ggTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCattatt                     <  1:5627306/93‑1 (MQ=255)

                                                                     gTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCCTTATTg                    >  2:3567240/1‑93 (MQ=255)

                                                                     gTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTg                    >  2:4286155/1‑93 (MQ=255)

                                                                     gTGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTg                    >  1:2893120/1‑93 (MQ=255)

                                                                      tGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGc                   <  2:1293113/93‑1 (MQ=255)

                                                                      tGACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGc                   <  2:3882855/93‑1 (MQ=255)

                                                                       gACAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGcc                  <  1:6041564/93‑1 (MQ=255)

                                                                        aCAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCcg                 >  1:610264/1‑93 (MQ=255)

                                                                        aCAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCcg                 <  1:2255917/93‑1 (MQ=255)

                                                                        aCAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCcg                 <  1:4444181/93‑1 (MQ=255)

                                                                         cAGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCcgc                >  1:896516/1‑93 (MQ=255)

                                                                          aGCGTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCcgcg               >  1:5706227/1‑93 (MQ=255)

                                                                             gTTTCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGgc            >  1:1902298/1‑93 (MQ=255)

                                                                              tttCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGgca           >  1:1959125/1‑93 (MQ=255)

                                                                              tttCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGgca           <  1:2807419/93‑1 (MQ=255)

                                                                              tttAAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGgca           <  1:538080/93‑1 (MQ=255)

                                                                                tCAGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGgcagc         <  2:3148291/93‑1 (MQ=255)

                                                                                  aGGCCTCGAAGCTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCaa       <  1:3203198/93‑1 (MQ=255)

                                                                                  aGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCaa       >  1:117073/1‑93 (MQ=255)

                                                                                  aGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCaa       <  1:4171015/93‑1 (MQ=255)

                                                                                  aGGCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCaa       <  2:3854922/93‑1 (MQ=255)

                                                                                   ggCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAg      >  1:2968109/1‑93 (MQ=255)

                                                                                   ggCCGCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAg      <  1:5891141/93‑1 (MQ=255)

                                                                                    gCCTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGc     <  2:4345111/93‑1 (MQ=255)

                                                                                      cTCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCg   >  1:5570010/1‑93 (MQ=255)

                                                                                       tCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCgg  >  2:3103481/1‑93 (MQ=255)

                                                                                       tCGAAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCgg  <  1:1351282/93‑1 (MQ=255)

GCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGG  >  NC_000913/260150‑260329

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