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%	552.0 / ‑5.7	156	A199A (GCC→GCT)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (75/80); total (76/80)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 165 total aligned reads displayed.

AGGGCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCCGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTG  >  NC_000913/4127921‑4128100

aGGGCTTTCTTCCCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTgg                                                                                         <  2:3973813/93‑1 (MQ=255)

aGGGCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCTGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTgg                                                                                         >  1:3520763/1‑93 (MQ=255)

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

  ggCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTg                                                                                       >  1:4516696/1‑93 (MQ=255)

  ggCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTg                                                                                       <  2:2578972/93‑1 (MQ=255)

  ggCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTg                                                                                       <  2:2987262/93‑1 (MQ=255)

    cTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaa                                                                                     >  1:3456655/1‑93 (MQ=255)

    cTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGaa                                                                                     <  1:1914044/93‑1 (MQ=255)

      ttcttcCCCGTTTCTGCAACAGGTGCTGGTGGAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAAc                                                                                   <  1:2355139/93‑1 (MQ=255)

      ttCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAAc                                                                                   >  2:1587670/1‑93 (MQ=255)

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

         ttACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGt                                                                                <  2:3482981/93‑1 (MQ=255)

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

          tACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGTAACGGtc                                                                               >  2:2560088/1‑92 (MQ=255)

          tACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtg                                                                               >  1:3562786/1‑93 (MQ=255)

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

            cccGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgct                                                                             <  2:1797204/93‑1 (MQ=255)

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

              cGTTGCTGCAACAGCTGCTGGTGCAACCTCCGGGCAAACGTCTGGTGGTGACCCGATTTATCAGCCGCAACAACGCTGGTGAAACGGtgctgc                                                                           <  1:1322822/93‑1 (MQ=255)

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

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

               gTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTgctgct                                                                          <  1:243166/93‑1 (MQ=255)

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

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

                ttGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGctgctg                                                                         >  1:2352179/1‑93 (MQ=255)

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

                 tgctgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCGGCTgg                                                                        >  1:875995/1‑93 (MQ=255)

                    tgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGTATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGc                                                                     <  2:2170393/93‑1 (MQ=255)

                    tgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGc                                                                     >  1:1820818/1‑93 (MQ=255)

                    tgcAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGc                                                                     <  1:1340877/93‑1 (MQ=255)

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

                      caccAGCTGCTGGTGCAACATCCGGGCAAGCGTCTGGTGGTGACCGGATTTATCCGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGt                                                                   <  2:1335278/90‑1 (MQ=255)

                       aaCAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTa                                                                  >  1:1119335/1‑93 (MQ=255)

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

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

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

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

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

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

                             tgctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGtt                                                            <  2:1816091/93‑1 (MQ=255)

                             tgctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGtt                                                            >  1:610701/1‑93 (MQ=255)

                              gctgGTGCAGCATCCGGGCAAACGTCTGGTGATGACCGGATTTATCAGCCGCACCAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTc                                                           <  2:2531058/93‑1 (MQ=255)

                              gctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGCAACGGTGCTGCTGGGGCGTAACGGTTc                                                           >  1:132471/1‑93 (MQ=255)

                              gctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTc                                                           <  1:594621/93‑1 (MQ=255)

                              gctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTc                                                           <  2:1565707/93‑1 (MQ=255)

                               ctgGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTcc                                                          <  2:1872844/93‑1 (MQ=255)

                                 ggTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTccta                                                        >  1:3077227/1‑91 (MQ=255)

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

                                   tGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACt                                                      >  1:873480/1‑93 (MQ=255)

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

                                       aCATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACACCGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTc                                                  <  2:3502868/93‑1 (MQ=255)

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

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

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

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

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

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

                                          tCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCCGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCcgc                                               >  1:1563067/1‑93 (MQ=255)

                                              gccAAATGTCGGGTGTTTACCTGAGTTATCAGCAGGAGCAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACAATTCCGCGtca                                           <  2:570075/91‑1 (MQ=255)

                                               gCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGacac                                          <  2:4750186/93‑1 (MQ=255)

                                                cAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGAcaca                                         <  2:4054139/93‑1 (MQ=255)

                                                   acttCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAAt                                      <  2:2679408/90‑1 (MQ=255)

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

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

                                                      tCTGGTGGTGCCCGGATTTATCAGCCGCAACAACTCTGGTGAAACGGTCCTGCTGGGGCGTAACGTTTCCGACTATTCCGCGACACAAATCgg                                   <  1:431774/93‑1 (MQ=255)

                                                      tCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCgg                                   >  2:109573/1‑93 (MQ=255)

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

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

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

                                                        tggtggGGACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTg                                 <  2:2797122/93‑1 (MQ=255)

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

                                                           tggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTgcgc                              >  2:3073064/1‑93 (MQ=255)

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

                                                            ggtgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCt                             >  1:46364/1‑93 (MQ=255)

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

                                                              tgACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTgg                           <  1:581218/93‑1 (MQ=255)

                                                               gACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGc                          <  2:1944754/93‑1 (MQ=255)

                                                               gACCGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGc                          <  2:1944755/93‑1 (MQ=255)

                                                                  cGGATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCggg                       <  1:1285573/93‑1 (MQ=255)

                                                                   ggATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCgggg                      >  1:480479/1‑92 (MQ=255)

                                                                   ggATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCgggg                      >  1:4245369/1‑92 (MQ=255)

                                                                   ggATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGgt                      >  1:4245368/1‑93 (MQ=255)

                                                                    gATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGgtg                     >  1:3186917/1‑93 (MQ=255)

                                                                    gATTTATCAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGgtg                     >  2:282790/1‑93 (MQ=255)

                                                                      tttATCAGCCGCAACAACGCTGGTGAAGCGGGTCTGCGGGGGCCTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGtt                   <  1:4333289/93‑1 (MQ=255)

                                                                        tATCAGCAGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTtc                 <  2:2864077/93‑1 (MQ=255)

                                                                           cAGCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTcg              <  2:32784/93‑1 (MQ=255)

                                                                             gCCGCAACAACGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTtgcg            >  1:920533/1‑93 (MQ=255)

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

                                                                                gCAACAACGCTGGTGAAACGGTGTTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGGTTCTCGCGTaa         >  1:4348077/1‑93 (MQ=255)

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

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

                                                                                  aacaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAAcc       <  1:2919138/93‑1 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      >  1:2580609/1‑93 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCa      >  2:158639/1‑93 (MQ=255)

                                                                                   acaacGCTGGTGAAACGGTACTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAAACa      >  2:3014403/1‑93 (MQ=255)

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

                                                                                     aacGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATc    >  1:2282636/1‑93 (MQ=255)

                                                                                      acGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCt   >  1:703483/1‑93 (MQ=255)

                                                                                      acGCTGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCt   >  2:581605/1‑93 (MQ=255)

AGGGCTTTCTTACCCGTTGCTGCAACAGCTGCTGGTGCAACATCCGGGCAAACGTCTGGTGGTGACCGGATTTATCAGCCGCAACAACGCCGGTGAAACGGTGCTGCTGGGGCGTAACGGTTCCGACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTG  >  NC_000913/4127921‑4128100

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