BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	865,377	C→T	G71R (GGG→AGG)	moeA ←	MoeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	865,377	0	C	T	96.4%	475.6 / 2.0	140	G71R (GGG→AGG)	moeA	MoeA
Reads supporting (aligned to +/- strand): ref base C (2/3); new base T (66/69); total (68/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 145 total aligned reads displayed.

ACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCT  >  NC_000913/865285‑865464

aCCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGccc                                                                                         >  2:4324068/1‑93 (MQ=255)

   ggTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGa                                                                                      >  1:3670101/1‑93 (MQ=255)

   ggTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGa                                                                                      >  1:1854116/1‑93 (MQ=255)

        taataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCaa                                                                                 <  1:488328/93‑1 (MQ=255)

          ataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAata                                                                               >  2:243124/1‑93 (MQ=255)

          ataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAata                                                                               <  2:1492622/93‑1 (MQ=255)

               aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGc                                                                          <  2:3426046/93‑1 (MQ=255)

               aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGc                                                                          <  1:4040345/93‑1 (MQ=255)

               aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTGACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGc                                                                          <  1:2064740/93‑1 (MQ=255)

                    aGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAAc                                                                     <  1:819291/93‑1 (MQ=255)

                    aGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAAc                                                                     <  1:743623/93‑1 (MQ=255)

                    aGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAAc                                                                     >  1:1833753/1‑93 (MQ=255)

                      gTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGc                                                                   >  1:1949131/1‑93 (MQ=255)

                        aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                 >  2:1969861/1‑93 (MQ=255)

                        aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                 >  1:4517678/1‑93 (MQ=255)

                         cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAcc                                                                <  2:2040125/93‑1 (MQ=255)

                         cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAcc                                                                <  1:1637255/93‑1 (MQ=255)

                          ccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcg                                                               <  2:1922454/93‑1 (MQ=255)

                           cGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgc                                                              <  1:762475/93‑1 (MQ=255)

                           cGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgc                                                              >  1:2780315/1‑93 (MQ=255)

                            gCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgcg                                                             >  1:3394893/1‑93 (MQ=255)

                             cAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGATGCCTGGAGGCAATATCGGCTAAACGCACCGCGt                                                            <  2:2975815/93‑1 (MQ=255)

                              aGGCCATTCACTATGGTATGGCTGACCGGCAATTGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                           <  1:3861472/93‑1 (MQ=255)

                              aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                           <  2:1415300/93‑1 (MQ=255)

                              aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                           <  1:1901278/93‑1 (MQ=255)

                              aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                           >  2:4541283/1‑93 (MQ=255)

                              aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                           <  2:2201402/93‑1 (MQ=255)

                               ggCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAg                                                          >  2:1084011/1‑93 (MQ=255)

                                gCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGc                                                         >  2:3851561/1‑93 (MQ=255)

                                 ccATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGcc                                                        <  2:1494819/93‑1 (MQ=255)

                                  cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                       >  1:1561630/1‑93 (MQ=255)

                                    ttCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTc                                                     >  1:993348/1‑93 (MQ=255)

                                    ttCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTc                                                     >  2:2764920/1‑93 (MQ=255)

                                    ttCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGCCGGCAATATCGGCTAAACGCACACCGTAGCCGTc                                                     >  1:993349/1‑93 (MQ=255)

                                       aCCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAt                                                  >  1:1608248/1‑93 (MQ=255)

                                       aCCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAt                                                  <  1:1479274/93‑1 (MQ=255)

                                       aCCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGACTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCTTCCAt                                                  >  2:1477669/1‑93 (MQ=255)

                                        ccATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAtt                                                 >  1:1768877/1‑93 (MQ=255)

                                           tGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCg                                              <  2:1190954/93‑1 (MQ=255)

                                            ggTATGGCTGACCGGCAAATGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCgg                                             <  1:2974973/93‑1 (MQ=255)

                                             gTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGa                                            <  2:3519515/93‑1 (MQ=255)

                                             gTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGa                                            >  1:580055/1‑93 (MQ=255)

                                               aTGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGt                                          >  1:429744/1‑93 (MQ=255)

                                               aTGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGt                                          <  2:63136/93‑1 (MQ=255)

                                                 ggCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTa                                        <  1:2818970/93‑1 (MQ=255)

                                                 ggCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTa                                        >  1:100852/1‑93 (MQ=255)

                                                  gCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCTGCTAAACGCACCGCGTAGCCGTCCATGGCGGAGTTAt                                       >  1:1501654/1‑93 (MQ=255)

                                                   gtGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                      <  2:3575778/92‑1 (MQ=255)

                                                   cTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                      >  2:3892452/1‑93 (MQ=255)

                                                   cTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                      <  2:1036752/93‑1 (MQ=255)

                                                    tGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCa                                     >  2:2337202/1‑93 (MQ=255)

                                                    tGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCa                                     >  1:3708881/1‑93 (MQ=255)

                                                        cGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcc                                 <  2:3027373/93‑1 (MQ=255)

                                                        cGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcc                                 <  2:762187/93‑1 (MQ=255)

                                                        cGGCAAACGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGTCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcc                                 <  2:3086162/93‑1 (MQ=255)

                                                          gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcccc                               <  2:3279557/93‑1 (MQ=255)

                                                          gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcccc                               <  2:3640988/93‑1 (MQ=255)

                                                          gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcccc                               <  2:3640989/93‑1 (MQ=255)

                                                            aaaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCgg                             >  1:1551935/1‑93 (MQ=255)

                                                            aaaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCgg                             >  2:4168994/1‑93 (MQ=255)

                                                            aaaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCgg                             <  1:2815264/93‑1 (MQ=255)

                                                            aaaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGGTATCACACCCCgg                             >  2:4168995/1‑93 (MQ=255)

                                                            aaaGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCCTCCATTGCGGAGTTATCAAACCCCgg                             >  1:123463/1‑93 (MQ=255)

                                                             aaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGa                            >  2:1875362/1‑93 (MQ=255)

                                                              aGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGaa                           <  2:1064841/93‑1 (MQ=255)

                                                               ggATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAAc                          <  2:1412658/93‑1 (MQ=255)

                                                                gATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACa                         <  2:1172274/93‑1 (MQ=255)

                                                                   ttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCa                      <  2:897547/93‑1 (MQ=255)

                                                                    tACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCaa                     <  2:1984224/93‑1 (MQ=255)

                                                                    tACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCaa                     <  1:674309/93‑1 (MQ=255)

                                                                      ccGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgc                   >  1:3125169/1‑93 (MQ=255)

                                                                      ccGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgc                   <  1:3456136/93‑1 (MQ=255)

                                                                      ccGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgc                   <  1:367683/93‑1 (MQ=255)

                                                                      ccGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgc                   <  2:3574556/93‑1 (MQ=255)

                                                                       cGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcg                  >  2:1866436/1‑93 (MQ=255)

                                                                        ggCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcgg                 <  1:1046253/93‑1 (MQ=255)

                                                                        ggCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcgg                 >  1:895673/1‑93 (MQ=255)

                                                                         gCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcggc                <  1:1705740/93‑1 (MQ=255)

                                                                         gCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcggc                >  1:407158/1‑93 (MQ=255)

                                                                          cAACGGGCAGCGGTTGCATGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCAATTGCGGAGTTATCAAACCCCGGAACATCAAgcggcg               <  1:2941726/93‑1 (MQ=255)

                                                                           aaCGGGCAGCGGTTGCCTGGAGTCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGa              <  1:4399368/93‑1 (MQ=255)

                                                                           aaCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGa              >  1:212659/1‑93 (MQ=255)

                                                                            aCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaa             >  1:4405665/1‑93 (MQ=255)

                                                                             cGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaaa            <  2:682369/93‑1 (MQ=255)

                                                                              gggCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAAc           >  2:154239/1‑93 (MQ=255)

                                                                              gggCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAAc           >  2:154240/1‑93 (MQ=255)

                                                                              gggCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAAc           <  2:3923387/93‑1 (MQ=255)

                                                                               ggCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACg          <  1:3582779/93‑1 (MQ=255)

                                                                                 cAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGAc        <  2:4691025/93‑1 (MQ=255)

                                                                                 cAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGAc        <  2:4691024/93‑1 (MQ=255)

                                                                                   gCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAt      <  2:2298687/93‑1 (MQ=255)

                                                                                   gCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAt      <  2:1927399/93‑1 (MQ=255)

                                                                                    cGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtc     <  1:2842240/93‑1 (MQ=255)

                                                                                    cGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtc     >  2:2827677/1‑93 (MQ=255)

                                                                                    cGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtc     <  2:3123888/93‑1 (MQ=255)

                                                                                     ggTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcg    >  1:2587003/1‑93 (MQ=255)

                                                                                      gTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgc   <  1:3063214/93‑1 (MQ=255)

                                                                                      gTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgc   >  1:1060438/1‑93 (MQ=255)

                                                                                      gTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgc   <  1:587464/93‑1 (MQ=255)

                                                                                       ttGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgct  <  1:4251359/93‑1 (MQ=255)

ACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCT  >  NC_000913/865285‑865464

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