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	99.0%	701.8 / ‑4.4	199	G71R (GGG→AGG)	moeA	MoeA
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (91/105); total (93/106)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 203 total aligned reads displayed.

CGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCTCG  >  NC_000913/865287‑865466

cGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTgg                                                                                         <  1:5136195/93‑1 (MQ=255)

  gTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAg                                                                                       >  2:817524/1‑93 (MQ=255)

  gTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAg                                                                                       <  1:6125226/93‑1 (MQ=255)

  gTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAg                                                                                       >  1:3868212/1‑93 (MQ=255)

  gTCATAATACGAATGCAGGAACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAg                                                                                       >  1:4809358/1‑93 (MQ=255)

   tCATAATACGAATGCAGGTCCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGGTGCCTGGAgg                                                                                      <  2:1176804/93‑1 (MQ=255)

   tCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAgg                                                                                      >  2:6489813/1‑93 (MQ=255)

   tCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAgg                                                                                      >  2:2787918/1‑93 (MQ=255)

   gcATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAgg                                                                                      <  1:853674/92‑1 (MQ=255)

     ataataCGACTGCAGGTACCCGCAGGCCATTCCCCATGGTATGGCTGACCGGCAAAGGATTTAACGGCGACGGGCAGCGGGTGCCTGGAGGcc                                                                                    >  2:3846989/1‑92 (MQ=255)

     ataataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCa                                                                                    >  1:6455909/1‑93 (MQ=255)

      taataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGGCCGGCAAAGGATTTCCCGGCAACGGGCAGCGGTTGCCTGGAGGCaa                                                                                   <  1:4629874/93‑1 (MQ=255)

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

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

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                                >  2:2424328/1‑93 (MQ=255)

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                                <  2:226421/93‑1 (MQ=255)

           cGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCg                                                                              <  2:3181667/93‑1 (MQ=255)

             aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGc                                                                            >  2:7911492/1‑93 (MQ=255)

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

                gCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaa                                                                         >  2:374082/1‑93 (MQ=255)

                 cAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaaa                                                                        <  1:6215494/93‑1 (MQ=255)

                 cAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaaa                                                                        <  2:6172594/93‑1 (MQ=255)

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

                   ggTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACg                                                                      >  1:3147939/1‑93 (MQ=255)

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

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                   <  1:7644147/93‑1 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                   <  1:7597405/93‑1 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                   <  2:505076/93‑1 (MQ=255)

                       cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAcc                                                                  >  2:1065377/1‑93 (MQ=255)

                       cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCAcc                                                                  >  1:5530348/1‑93 (MQ=255)

                        ccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcg                                                                 >  2:901069/1‑93 (MQ=255)

                         cGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGTAATATCGGCTAAACGCACcgc                                                                <  2:2991856/93‑1 (MQ=255)

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

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

                             ggCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGTAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAAAGCACCGCGTAg                                                            >  1:3764091/1‑93 (MQ=255)

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

                                cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                         >  2:3312029/1‑93 (MQ=255)

                                cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                         <  2:3589889/93‑1 (MQ=255)

                                cATTCACCATGGTATGGCTGAACCGCAAAGGAATTACCGGGAACGGGCAGCGGTTGCCTGGAGGCAATATTAGCTAAACCCACCGCGCAGCCg                                                         >  1:5600868/1‑93 (MQ=255)

                                 aTTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGt                                                        >  1:948330/1‑93 (MQ=255)

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

                                      ccATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAtt                                                   <  1:7864321/93‑1 (MQ=255)

                                      ccATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAtt                                                   <  2:6869283/93‑1 (MQ=255)

                                        aTGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGc                                                 >  2:3632687/1‑93 (MQ=255)

                                        aTGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGc                                                 <  2:1683465/93‑1 (MQ=255)

                                         tGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCg                                                >  1:6728142/1‑93 (MQ=255)

                                         tGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCg                                                >  1:7464057/1‑93 (MQ=255)

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

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

                                            tATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAg                                             >  1:6860763/1‑93 (MQ=255)

                                            tATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGCCCATTGCGGAg                                             >  2:5366495/1‑93 (MQ=255)

                                            tATGGCAGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAg                                             >  2:7679818/1‑93 (MQ=255)

                                             aTGTCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGt                                            <  1:7304844/93‑1 (MQ=255)

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

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

                                                gCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTAt                                         <  1:293007/93‑1 (MQ=255)

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

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

                                                   gACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaa                                      <  2:5376540/93‑1 (MQ=255)

                                                    aCCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaaa                                     <  1:6079786/93‑1 (MQ=255)

                                                        gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATGAAAcccg                                 >  1:1442389/1‑92 (MQ=255)

                                                        gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcccc                                 >  2:1840589/1‑93 (MQ=255)

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

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

                                                           aaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGa                              <  1:2857683/93‑1 (MQ=255)

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

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

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

                                                           aaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGGTATCAAACCCCGGa                              >  2:4498591/1‑93 (MQ=255)

                                                             ggATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAAc                            >  2:2037417/1‑93 (MQ=255)

                                                              gATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACa                           >  1:6479681/1‑93 (MQ=255)

                                                                tttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATc                         <  2:7539552/93‑1 (MQ=255)

                                                                tttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATc                         >  1:114528/1‑93 (MQ=255)

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

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

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

                                                                   aCCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAg                      <  2:32852/93‑1 (MQ=255)

                                                                     cGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcg                    <  2:4721959/93‑1 (MQ=255)

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

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

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

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

                                                                         aaCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGCCGa                >  2:1243340/1‑93 (MQ=255)

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

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

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

                                                                           agggCAGCGGTTGCCGGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaaa              <  1:2299656/92‑1 (MQ=255)

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

                                                                             ggCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACg            >  2:231324/1‑93 (MQ=255)

                                                                               cAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGAc          >  2:5895847/1‑93 (MQ=255)

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

                                                                                aGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACa         <  2:3326643/93‑1 (MQ=255)

                                                                                aGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACa         <  1:2761703/93‑1 (MQ=255)

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

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

                                                                                  cGGTTTCCTGGAGACAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAACCCCCGGAACATCAAGCGGCGAAACGACAtc       <  2:6652634/93‑1 (MQ=255)

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

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

                                                                                      tGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgctc   <  1:1976308/93‑1 (MQ=255)

                                                                                       gCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGCAACGACAtcgctcg  >  2:442825/1‑93 (MQ=255)

CGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCTCG  >  NC_000913/865287‑865466

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