BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,011	C→T	A199A (GCC→GCT)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,011	0	C	T	99.4%	1127.4 / ‑4.7	312	A199A (GCC→GCT)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (134/176); total (136/176)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 320 total aligned reads displayed.

AGGGCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCCGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTG  >  NC_000913/4127921‑4128100

 gggCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGt                                                                                        <  1:6088230/93‑1 (MQ=255)

   gCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGa                                                                                      >  2:3593675/1‑93 (MQ=255)

   gCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGa                                                                                      <  1:3058482/93‑1 (MQ=255)

     tttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGTTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    <  2:6082709/93‑1 (MQ=255)

     tttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    <  2:3006341/93‑1 (MQ=255)

     tttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    >  2:4664507/1‑93 (MQ=255)

     tttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    >  1:7065835/1‑93 (MQ=255)

     tttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    >  1:705225/1‑93 (MQ=255)

      ttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAAc                                                                                   <  2:5505422/93‑1 (MQ=255)

       tCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACg                                                                                  <  2:2689144/93‑1 (MQ=255)

        cTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaaa                                                                                    >  1:6779722/1‑90 (MQ=255)

        cTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACgg                                                                                 <  1:6957131/93‑1 (MQ=255)

        cTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACgg                                                                                 <  2:3183249/93‑1 (MQ=255)

         ttACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGt                                                                                >  1:2243206/1‑93 (MQ=255)

           aCCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgc                                                                              <  2:5357655/93‑1 (MQ=255)

             ccGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgctg                                                                            >  1:7450740/1‑93 (MQ=255)

              cGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgctgc                                                                           <  1:1640342/93‑1 (MQ=255)

              cGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgctgc                                                                           >  1:1048337/1‑93 (MQ=255)

               gTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTgctgct                                                                          >  1:1994711/1‑93 (MQ=255)

                gtgctgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGctgctg                                                                         <  1:234252/92‑1 (MQ=255)

                 tgctgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTgg                                                                        <  1:6532219/93‑1 (MQ=255)

                     gcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCg                                                                    >  2:5852223/1‑93 (MQ=255)

                     gcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCg                                                                    >  1:7848882/1‑93 (MQ=255)

                      cAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGt                                                                   <  2:7704135/93‑1 (MQ=255)

                        aCAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTaa                                                                 >  2:2548870/1‑93 (MQ=255)

                        aCAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTaa                                                                 >  1:3413610/1‑93 (MQ=255)

                        aCAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTaa                                                                 >  1:724731/1‑93 (MQ=255)

                        aCAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTaa                                                                 <  2:4287427/93‑1 (MQ=255)

                          aGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACg                                                               <  2:6495210/93‑1 (MQ=255)

                           gctgctGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACgg                                                              <  2:488395/93‑1 (MQ=255)

                            ctgctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGt                                                             <  2:3996970/93‑1 (MQ=255)

                                tgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCg                                                         <  2:7624636/93‑1 (MQ=255)

                                   tGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACt                                                      <  2:580219/93‑1 (MQ=255)

                                    gCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTa                                                     <  2:1190116/93‑1 (MQ=255)

                                    gCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTa                                                     <  2:5063910/93‑1 (MQ=255)

                                     cAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAt                                                    <  2:3502982/93‑1 (MQ=255)

                                     cAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAt                                                    <  1:6906128/93‑1 (MQ=255)

                                     cAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAt                                                    >  1:7143529/1‑93 (MQ=255)

                                     cAACATCCGGCCAAACGTCTGGTGGTGGCCTGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAt                                                    <  1:772644/93‑1 (MQ=255)

                                      aaCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAtt                                                   >  2:1027250/1‑93 (MQ=255)

                                      aaCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAtt                                                   >  2:7362449/1‑93 (MQ=255)

                                      aaCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTAtt                                                   >  1:1506853/1‑93 (MQ=255)

                                       aCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTc                                                  >  1:3238327/1‑93 (MQ=255)

                                       aCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTc                                                  <  2:1698378/93‑1 (MQ=255)

                                       aCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTc                                                  >  2:4728334/1‑93 (MQ=255)

                                       aCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTc                                                  <  1:2912385/93‑1 (MQ=255)

                                         aTCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCcg                                                >  2:5132943/1‑93 (MQ=255)

                                          tCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCcgc                                               >  1:7332552/1‑93 (MQ=255)

                                              ggCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGaca                                           <  1:2694183/93‑1 (MQ=255)

                                                 aaaCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGATGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACaa                                        <  1:261045/93‑1 (MQ=255)

                                                  aaCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACaaa                                       >  1:3890639/1‑93 (MQ=255)

                                                  aaCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACaaa                                       >  2:4416772/1‑93 (MQ=255)

                                                  aaCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACaaa                                       >  2:6880050/1‑93 (MQ=255)

                                                   aCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAAt                                      <  1:1447949/93‑1 (MQ=255)

                                                   aCGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAAt                                      <  2:1940460/93‑1 (MQ=255)

                                                    cGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATc                                     <  1:4450108/93‑1 (MQ=255)

                                                    cGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATc                                     <  1:7906931/93‑1 (MQ=255)

                                                    cGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATc                                     <  1:664547/93‑1 (MQ=255)

                                                     gTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCg                                    <  1:2038499/93‑1 (MQ=255)

                                                      tCTGGTGTTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCgg                                   <  1:6355582/93‑1 (MQ=255)

                                                      tCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCgg                                   <  1:427643/93‑1 (MQ=255)

                                                       ttggtggTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGt                                  <  1:5367632/92‑1 (MQ=255)

                                                        tggtggTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTg                                 <  2:7496109/93‑1 (MQ=255)

                                                          gtggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTgcg                               >  2:6210605/1‑93 (MQ=255)

                                                          gtggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTgcg                               <  2:7209021/93‑1 (MQ=255)

                                                           tggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTgcgc                              <  1:572698/93‑1 (MQ=255)

                                                            ggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTgcgcg                             >  2:6365497/1‑92 (MQ=255)

                                                            ggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCt                             <  2:2384018/93‑1 (MQ=255)

                                                             gtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTg                            >  1:1203853/1‑93 (MQ=255)

                                                             gtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTg                            <  2:3037165/93‑1 (MQ=255)

                                                             gtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTg                            >  1:592069/1‑93 (MQ=255)

                                                              tgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTgg                           >  1:4114606/1‑93 (MQ=255)

                                                              tgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTgg                           >  1:7781769/1‑93 (MQ=255)

                                                              tgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTgg                           >  2:5261446/1‑93 (MQ=255)

                                                               gACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGc                          >  1:3795985/1‑93 (MQ=255)

                                                                aCCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCg                         >  1:5265208/1‑93 (MQ=255)

                                                                   tgATTTATCAGCCGCAACAACGCTGGTGAACCGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGgt                      <  1:4720657/92‑1 (MQ=255)

                                                                      tttATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGtt                   <  1:698978/93‑1 (MQ=255)

                                                                       ttATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGttt                  <  1:5757322/93‑1 (MQ=255)

                                                                       ttATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGttt                  <  1:5504924/93‑1 (MQ=255)

                                                                        tATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtc                 >  1:5170118/1‑93 (MQ=255)

                                                                         aTCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtct                <  2:707910/93‑1 (MQ=255)

                                                                         aTCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtct                <  1:4422892/93‑1 (MQ=255)

                                                                         aTCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtct                <  1:7269175/93‑1 (MQ=255)

                                                                          tCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtctc               <  1:6023324/93‑1 (MQ=255)

                                                                          tCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtctc               <  2:1399107/93‑1 (MQ=255)

                                                                             gCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTcgcg            >  1:5440124/1‑93 (MQ=255)

                                                                             gCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTcgcg            >  2:743693/1‑93 (MQ=255)

                                                                             gCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTcgcg            >  2:3912030/1‑93 (MQ=255)

                                                                               cGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGGGTTTCTCGCGTa          >  1:4317475/1‑93 (MQ=255)

                                                                                gCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTaa         <  1:3620537/93‑1 (MQ=255)

                                                                                 caacaaCGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAAc        <  1:6237619/93‑1 (MQ=255)

                                                                                 caacaaCGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGGGTTTCTCGCGTAAc        >  1:2132336/1‑93 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      <  1:1678170/93‑1 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      <  1:3175274/93‑1 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      <  1:3966454/93‑1 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTACGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      <  2:6753898/93‑1 (MQ=255)

                                                                                    caacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCAt     >  1:5051375/1‑93 (MQ=255)

                                                                                     aacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATc    <  1:2044537/93‑1 (MQ=255)

                                                                                      acGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCt   <  1:7835209/93‑1 (MQ=255)

AGGGCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCCGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTG  >  NC_000913/4127921‑4128100

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