BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,434,055	A→G	P404P (CCA→CCG)	sun →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,434,055	0	A	G	99.8%	1342.8 / ‑5.1	413	P404P (CCA→CCG)	sun	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (195/216); total (195/218)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

Only 100 of 428 total aligned reads displayed.

CTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGG  >  NC_000913/3433965‑3434144

cTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCgga                                                                                         <  1:920219/93‑1 (MQ=255)

cTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCgga                                                                                         <  1:8968122/93‑1 (MQ=255)

cTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCgga                                                                                         <  1:6367033/93‑1 (MQ=255)

 tGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCggag                                                                                        <  2:9728311/93‑1 (MQ=255)

  gtccGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGc                                                                                       <  2:1356175/90‑1 (MQ=255)

  gTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGc                                                                                       >  1:9373718/1‑93 (MQ=255)

  gTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGc                                                                                       <  2:9158614/93‑1 (MQ=255)

   ttCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCa                                                                                      >  1:2632866/1‑93 (MQ=255)

   ttCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCa                                                                                      >  2:5469100/1‑93 (MQ=255)

   ttCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCa                                                                                      >  2:5895357/1‑93 (MQ=255)

    tCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGCACACCGGAGCaa                                                                                     >  2:3931229/1‑93 (MQ=255)

    tCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCaa                                                                                     <  1:65743/93‑1 (MQ=255)

     cGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAAc                                                                                    <  2:3694576/93‑1 (MQ=255)

       gtgtTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCg                                                                                  >  1:1045209/1‑93 (MQ=255)

         gtTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCggg                                                                                <  1:3035503/93‑1 (MQ=255)

          ttACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGt                                                                               >  2:1106608/1‑93 (MQ=255)

          ttACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGt                                                                               >  2:8984059/1‑93 (MQ=255)

           tACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTa                                                                              <  2:3288798/93‑1 (MQ=255)

               ggAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACa                                                                          >  1:7245812/1‑93 (MQ=255)

                gAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACaa                                                                         >  1:4871588/1‑93 (MQ=255)

                  agagAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACaaaa                                                                       >  2:6239797/1‑93 (MQ=255)

                   gagaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAAt                                                                      <  2:7607177/93‑1 (MQ=255)

                   gagaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAAt                                                                      <  2:3420487/93‑1 (MQ=255)

                      aaTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTa                                                                   >  1:1621390/1‑93 (MQ=255)

                       aTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAc                                                                  <  1:2719123/93‑1 (MQ=255)

                        tAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAcc                                                                 <  2:4394855/93‑1 (MQ=255)

                        tAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAcc                                                                 <  1:7668582/93‑1 (MQ=255)

                         aGCCTGCAGATTAAAGCCTTTTTGCAACGTCCCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCt                                                                <  2:2141674/93‑1 (MQ=255)

                          gCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTg                                                               >  1:8563200/1‑93 (MQ=255)

                           ccTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTgg                                                              >  2:1514138/1‑93 (MQ=255)

                            ctccAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGt                                                             <  1:8646540/90‑1 (MQ=255)

                            cTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGt                                                             <  2:502652/93‑1 (MQ=255)

                            cTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGt                                                             >  2:5320070/1‑93 (MQ=255)

                             tGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTg                                                            >  1:823080/1‑93 (MQ=255)

                                 gATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGa                                                        >  1:7453504/1‑93 (MQ=255)

                                  aTTAAAGCCTTTTTGCAACGTACCGCTGATGCCGACCTTTGCGAAACAGGAACGCCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGaa                                                       <  1:2262846/93‑1 (MQ=255)

                                    tAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAaga                                                     <  1:2037654/93‑1 (MQ=255)

                                    tAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAaga                                                     <  2:8191359/93‑1 (MQ=255)

                                    tAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGTTAAACAAAATCTACCTGGTGCCGAaga                                                     <  1:7080679/93‑1 (MQ=255)

                                     aaaGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAagag                                                    <  1:8038310/93‑1 (MQ=255)

                                     aaaGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAagag                                                    <  1:8323046/93‑1 (MQ=255)

                                       aGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAggg                                                  <  2:1263390/93‑1 (MQ=255)

                                        gCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGc                                                 >  1:427766/1‑93 (MQ=255)

                                        accTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGc                                                 <  1:3556525/92‑1 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAc                                              >  1:5878575/1‑93 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAc                                              >  1:1496251/1‑93 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAc                                              >  2:4494584/1‑93 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAAAAGGAACACCGGAGAACCGGGGGAAACAAATTCTACCTGGTGCCGAAGAGGgagag                                              >  1:2483185/1‑92 (MQ=255)

                                            ttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACg                                             <  2:2541470/93‑1 (MQ=255)

                                            ttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACg                                             <  1:4461162/93‑1 (MQ=255)

                                               tGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAACCAAAATCTACCTGGTGCCGAAGAGGGCGACGGct                                          <  2:3538229/93‑1 (MQ=255)

                                                gCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGctt                                         <  2:4192310/93‑1 (MQ=255)

                                                  aaCGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGcttct                                       >  1:1360656/1‑93 (MQ=255)

                                                  aaCGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGcttct                                       >  2:6946378/1‑93 (MQ=255)

                                                    cGTCCCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGTGTAAACAAAATCTACCTGGTGTCGAAGAGGGCGACGGCTTCttt                                     <  1:6034503/93‑1 (MQ=255)

                                                     gTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCtttt                                    <  2:3029251/93‑1 (MQ=255)

                                                      tACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTa                                   <  2:56734/93‑1 (MQ=255)

                                                      tACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTa                                   <  2:7818533/93‑1 (MQ=255)

                                                      tACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTa                                   >  1:9646800/1‑93 (MQ=255)

                                                      tACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTa                                   <  2:8792819/93‑1 (MQ=255)

                                                       aCCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAc                                  <  2:1963450/93‑1 (MQ=255)

                                                        ccGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACg                                 >  2:4715483/1‑93 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                >  2:2787012/1‑93 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                <  1:5547758/93‑1 (MQ=255)

                                                           cTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAACATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAAgcta                              >  2:1858856/1‑93 (MQ=255)

                                                            tGATGCCGAACTTTGCGTAGCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaa                             <  2:8955016/93‑1 (MQ=255)

                                                            tGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaa                             <  2:8341061/93‑1 (MQ=255)

                                                            tGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaa                             <  1:6900311/93‑1 (MQ=255)

                                                            tGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaa                             <  1:2441655/93‑1 (MQ=255)

                                                             gATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaag                            <  1:7737593/93‑1 (MQ=255)

                                                             gATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaag                            >  2:796354/1‑93 (MQ=255)

                                                             gATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaag                            >  1:1805728/1‑93 (MQ=255)

                                                              aTGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagc                           >  2:5763044/1‑93 (MQ=255)

                                                                gCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagcta                         <  1:4680706/93‑1 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        <  1:6433503/93‑1 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        <  1:3406374/93‑1 (MQ=255)

                                                                  cGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAAt                       <  2:6478240/93‑1 (MQ=255)

                                                                   gAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATc                      <  1:1699311/93‑1 (MQ=255)

                                                                    aaCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCa                     >  1:5177055/1‑93 (MQ=255)

                                                                     aCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                    <  1:5404602/93‑1 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  <  1:4030781/93‑1 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  >  2:4269023/1‑93 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  >  2:7328686/1‑93 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  <  2:6621166/93‑1 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  <  1:3760866/93‑1 (MQ=255)

                                                                        ttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaaa                 <  2:9459955/93‑1 (MQ=255)

                                                                            gAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtga             <  2:4825768/93‑1 (MQ=255)

                                                                             aaaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgat            <  1:9917595/93‑1 (MQ=255)

                                                                              aaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatg           <  1:3990328/93‑1 (MQ=255)

                                                                              aaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatg           <  2:8688719/93‑1 (MQ=255)

                                                                               aCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatga          <  1:2897015/93‑1 (MQ=255)

                                                                                cAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATgag         >  2:8613776/1‑93 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       >  1:9464938/1‑93 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       <  2:5584716/93‑1 (MQ=255)

                                                                                    aaCACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATaa     >  1:9823098/1‑93 (MQ=255)

                                                                                     acacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAAc    >  2:7281206/1‑93 (MQ=255)

                                                                                      cacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACg   <  2:8876856/93‑1 (MQ=255)

                                                                                      cacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACg   <  1:2112464/93‑1 (MQ=255)

                                                                                       acCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACgg  >  1:8791589/1‑93 (MQ=255)

                                                                                       acCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACgg  <  2:8033685/93‑1 (MQ=255)

CTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGG  >  NC_000913/3433965‑3434144

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