BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,226,344	C→T	A118T (GCC→ACC)	ycgK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,226,344	0	C	T	99.5%	730.5 / NA	206	A118T (GCC→ACC)	ycgK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (110/95); minor base A (1/0); total (111/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATC  >  NC_000913/1226255‑1226433

ttAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGtgtc                                                                                        <  1:4656431/93‑1 (MQ=255)

  aaCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGtgtcgt                                                                                      <  1:1289123/93‑1 (MQ=255)

  aaCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGtgtcgt                                                                                      >  2:3613494/1‑93 (MQ=255)

   aCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGtt                                                                                     >  1:3003708/1‑93 (MQ=255)

   aCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGtt                                                                                     <  1:235767/93‑1 (MQ=255)

   aCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGCTGTCGtt                                                                                     >  1:3378273/1‑93 (MQ=255)

    cGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATCTTTTTTAGTTTTGTTTTTACGGGTGTCGttt                                                                                    >  2:6644538/1‑93 (MQ=255)

     gTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTc                                                                                   >  1:853077/1‑93 (MQ=255)

          tGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGtt                                                                              >  1:5697099/1‑93 (MQ=255)

          tGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGtt                                                                              <  1:445749/93‑1 (MQ=255)

            aCGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTg                                                                            >  1:3686932/1‑93 (MQ=255)

            aCGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTg                                                                            <  1:3266092/93‑1 (MQ=255)

            aCGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTg                                                                            >  2:6118176/1‑93 (MQ=255)

             cGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGc                                                                           >  1:6064151/1‑93 (MQ=255)

               gggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCaa                                                                         <  2:5457368/93‑1 (MQ=255)

                  agagTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAacac                                                                      <  2:6420828/93‑1 (MQ=255)

                   gagTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACt                                                                     <  2:76184/93‑1 (MQ=255)

                   gagTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACt                                                                     <  2:2203861/93‑1 (MQ=255)

                    agTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTc                                                                    <  1:3904358/93‑1 (MQ=255)

                     gTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCg                                                                   >  1:1957081/1‑93 (MQ=255)

                       ccccGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCa                                                                 >  1:2840738/1‑93 (MQ=255)

                        cccGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAg                                                                >  2:915102/1‑93 (MQ=255)

                        cccGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAg                                                                >  1:4270226/1‑93 (MQ=255)

                          cGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGGAGCt                                                              >  1:3775461/1‑93 (MQ=255)

                          cGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCt                                                              >  1:134683/1‑93 (MQ=255)

                           gTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTc                                                             <  1:6442953/93‑1 (MQ=255)

                            tCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCa                                                            <  2:4176710/93‑1 (MQ=255)

                            tCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGGTTCGGGTTTGCAACACTCGCAGCTCa                                                            >  2:3142463/1‑93 (MQ=255)

                             cGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCat                                                           <  1:1776262/93‑1 (MQ=255)

                             cGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCat                                                           >  2:6444385/1‑93 (MQ=255)

                              ggTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCata                                                          <  1:4988029/93‑1 (MQ=255)

                                 tGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATAttt                                                       >  2:1722087/1‑93 (MQ=255)

                                  gACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTg                                                      >  2:1604850/1‑93 (MQ=255)

                                   aCATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGc                                                     >  1:2083630/1‑93 (MQ=255)

                                    cATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGcc                                                    <  1:5866779/93‑1 (MQ=255)

                                    cATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGcc                                                    <  2:6268869/93‑1 (MQ=255)

                                     aTTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCa                                                   >  2:1091715/1‑93 (MQ=255)

                                     aTTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCa                                                   <  1:307339/93‑1 (MQ=255)

                                      tttACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAg                                                  >  2:5003451/1‑93 (MQ=255)

                                      tttACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAg                                                  >  2:1932713/1‑93 (MQ=255)

                                        tACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGaa                                                <  2:5718580/93‑1 (MQ=255)

                                        tACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGaa                                                <  2:5362213/93‑1 (MQ=255)

                                         aCTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCTAGAAg                                               >  1:764269/1‑93 (MQ=255)

                                         aCTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAg                                               >  1:4394304/1‑93 (MQ=255)

                                         aCTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAg                                               >  2:3912663/1‑93 (MQ=255)

                                         aCTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAg                                               >  2:342010/1‑93 (MQ=255)

                                           ttAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCa                                             >  1:5576079/1‑93 (MQ=255)

                                            tAATTTGAATATCGACGTTATATTTTTTATTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAg                                            <  1:2921127/93‑1 (MQ=255)

                                            tAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAg                                            >  1:1823833/1‑93 (MQ=255)

                                            tAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAg                                            >  1:4342046/1‑93 (MQ=255)

                                            tAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAg                                            >  2:2871564/1‑93 (MQ=255)

                                             aaTTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAgg                                           >  1:3055036/1‑93 (MQ=255)

                                             aaTTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAgg                                           >  1:3408913/1‑93 (MQ=255)

                                             aaTTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAgg                                           >  2:1144122/1‑93 (MQ=255)

                                              aTTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAggg                                          >  1:6556376/1‑93 (MQ=255)

                                               tttGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGga                                         >  2:4345101/1‑93 (MQ=255)

                                               tttGAATAGCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAgggg                                         >  1:6568322/1‑92 (MQ=255)

                                                 tGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagt                                       <  2:64868/93‑1 (MQ=255)

                                                    atatCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtgag                                    >  2:6659893/1‑93 (MQ=255)

                                                     tatCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgcgtgagt                                   >  2:4421545/1‑93 (MQ=255)

                                                     tatCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtgagt                                   >  2:628062/1‑93 (MQ=255)

                                                      atCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTa                                  <  2:1282059/93‑1 (MQ=255)

                                                       tCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGGAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGCGAGTAt                                 >  2:4770721/1‑93 (MQ=255)

                                                       tCGACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGGTTGCAACACTCGCAGCTCATATTTACCAGAAGCAGGGAGTGAGTAt                                 >  1:4814377/1‑93 (MQ=255)

                                                         gACGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTg                               >  1:3532426/1‑93 (MQ=255)

                                                          aCGTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGt                              >  2:6290434/1‑93 (MQ=255)

                                                            gTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTcc                            >  1:3183775/1‑93 (MQ=255)

                                                            gTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTcc                            <  1:1691183/93‑1 (MQ=255)

                                                            gTTATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGGTTGCAACACTCGCAGCTCATATTTGACAGAAGCAGGGAGTGAGTATTGtgc                            >  1:6125448/1‑91 (MQ=255)

                                                             ttATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCg                           >  2:5809252/1‑93 (MQ=255)

                                                             ttATATTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCg                           >  1:1222536/1‑93 (MQ=255)

                                                             ttATATTTTTTAGTTTTGTTTTTACGGGGGTCGGTTCGGGTTTGCGACACTCGCAGCTCATATTTGGCAGAACCAGGGAGTGAGTATTGTcct                           >  2:6039724/1‑92 (MQ=255)

                                                                 atTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc                       <  2:2676288/93‑1 (MQ=255)

                                                                 atTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc                       <  2:4648500/93‑1 (MQ=255)

                                                                 atTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc                       <  2:1489894/93‑1 (MQ=255)

                                                                 atTTTTTAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc                       >  2:5657773/1‑93 (MQ=255)

                                                                    ttttAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTAt                    <  1:1013793/93‑1 (MQ=255)

                                                                     tttAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATc                   <  1:5391239/93‑1 (MQ=255)

                                                                     tttAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATc                   <  2:1330918/93‑1 (MQ=255)

                                                                     tttAGTTTTGTTTTTACGGGAGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATc                   >  1:1625731/1‑93 (MQ=255)

                                                                      ttAGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCg                  >  1:2029874/1‑93 (MQ=255)

                                                                       tAGTTTTGTTTTTACGGGTGTCGTTTCGGATTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGAGACTATCGa                 <  2:3285443/93‑1 (MQ=255)

                                                                        aGTTTTGTTTTTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAg                >  2:5168563/1‑93 (MQ=255)

                                                                           tttgttttTTCGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTc             <  2:4818605/93‑1 (MQ=255)

                                                                           tttgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTc             <  2:5622689/93‑1 (MQ=255)

                                                                           tttgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTc             <  2:5589444/93‑1 (MQ=255)

                                                                            ttgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCt            >  1:2581156/1‑93 (MQ=255)

                                                                            ttgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCt            <  2:448738/93‑1 (MQ=255)

                                                                             tgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTg           <  2:1624181/93‑1 (MQ=255)

                                                                             tgttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTg           <  1:669161/93‑1 (MQ=255)

                                                                              gttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTgg          <  1:2510264/93‑1 (MQ=255)

                                                                              gttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTgg          <  2:3785887/93‑1 (MQ=255)

                                                                              gttttTACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTgg          <  1:5490963/93‑1 (MQ=255)

                                                                                  ttACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGaa      >  2:4583909/1‑93 (MQ=255)

                                                                                  ttACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGaa      >  2:722495/1‑93 (MQ=255)

                                                                                   tACGGGTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAat     <  2:5258270/93‑1 (MQ=255)

                                                                                    aCGGGTGTCGTTTCGGGTTTGCAACGCTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAata    <  1:3138851/93‑1 (MQ=255)

                                                                                     cGGGTGGCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAatat   <  1:976150/93‑1 (MQ=255)

                                                                                      gggTGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATc  <  2:5189253/93‑1 (MQ=255)

                                                                                      gggTGTCGTTGCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATc  <  1:4050572/93‑1 (MQ=255)

TTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATC  >  NC_000913/1226255‑1226433

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