BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,327,579	T→C	Y80C (TAC→TGC)	b2226 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,327,579	0	T	C	83.4%	618.5 / 111.9	277	Y80C (TAC→TGC)	b2226	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (18/28); new base C (117/114); total (135/142)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAAT  >  NC_000913/2327487‑2327671

  ccACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAg                                                                                            <  1:6175002/93‑1 (MQ=255)

    acaGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAggg                                                                                          <  1:7049229/93‑1 (MQ=255)

      aGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgct                                                                                        >  2:2832798/1‑93 (MQ=255)

      aGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgct                                                                                        >  2:2517626/1‑93 (MQ=255)

      aGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgcg                                                                                        >  1:4426673/1‑92 (MQ=255)

       ggAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctg                                                                                       >  1:362/1‑93 (MQ=255)

        gAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGgctgc                                                                                      <  2:3214607/93‑1 (MQ=255)

        gAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGgctgc                                                                                      >  2:545048/1‑93 (MQ=255)

             tACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCgg                                                                                 >  1:1890840/1‑93 (MQ=255)

             tACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCgg                                                                                 >  2:5012604/1‑93 (MQ=255)

                cAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTg                                                                              <  1:7438684/93‑1 (MQ=255)

                cAGCGCGGTATCGCGCGTTATGCCCT‑CCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGa                                                                             <  1:4709282/93‑1 (MQ=255)

                  gcgcgcTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGAt                                                                            <  2:6788935/93‑1 (MQ=255)

                  gcgcGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGAt                                                                            <  2:1584947/93‑1 (MQ=255)

                   cgcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATg                                                                           <  1:6441586/93‑1 (MQ=255)

                   cgcgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATg                                                                           >  1:5692917/1‑93 (MQ=255)

                     cgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATggg                                                                         <  1:3534218/93‑1 (MQ=255)

                     cgGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATggg                                                                         <  2:4496375/93‑1 (MQ=255)

                       gTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTa                                                                       >  1:1309744/1‑93 (MQ=255)

                       gTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTa                                                                       <  2:5745694/93‑1 (MQ=255)

                        tATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTAt                                                                      >  2:7318805/1‑93 (MQ=255)

                         aTCGCGCGTTATGCTCTCCCCAGCCTGATAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTAtc                                                                     >  2:1076385/1‑93 (MQ=255)

                         aTCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTAtc                                                                     <  1:3810843/93‑1 (MQ=255)

                              gcgTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTcca                                                                <  1:4852053/93‑1 (MQ=255)

                                 ttATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTccagcc                                                             >  1:5861388/1‑93 (MQ=255)

                                    tGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAgg                                                          <  2:2353187/93‑1 (MQ=255)

                                     gCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGc                                                         <  1:1061622/93‑1 (MQ=255)

                                      cccaccccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCg                                                        <  1:6612158/89‑1 (MQ=255)

                                        ctcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTa                                                      >  1:6946080/1‑93 (MQ=255)

                                         tcCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTaa                                                     >  2:2263204/1‑93 (MQ=255)

                                            ccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACa                                                  >  2:6082688/1‑93 (MQ=255)

                                            ccAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGGTGATGGGGATCCTCCAGCCAGGCGTAAACa                                                  >  1:1215262/1‑93 (MQ=255)

                                             cAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACAt                                                 <  1:1402774/93‑1 (MQ=255)

                                               gCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATAc                                               <  1:1727865/93‑1 (MQ=255)

                                                ccTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACa                                              <  2:3399760/93‑1 (MQ=255)

                                                ccTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACa                                              <  1:1618001/93‑1 (MQ=255)

                                                  tGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGt                                            <  1:5117115/93‑1 (MQ=255)

                                                  tGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGt                                            >  1:729263/1‑93 (MQ=255)

                                                  tGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGt                                            <  2:6211646/93‑1 (MQ=255)

                                                   gACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGtt                                           >  1:2513751/1‑93 (MQ=255)

                                                    acaacaATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGttt                                          >  2:4421754/1‑93 (MQ=255)

                                                    acaacaATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGttt                                          <  2:4084919/93‑1 (MQ=255)

                                                      aacaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTcc                                        <  2:5366199/93‑1 (MQ=255)

                                                      aacaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTcc                                        >  1:2715210/1‑93 (MQ=255)

                                                      aacaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTcc                                        >  1:839038/1‑93 (MQ=255)

                                                       acaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCca                                       <  1:2230358/93‑1 (MQ=255)

                                                       acaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCca                                       <  2:5411275/93‑1 (MQ=255)

                                                       acaaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCca                                       <  2:6561658/93‑1 (MQ=255)

                                                        caaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcac                                      >  2:2634238/1‑93 (MQ=255)

                                                        caaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcac                                      >  1:1135052/1‑93 (MQ=255)

                                                         aaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcaca                                     >  2:4244370/1‑93 (MQ=255)

                                                         aaTAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCcaca                                     >  2:5547401/1‑93 (MQ=255)

                                                           tAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACaaa                                   <  2:1888219/93‑1 (MQ=255)

                                                           tAAACCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACaaa                                   <  1:4782724/93‑1 (MQ=255)

                                                            aaaCCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAAc                                  >  1:3191982/1‑93 (MQ=255)

                                                            aaaCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAAc                                  <  1:945082/93‑1 (MQ=255)

                                                              aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                                >  1:6031373/1‑93 (MQ=255)

                                                              aCCTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTc                                <  2:7160485/93‑1 (MQ=255)

                                                               ccTCTTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTcc                               >  1:3053358/1‑93 (MQ=255)

                                                                 tctTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTCTCCACAAACTCCCt                             <  2:6453370/93‑1 (MQ=255)

                                                                 tctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCt                             >  1:7180941/1‑93 (MQ=255)

                                                                 tctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCt                             <  2:691587/93‑1 (MQ=255)

                                                                 tctGCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCt                             <  1:6339563/93‑1 (MQ=255)

                                                                  ctTCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTa                            <  2:3980399/93‑1 (MQ=255)

                                                                   ttCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTaa                           <  2:132264/93‑1 (MQ=255)

                                                                   ttCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTaa                           >  2:3373477/1‑93 (MQ=255)

                                                                    tCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAAc                          <  2:1412476/93‑1 (MQ=255)

                                                                    tCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAAc                          <  1:1038914/93‑1 (MQ=255)

                                                                    tCCGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAAc                          <  1:6719434/93‑1 (MQ=255)

                                                                     ccGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACg                         <  2:1780452/93‑1 (MQ=255)

                                                                      cGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACgg                        <  1:609894/93‑1 (MQ=255)

                                                                      cGGGGAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACgg                        >  1:2059317/1‑93 (MQ=255)

                                                                       ggggAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGc                       >  1:5296669/1‑93 (MQ=255)

                                                                       gcggAGTTTCCGCGACACTGGTAGGGCCGCTCCGGTTGATGGGTGTCCTTCAGCCATGCGTAAACAAACAGTTTCCACAACCTCCCTAACGGc                       >  1:5611639/3‑93 (MQ=255)

                                                                         ggAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                     <  1:4211273/93‑1 (MQ=255)

                                                                         ggAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                     <  2:6876336/93‑1 (MQ=255)

                                                                         ggAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgt                     >  2:4706110/1‑93 (MQ=255)

                                                                          gAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtg                    >  2:4637156/1‑93 (MQ=255)

                                                                          gAGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtg                    >  2:6051413/1‑93 (MQ=255)

                                                                           aGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtgg                   >  2:2102981/1‑93 (MQ=255)

                                                                           aGTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtgg                   >  2:4839071/1‑93 (MQ=255)

                                                                            gTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtggt                  <  2:857044/93‑1 (MQ=255)

                                                                            gTTTCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCgtggt                  <  2:758366/93‑1 (MQ=255)

                                                                              ttCCGCGACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGa                <  2:3660331/93‑1 (MQ=255)

                                                                                 cgcgACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCg             <  2:7433137/93‑1 (MQ=255)

                                                                                  gcgACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGa            >  1:7009918/1‑93 (MQ=255)

                                                                                  gcgACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGa            >  2:3897232/1‑93 (MQ=255)

                                                                                  gcgACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGa            >  1:3291311/1‑93 (MQ=255)

                                                                                    gACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt          <  1:4282182/93‑1 (MQ=255)

                                                                                    gACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt          >  2:3541082/1‑93 (MQ=255)

                                                                                    gACACTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt          <  2:2997221/93‑1 (MQ=255)

                                                                                    gACACTGGCAGGGCTGCTCCGGATGATGGGTATCCTCCAGGCAGGCGTAAACATACAATTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGt          >  1:2913760/1‑93 (MQ=255)

                                                                                     acacTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTc         >  2:6996562/1‑93 (MQ=255)

                                                                                     acacTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAAGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTc         >  1:1339224/1‑93 (MQ=255)

                                                                                      cacTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCa        >  2:6365881/1‑93 (MQ=255)

                                                                                       acTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAg       <  2:4389074/93‑1 (MQ=255)

                                                                                       acTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAg       >  2:5837667/1‑93 (MQ=255)

                                                                                        cTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAgg      >  2:6601697/1‑93 (MQ=255)

                                                                                        cTGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAgg      <  2:7099270/93‑1 (MQ=255)

                                                                                         tGGCAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGt     >  1:4972552/1‑93 (MQ=255)

GCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTCCGCGACACTGGTAGGGCTGCTCCGGTTGATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAAT  >  NC_000913/2327487‑2327671

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