BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,055	Δ1 bp	coding (509/1458 nt)	ascF →	PTS system enzyme IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,055	0	G	.	99.7%	1344.8 / ‑3.6	292	coding (509/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (147/144); minor base C (0/1); total (147/145)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 300 total aligned reads displayed.

CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑G‑CATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGC  >  NC_000913/2837968‑2838142

cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                                    <  2:1922490/93‑1 (MQ=255)

cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                                    >  1:8697441/1‑93 (MQ=255)

 aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCaa                                                                                   >  2:2990787/1‑93 (MQ=255)

   aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAAt                                                                                 <  2:134985/93‑1 (MQ=255)

     ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAAttt                                                                               >  1:3593074/1‑93 (MQ=255)

       aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaa                                                                             >  1:1213849/1‑93 (MQ=255)

       aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaa                                                                             >  2:4002200/1‑93 (MQ=255)

       aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaa                                                                             >  1:6650269/1‑93 (MQ=255)

        tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaaa                                                                            <  2:6382270/93‑1 (MQ=255)

         cGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaaaa                                                                           <  2:6337713/93‑1 (MQ=255)

         cGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑AATCTGCCGC‑CATCAAATTTaaaa                                                                           >  1:1476074/1‑93 (MQ=255)

          gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAAc                                                                          <  2:5689629/93‑1 (MQ=255)

           ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAAcc                                                                         <  1:8382468/93‑1 (MQ=255)

            cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                                        <  2:6517026/93‑1 (MQ=255)

            cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                                        >  2:5757413/1‑93 (MQ=255)

            cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                                        <  2:2615068/93‑1 (MQ=255)

             gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCaa                                                                       >  2:592668/1‑93 (MQ=255)

             gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCaa                                                                       >  1:1735755/1‑93 (MQ=255)

             gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCaa                                                                       <  1:3768169/93‑1 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     <  1:5232351/93‑1 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     <  1:551480/93‑1 (MQ=255)

               ccTTAACCATTCTGAATGTGATTGGCGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     <  1:1066316/93‑1 (MQ=255)

               ccTTAACCATTCTGAAAGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                                     >  2:1264348/1‑93 (MQ=255)

                cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAatat                                                                    >  2:6704110/1‑93 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   <  2:3756550/93‑1 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   >  2:8014619/1‑93 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   <  2:8417059/93‑1 (MQ=255)

                 ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                                   >  2:1061272/1‑93 (MQ=255)

                  tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGt                                                                  >  1:2883290/1‑93 (MQ=255)

                   aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                                 <  2:4824836/93‑1 (MQ=255)

                   aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                                 <  2:439317/93‑1 (MQ=255)

                    aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                                >  1:4906138/1‑93 (MQ=255)

                     ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                               <  1:3116826/93‑1 (MQ=255)

                     ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                               >  2:4040382/1‑93 (MQ=255)

                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                              <  2:5473927/93‑1 (MQ=255)

                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                              >  1:7759249/1‑93 (MQ=255)

                      cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                              <  2:1570939/93‑1 (MQ=255)

                       aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                             >  1:8057779/1‑93 (MQ=255)

                        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                            >  1:3014337/1‑93 (MQ=255)

                        ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                            >  1:859163/1‑93 (MQ=255)

                        ttCTGAATGTGATTGGTGACGGTGCTTTATTCTTCCTTCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAATACCAAGATGTCGCTgg                                                            >  2:6245046/1‑93 (MQ=255)

                         tCTGAATGTTATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGc                                                           <  1:6531595/93‑1 (MQ=255)

                         tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGc                                                           <  1:5275905/93‑1 (MQ=255)

                            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                                        <  1:7495778/93‑1 (MQ=255)

                            gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                                        >  2:8753904/1‑93 (MQ=255)

                             aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                                       <  2:8227446/93‑1 (MQ=255)

                               tgtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAACCCAATATGTCGCTGGCGATTGc                                                     >  1:7348129/1‑93 (MQ=255)

                                 tgATTGGTGACGGTGTTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                                   <  2:2997419/93‑1 (MQ=255)

                                 tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                                   >  1:5339197/1‑93 (MQ=255)

                                  gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                                  >  2:5558085/1‑93 (MQ=255)

                                  gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGCCGCTGGCGATTGCCAt                                                  <  2:8363624/93‑1 (MQ=255)

                                   aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                                 >  1:2190477/1‑93 (MQ=255)

                                   aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                                 <  1:6776553/93‑1 (MQ=255)

                                    ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                                <  2:1793013/93‑1 (MQ=255)

                                      ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                                              <  1:727553/93‑1 (MQ=255)

                                      ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                                              <  1:2159668/93‑1 (MQ=255)

                                        tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                            <  2:3294497/93‑1 (MQ=255)

                                         gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                                           >  1:4339545/1‑93 (MQ=255)

                                           cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                                         <  2:1470659/93‑1 (MQ=255)

                                           cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                                         >  2:7724658/1‑93 (MQ=255)

                                           cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgg                                         >  2:7098117/1‑92 (MQ=255)

                                             gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                                       >  1:1657192/1‑93 (MQ=255)

                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                                   >  2:5991785/1‑93 (MQ=255)

                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                                   >  1:5452529/1‑93 (MQ=255)

                                                 tttCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                                   >  1:6066502/1‑93 (MQ=255)

                                                  ttcttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                                  <  1:5772252/93‑1 (MQ=255)

                                                   tcttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                                 <  2:2315585/93‑1 (MQ=255)

                                                   tcttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                                 >  2:6745498/1‑93 (MQ=255)

                                                     ttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                               >  2:3570107/1‑93 (MQ=255)

                                                     ttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                               >  1:4387316/1‑93 (MQ=255)

                                                      tcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                              <  1:2417333/93‑1 (MQ=255)

                                                      tcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                              >  1:7048621/1‑93 (MQ=255)

                                                       ctcccTGCCGCTGATGGTCGC‑GTCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                             <  2:4524235/90‑1 (MQ=255)

                                                        ttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                            <  1:3862016/93‑1 (MQ=255)

                                                         tcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                           <  2:5852168/93‑1 (MQ=255)

                                                         tcCTGCCGATGATGGCCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                           <  2:4310244/93‑1 (MQ=255)

                                                          ccTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg                          <  1:6194311/93‑1 (MQ=255)

                                                          ccTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg                          <  1:8192887/93‑1 (MQ=255)

                                                           cTGCCGCTGATGTTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt                         <  2:4986881/93‑1 (MQ=255)

                                                              ccGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt                      <  1:1899303/93‑1 (MQ=255)

                                                              ccGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt                      >  1:3785880/1‑93 (MQ=255)

                                                              ccGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt                      >  2:923378/1‑93 (MQ=255)

                                                               cGCTGGTGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa                     <  1:6912030/93‑1 (MQ=255)

                                                                gCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt                    >  2:7562668/1‑93 (MQ=255)

                                                                         cgcgg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGAt           <  2:6147599/89‑1 (MQ=255)

                                                                         cgcgg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGAt           <  2:5051293/89‑1 (MQ=255)

                                                                          gcgg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATg          >  2:3888031/4‑93 (MQ=255)

                                                                           cgg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGCTTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATgg         >  1:2696465/3‑93 (MQ=255)

                                                                            gg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGc        <  2:1711776/92‑1 (MQ=255)

                                                                            gg‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGc        <  2:275640/92‑1 (MQ=255)

                                                                             g‑cATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCg       >  1:6420728/1‑93 (MQ=255)

                                                                                aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATNCGNNTTTTATTGAGCNGATGGCNag     >  1:465812/1‑92 (MQ=255)

                                                                                aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaa     >  1:222323/1‑93 (MQ=255)

                                                                                aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaa     >  1:4809229/1‑93 (MQ=255)

                                                                                aTCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGaa     <  2:797646/93‑1 (MQ=255)

                                                                                  cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAg   <  1:5034502/93‑1 (MQ=255)

                                                                                  cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAg   >  1:8699330/1‑93 (MQ=255)

                                                                                  cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAg   <  1:415495/93‑1 (MQ=255)

                                                                                  cTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAg   <  2:21307/93‑1 (MQ=255)

                                                                                   tGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGGGTGCTGGTACATCCGAGCTTTATTGAACTGATGGCGAAAga  >  1:2313467/1‑92 (MQ=255)

CAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑G‑CATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGATGGCGAAAGC  >  NC_000913/2837968‑2838142

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