BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,799	T→C	V462A (GTA→GCA)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,799	0	T	C	98.8%	525.6 / ‑3.2	164	V462A (GTA→GCA)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (82/80); total (83/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 = 7.95e-01

Only 100 of 168 total aligned reads displayed.

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         <  2:2228588/93‑1 (MQ=255)

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         <  2:2228589/93‑1 (MQ=255)

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         <  1:2084481/93‑1 (MQ=255)

tggtggAAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTACGACTTCCGCGCAGAACAAAGCATCGGCCTGGCAt                                                                                         >  2:4591795/1‑93 (MQ=255)

 ggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAtt                                                                                        <  1:860979/93‑1 (MQ=255)

 ggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAtt                                                                                        <  2:700980/93‑1 (MQ=255)

 ggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAtt                                                                                        <  2:3848886/93‑1 (MQ=255)

      cAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCg                                                                                   <  1:1088624/93‑1 (MQ=255)

       aGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCgg                                                                                  <  2:4381303/93‑1 (MQ=255)

       aGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCgg                                                                                  >  1:4614237/1‑93 (MQ=255)

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 <  2:1024337/93‑1 (MQ=255)

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 <  2:1024336/93‑1 (MQ=255)

         tACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTa                                                                                >  2:4729572/1‑93 (MQ=255)

          aCTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTaa                                                                               <  1:4611431/93‑1 (MQ=255)

          aCTGCGCACCGGCCAGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTaa                                                                               <  1:692052/93‑1 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             <  2:3604108/93‑1 (MQ=255)

             gcgcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAggg                                                                            >  2:476076/1‑93 (MQ=255)

             gcgcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAggg                                                                            <  2:1813280/93‑1 (MQ=255)

             gcgcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAggg                                                                            <  2:1280742/93‑1 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          >  2:3770289/1‑93 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          >  2:939560/1‑93 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          <  1:4682720/93‑1 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        >  1:920818/1‑93 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  1:4187717/93‑1 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  1:403929/93‑1 (MQ=255)

                   cGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAatat                                                                      <  1:500322/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCGTCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  1:4107708/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     >  2:1556318/1‑93 (MQ=255)

                    agCCCGACCGAAAGCCTGATTCAGGCGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACCTATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  2:1782381/92‑1 (MQ=255)

                         gaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTc                                                                >  1:4511402/1‑93 (MQ=255)

                         gaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTc                                                                >  2:1181678/1‑93 (MQ=255)

                          accgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTcc                                                               >  1:4083841/1‑93 (MQ=255)

                          accgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGATTcc                                                               >  1:4452771/1‑93 (MQ=255)

                            cgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGGTCCCg                                                             >  2:2004472/1‑93 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            >  1:3684545/1‑93 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            <  1:3679186/93‑1 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            <  2:348210/93‑1 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            <  1:949362/93‑1 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            <  1:1378277/93‑1 (MQ=255)

                               aaGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGCCAATATCGGTTCCCGCTg                                                          >  1:2117153/1‑93 (MQ=255)

                                aGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTgg                                                         >  1:3762346/1‑93 (MQ=255)

                                 gCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGc                                                        >  1:426679/1‑93 (MQ=255)

                                 gCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGc                                                        >  2:4365646/1‑93 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      <  2:2380199/93‑1 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      >  2:2415610/1‑93 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      <  1:982523/93‑1 (MQ=255)

                                    tGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGCAAGGGCAATATCGGTTCCCGTTGGCTgg                                                     >  1:1002532/1‑93 (MQ=255)

                                      aTTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGaa                                                   >  1:333920/1‑93 (MQ=255)

                                       ttCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATAGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAc                                                  <  2:2157269/93‑1 (MQ=255)

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              >  1:2051836/1‑93 (MQ=255)

                                            gggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTc                                             <  2:2047264/93‑1 (MQ=255)

                                              gCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGc                                           <  2:367090/93‑1 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  1:2767024/1‑93 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  2:4010372/1‑93 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  1:2012651/1‑93 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCg                                        <  1:3844969/93‑1 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCg                                        >  2:2963431/1‑93 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCg                                        <  1:3434733/93‑1 (MQ=255)

                                                  cATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGt                                       <  1:2898337/93‑1 (MQ=255)

                                                  cATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGt                                       >  2:1875188/1‑93 (MQ=255)

                                                   aTCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTg                                      >  2:1469506/1‑93 (MQ=255)

                                                   aTCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTg                                      >  2:2321576/1‑93 (MQ=255)

                                                    tCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCACTATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     >  1:2282487/1‑93 (MQ=255)

                                                    tCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     >  1:2597902/1‑93 (MQ=255)

                                                    tCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     >  1:2282486/1‑93 (MQ=255)

                                                       gTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagca                                  >  1:2412814/1‑93 (MQ=255)

                                                        tCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcag                                 >  1:279878/1‑93 (MQ=255)

                                                         ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcaga                                >  1:4251269/1‑93 (MQ=255)

                                                           gTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagc                              >  1:2731811/1‑93 (MQ=255)

                                                            tCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagca                             <  1:219438/93‑1 (MQ=255)

                                                            tCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagca                             >  2:1279209/1‑93 (MQ=255)

                                                              ttCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACg                           >  2:779951/1‑93 (MQ=255)

                                                              ttCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACg                           >  2:2922498/1‑93 (MQ=255)

                                                               tCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGc                          <  1:2518506/93‑1 (MQ=255)

                                                                 cgcgCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCtt                        <  1:30325/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGGACTGTTCGCCCGTGAGCAGAGCACGCttt                       >  1:3820821/1‑93 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       <  1:656503/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       <  2:3211800/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       <  1:445634/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTACCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       >  2:2584147/1‑93 (MQ=255)

                                                                   cgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTc                      <  1:4422925/93‑1 (MQ=255)

                                                                   cgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTc                      <  1:1941879/93‑1 (MQ=255)

                                                                     cAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCgg                    >  1:123290/1‑93 (MQ=255)

                                                                      aGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGc                   >  1:251913/1‑93 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                >  2:580779/1‑93 (MQ=255)

                                                                          aaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGt               >  1:2864238/1‑93 (MQ=255)

                                                                            aCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTAc             >  2:268955/1‑93 (MQ=255)

                                                                             cGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTAcc            <  2:1427820/93‑1 (MQ=255)

                                                                               cATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCgg          <  2:3987865/93‑1 (MQ=255)

                                                                                aTCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGc         <  2:4440755/93‑1 (MQ=255)

                                                                                aTCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGc         >  1:1817829/1‑93 (MQ=255)

                                                                                aTCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGc         >  2:630885/1‑93 (MQ=255)

                                                                                  cGGCCTTGGATTGTTCGGTAAGGGCCATATCCGTTCCCGTTGGCTGGAACTGTTCTCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCtt       <  1:1543977/93‑1 (MQ=255)

                                                                                  cGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCtt       >  2:480559/1‑93 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTGGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     >  1:3197306/1‑93 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCTGCTTTg     >  2:3458310/1‑93 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     >  2:4642525/1‑93 (MQ=255)

                                                                                      cTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAg   >  2:1391952/1‑93 (MQ=255)

                                                                                      cTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAg   >  2:4315117/1‑93 (MQ=255)

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  >  1:3373555/1‑93 (MQ=255)

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

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