BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	411,562	T→C	A76A (GCA→GCG)	araJ ←	AraJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	411,562	0	T	C	97.4%	459.2 / 1.3	152	A76A (GCA→GCG)	araJ	AraJ
Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (87/61); total (88/64)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 161 total aligned reads displayed.

GCACAACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAATTC  >  NC_000913/411472‑411651

 cacaACGCCACCAGAATTAACAAGATATGTTTGAGTGAGTAGCGGCTGGTAAAGAGTGCGATGATTGGCGCACCGACCACTATCCCCAGCGCa                                                                                        >  2:4212242/1‑93 (MQ=255)

 cacaACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCa                                                                                        >  1:4675157/1‑93 (MQ=255)

   caACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCata                                                                                      >  2:1887213/1‑93 (MQ=255)

     aCGCTACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGGAAAGACGGCGATGAATTGAGCAACGACCACCACCCCCAGCGCataat                                                                                    >  2:2474867/1‑93 (MQ=255)

       gCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATAc                                                                                  <  2:1415115/93‑1 (MQ=255)

       gCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATAc                                                                                  <  1:2070564/93‑1 (MQ=255)

         caccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGATCACCACCCCCAGCGCATAATACga                                                                                >  2:473219/1‑93 (MQ=255)

         caccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACga                                                                                >  2:2626009/1‑93 (MQ=255)

          accaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACgag                                                                               >  1:608837/1‑93 (MQ=255)

           ccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGGATAATACgaga                                                                              >  2:3128952/1‑93 (MQ=255)

           ccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACgaga                                                                              <  1:1445076/93‑1 (MQ=255)

            caGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGAt                                                                             >  2:3950065/1‑93 (MQ=255)

             aGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATc                                                                            <  1:3123913/93‑1 (MQ=255)

              gAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCa                                                                           <  1:3451578/93‑1 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCCCCCCCAGCGCATAATACGAGATCata                                                                         <  2:3455659/93‑1 (MQ=255)

                 acaacaAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCatat                                                                        >  2:1616711/1‑93 (MQ=255)

                 acaacaAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCatat                                                                        <  1:1088369/93‑1 (MQ=255)

                   aacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGc                                                                      >  2:2495052/1‑93 (MQ=255)

                   aacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGc                                                                      >  2:408102/1‑93 (MQ=255)

                   aacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGc                                                                      >  1:2554504/1‑93 (MQ=255)

                   aacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGc                                                                      >  1:2940862/1‑93 (MQ=255)

                    acaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGcc                                                                     <  1:1870977/93‑1 (MQ=255)

                     caaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGccc                                                                    >  2:3420578/1‑93 (MQ=255)

                     caaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGccc                                                                    <  2:3764575/93‑1 (MQ=255)

                       aGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCCCCACCCCCAGCGCATAATACGAGATCATATGCCcgg                                                                  <  2:3175827/93‑1 (MQ=255)

                        gATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGGGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATTATATGCCcggc                                                                 >  1:1818946/1‑93 (MQ=255)

                         atatGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggcg                                                                >  1:1741390/1‑93 (MQ=255)

                          tatGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggcgg                                                               >  1:4716190/1‑93 (MQ=255)

                           atGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCggcggc                                                              >  1:858347/1‑93 (MQ=255)

                           atGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCggcggc                                                              <  1:2791344/93‑1 (MQ=255)

                            tGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCa                                                             <  2:830653/93‑1 (MQ=255)

                            tGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCa                                                             <  2:3527676/93‑1 (MQ=255)

                            tGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCa                                                             >  1:778049/1‑93 (MQ=255)

                             gTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAg                                                            <  2:3562215/93‑1 (MQ=255)

                              tttGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAgg                                                           <  1:3339715/93‑1 (MQ=255)

                                tgagtgagTAGCGGCTGGAAAAGAGGGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAggga                                                         >  2:1069825/1‑91 (MQ=255)

                                  agtgagtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATc                                                       <  2:1563280/93‑1 (MQ=255)

                                        tAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAAtt                                                 >  1:1510289/1‑93 (MQ=255)

                                         aGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTc                                                >  2:4448596/1‑93 (MQ=255)

                                          gCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTcc                                               <  1:1905731/93‑1 (MQ=255)

                                           cGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCCCCCCCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCt                                              <  2:1180513/93‑1 (MQ=255)

                                            ggCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTa                                             <  2:4039383/93‑1 (MQ=255)

                                            cgCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTa                                             <  1:2131128/92‑1 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            <  2:4031750/93‑1 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            >  1:3893805/1‑93 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            >  1:4170275/1‑93 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            <  2:2230920/93‑1 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            <  1:2287856/93‑1 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            >  1:2075257/1‑93 (MQ=255)

                                             gCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTAc                                            >  2:2747337/1‑93 (MQ=255)

                                               tGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGt                                          <  2:2063973/93‑1 (MQ=255)

                                               tGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGt                                          <  2:1455950/93‑1 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                        >  1:1740078/1‑93 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                        <  1:2332589/93‑1 (MQ=255)

                                                   aaaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATg                                      >  1:2001902/1‑93 (MQ=255)

                                                    aaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgc                                     <  2:3082500/93‑1 (MQ=255)

                                                    aaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgc                                     <  1:2499504/93‑1 (MQ=255)

                                                    aaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgc                                     >  2:330007/1‑93 (MQ=255)

                                                      gagTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgcgc                                   >  2:3423577/1‑93 (MQ=255)

                                                        gTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCa                                 <  1:663776/93‑1 (MQ=255)

                                                        gTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCa                                 >  2:2278445/1‑93 (MQ=255)

                                                          gCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGc                               >  1:3194119/1‑93 (MQ=255)

                                                          gCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGc                               >  1:557371/1‑93 (MQ=255)

                                                             atgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTcc                            >  2:3584350/1‑93 (MQ=255)

                                                             atgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTcc                            <  2:1140790/93‑1 (MQ=255)

                                                              tgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCg                           <  1:3892069/93‑1 (MQ=255)

                                                              tgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCg                           >  2:2706528/1‑93 (MQ=255)

                                                              tgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCg                           <  1:1930074/93‑1 (MQ=255)

                                                               gatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGt                          >  1:3541394/1‑93 (MQ=255)

                                                               gatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGt                          <  2:1086332/93‑1 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                        >  1:2194700/1‑93 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                        <  2:4558759/93‑1 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                        <  1:1041400/93‑1 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                        >  2:29595/1‑93 (MQ=255)

                                                                   ggCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGc                      >  1:3620079/1‑93 (MQ=255)

                                                                    gcgcACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCa                     <  2:868894/93‑1 (MQ=255)

                                                                     cgcACCGCCCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCAc                    <  1:2893491/93‑1 (MQ=255)

                                                                     cgcACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCAc                    >  2:4322590/1‑93 (MQ=255)

                                                                        aCCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGcc                 >  2:1575629/1‑93 (MQ=255)

                                                                         ccGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGccc                >  2:4146240/1‑93 (MQ=255)

                                                                         ccGACCACCACCCCCAGCGCATAATACGAGATCATATGACCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCGAGCTCCGTGAGCACGccc                >  2:1881765/1‑93 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa               >  2:884290/1‑93 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa               >  1:2232885/1‑93 (MQ=255)

                                                                            accaccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATa             >  2:649290/1‑93 (MQ=255)

                                                                            accaccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATa             >  1:561282/1‑93 (MQ=255)

                                                                             ccaccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATaa            <  2:4682313/93‑1 (MQ=255)

                                                                              caccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGCGCACGCCCATAAt           >  1:3171674/1‑93 (MQ=255)

                                                                                ccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGc         >  1:313026/1‑93 (MQ=255)

                                                                                  ccccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCa       <  1:1257733/92‑1 (MQ=255)

                                                                                  accCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGCGCACGCCCATAATGCCa       >  1:1736705/1‑93 (MQ=255)

                                                                                  accCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCTCGCCCATACTGCCa       >  2:3749249/1‑93 (MQ=255)

                                                                                   cccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaa      >  2:3914316/1‑93 (MQ=255)

                                                                                    ccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaaa     <  1:1638198/93‑1 (MQ=255)

                                                                                    ccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaaa     <  1:4469880/93‑1 (MQ=255)

                                                                                    ccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaaa     >  2:1777352/1‑93 (MQ=255)

                                                                                    ccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaaa     >  1:1646428/1‑93 (MQ=255)

                                                                                     cccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAAt    <  1:3782006/93‑1 (MQ=255)

                                                                                     cccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAAt    >  1:3247871/1‑93 (MQ=255)

                                                                                     cccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGCGAGCACGCCCATAATGCCAAAt    >  1:446689/1‑93 (MQ=255)

                                                                                       cAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAATTc  <  1:2177979/93‑1 (MQ=255)

GCACAACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAATTC  >  NC_000913/411472‑411651

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