BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,383	T→C	S475S (AGT→AGC)	ytfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,383	0	T	C	97.4%	479.4 / ‑0.3	156	S475S (AGT→AGC)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (64/88); total (66/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 161 total aligned reads displayed.

TGTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGTATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTGG  >  NC_000913/4441293‑4441471

tGTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCAt                                                                                        <  1:3158091/93‑1 (MQ=255)

 gTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCAtt                                                                                       <  2:1521421/93‑1 (MQ=255)

  ttACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTc                                                                                      >  1:576614/1‑93 (MQ=255)

     cacGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                   >  1:4318293/1‑93 (MQ=255)

     cacGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTg                                                                                   >  1:2727068/1‑93 (MQ=255)

      acGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTgg                                                                                  <  2:2588354/93‑1 (MQ=255)

      acGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTgg                                                                                  >  2:2091031/1‑93 (MQ=255)

       cgcgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGc                                                                                 <  1:3314033/93‑1 (MQ=255)

        gcgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCt                                                                                <  2:2635710/93‑1 (MQ=255)

        gcgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCt                                                                                <  1:3123327/93‑1 (MQ=255)

         cgGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCta                                                                               >  2:553387/1‑93 (MQ=255)

            cACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaa                                                                            >  1:2212848/1‑93 (MQ=255)

              cGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaat                                                                          <  2:4010838/93‑1 (MQ=255)

               gctgggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaata                                                                         <  2:3658483/93‑1 (MQ=255)

               gctgggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaata                                                                         <  1:2440518/93‑1 (MQ=255)

               cctgggctggATTGAAACCCGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaata                                                                         <  2:3658484/92‑1 (MQ=255)

                 tgggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaataca                                                                       <  2:281612/93‑1 (MQ=255)

                   ggctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCtacaaatacaaa                                                                     <  1:2679721/93‑1 (MQ=255)

                     ctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATc                                                                   >  2:619526/1‑93 (MQ=255)

                     ctggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCCTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATc                                                                   <  2:150924/93‑1 (MQ=255)

                       ggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTa                                                                 <  2:438591/93‑1 (MQ=255)

                       ggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTa                                                                 >  1:2574961/1‑93 (MQ=255)

                       ggATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTa                                                                 >  2:1074345/1‑93 (MQ=255)

                        gATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTAt                                                                >  1:1722047/1‑93 (MQ=255)

                        gATGGAAACCCGTCATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATTCAAATCTAt                                                                <  2:2507633/93‑1 (MQ=255)

                         aTTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATc                                                               <  1:4589455/93‑1 (MQ=255)

                         aTTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATc                                                               >  2:144861/1‑93 (MQ=255)

                                   gTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAAt                                                     >  1:2067080/1‑93 (MQ=255)

                                     gATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATAc                                                   >  1:903905/1‑93 (MQ=255)

                                     gATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATAc                                                   <  1:3608201/93‑1 (MQ=255)

                                      aTTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACg                                                  <  2:1718235/93‑1 (MQ=255)

                                       tttCGACAAAGTACCGCCGGATGTGCGGTTCTTCGCCGGGGGCGCCCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGc                                                 <  2:598127/93‑1 (MQ=255)

                                       tttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGTGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGc                                                 <  1:2482444/93‑1 (MQ=255)

                                       tttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGAGACCGCAGCATTCGTAGCTCCAAATTCTAATCTACCGCTCCGACAttatc                                                 >  2:2154924/1‑89 (MQ=255)

                                        ttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGcc                                                <  2:3166321/93‑1 (MQ=255)

                                        ttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGcc                                                <  2:2108623/93‑1 (MQ=255)

                                        ttCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGcc                                                <  2:3896958/93‑1 (MQ=255)

                                         tCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCa                                               <  2:3025304/93‑1 (MQ=255)

                                         tCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCa                                               >  1:4167050/1‑93 (MQ=255)

                                          cGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCaa                                              <  2:469495/93‑1 (MQ=255)

                                          cGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCaa                                              >  2:2709132/1‑93 (MQ=255)

                                           gACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAAc                                             <  2:1057863/93‑1 (MQ=255)

                                           gACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAAc                                             >  2:2686345/1‑93 (MQ=255)

                                            aCAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACg                                            <  2:1964866/93‑1 (MQ=255)

                                              aaaGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGt                                          >  2:2348999/1‑93 (MQ=255)

                                              aaaGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGt                                          >  1:4601139/1‑93 (MQ=255)

                                               aaGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTg                                         <  1:4434163/93‑1 (MQ=255)

                                                aGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGa                                        >  1:969026/1‑93 (MQ=255)

                                                 gTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGAc                                       >  1:1582342/1‑93 (MQ=255)

                                                  tACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGAcc                                      >  2:2180290/1‑93 (MQ=255)

                                                   aCCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCt                                     <  1:3791396/93‑1 (MQ=255)

                                                   aCCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCt                                     >  1:1472462/1‑93 (MQ=255)

                                                    ccgccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTg                                    <  2:4530531/93‑1 (MQ=255)

                                                    ccgccgGATCTGCGTTTCTACGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCAAACGGTGACCTg                                    <  2:2156773/93‑1 (MQ=255)

                                                     cgccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGa                                   >  1:4029156/1‑93 (MQ=255)

                                                       ccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGaaa                                 >  1:817976/1‑93 (MQ=255)

                                                       ccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGaaa                                 >  1:2555929/1‑93 (MQ=255)

                                                       ccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGaaa                                 <  1:1943563/93‑1 (MQ=255)

                                                       ccgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGaaa                                 <  1:166653/93‑1 (MQ=255)

                                                        cgGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAg                                >  1:170314/1‑93 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                               <  2:3287655/93‑1 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                               >  2:1545455/1‑93 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgg                               <  2:804680/93‑1 (MQ=255)

                                                         ggATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCACTCCGAAATACGCCAACGGTGACCTGAAAgg                               >  2:1545456/1‑93 (MQ=255)

                                                          gATCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAggg                              >  2:784517/1‑93 (MQ=255)

                                                          gATCTGCGTTTCTGCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAggg                              <  1:3335999/93‑1 (MQ=255)

                                                           aTCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAgggg                             >  1:864962/1‑93 (MQ=255)

                                                            tCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGc                            >  2:941331/1‑93 (MQ=255)

                                                            tCTGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGc                            >  1:2235471/1‑93 (MQ=255)

                                                              tGCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCt                          >  2:3355903/1‑93 (MQ=255)

                                                               gCGTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTc                         <  1:2172698/93‑1 (MQ=255)

                                                                 gTTTCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGa                       <  1:2119724/93‑1 (MQ=255)

                                                                  tttCTTCGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGaa                      >  1:612307/1‑93 (MQ=255)

                                                                     tttcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGtt                   <  1:1442335/92‑1 (MQ=255)

                                                                      ttcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTg                  >  2:1824773/1‑93 (MQ=255)

                                                                      ttcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGAGGCCTCGAAGTTg                  >  1:3409506/1‑93 (MQ=255)

                                                                       tcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGa                 <  2:2502180/93‑1 (MQ=255)

                                                                       tcGCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGa                 <  1:2067121/93‑1 (MQ=255)

                                                                         gCCGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCTACGGTGACCTGAAAGGGGCCTCGAAGTTGATa               <  2:3942837/93‑1 (MQ=255)

                                                                          ccGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCTAAGCTCATaa              >  1:4091659/1‑93 (MQ=255)

                                                                          ccGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATaa              <  1:3880207/93‑1 (MQ=255)

                                                                           cGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAAc             <  1:3573463/93‑1 (MQ=255)

                                                                           cGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAAc             <  1:1187296/93‑1 (MQ=255)

                                                                           cGGGGGCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAAc             >  1:2267015/1‑93 (MQ=255)

                                                                            gggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAAcc            <  1:2418659/93‑1 (MQ=255)

                                                                             ggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCg           <  2:1322645/93‑1 (MQ=255)

                                                                             ggggCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCg           >  2:3561105/1‑93 (MQ=255)

                                                                               gacGACCGCAGCATTCGTGGCTACAAGTCCAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGa         <  1:4161656/91‑1 (MQ=255)

                                                                                gCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGAt        <  1:308837/93‑1 (MQ=255)

                                                                                gCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGAt        <  1:1795813/93‑1 (MQ=255)

                                                                                gCGACCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGAt        >  2:3447964/1‑93 (MQ=255)

                                                                                 cGACCGCAGTATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATc       <  2:4584379/93‑1 (MQ=255)

                                                                                   gccGGAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc     <  1:548631/92‑1 (MQ=255)

                                                                                   aCCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc     >  1:2661153/1‑93 (MQ=255)

                                                                                   aCCGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGc     <  2:1355401/93‑1 (MQ=255)

                                                                                     cGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTg   >  2:2839013/1‑93 (MQ=255)

                                                                                     cGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTg   >  2:3802181/1‑93 (MQ=255)

                                                                                     cGCAGCATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTg   <  2:18523/93‑1 (MQ=255)

                                                                                      gcagcaTTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTgg  >  1:643189/1‑93 (MQ=255)

                                                                                      gcagcaTTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTgg  >  1:1035585/1‑93 (MQ=255)

TGTTACACGCGGCACGCTGGGCTGGATTGAAACCGGTGATTTCGACAAAGTACCGCCGGATCTGCGTTTCTTCGCCGGGGGCGACCGCAGTATTCGTGGCTACAAATACAAATCTATCGCTCCGAAATACGCCAACGGTGACCTGAAAGGGGCCTCGAAGTTGATAACCGGATCGCTGG  >  NC_000913/4441293‑4441471

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