BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,385,769	T→C	N342N (AAT→AAC)	tyrR →	TyrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,385,769	0	T	C	99.5%	705.9 / ‑4.7	218	N342N (AAT→AAC)	tyrR	TyrR
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (104/113); total (105/113)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 223 total aligned reads displayed.

GCGTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAATCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCCGCCGCTACGTGA  >  NC_000913/1385679‑1385858

gCGTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCt                                                                                         <  2:4466922/93‑1 (MQ=255)

gCGTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCt                                                                                         >  2:3372776/1‑93 (MQ=255)

 cGTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTg                                                                                        <  2:1776281/93‑1 (MQ=255)

  gTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTgg                                                                                       >  2:2828900/1‑93 (MQ=255)

   tttCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAATCTggc                                                                                      >  1:4762858/1‑92 (MQ=255)

    ttCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTc                                                                                     <  2:3745295/93‑1 (MQ=255)

       cTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGaa                                                                                  >  1:6046186/1‑93 (MQ=255)

        ttAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAAc                                                                                 >  1:2016142/1‑93 (MQ=255)

        ttAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAAc                                                                                 <  1:5004742/93‑1 (MQ=255)

        ttAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAAc                                                                                 >  2:1910794/1‑93 (MQ=255)

        ttAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAAc                                                                                 >  1:1333995/1‑93 (MQ=255)

         tAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACt                                                                                <  1:5026946/93‑1 (MQ=255)

         tAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACt                                                                                >  2:601443/1‑93 (MQ=255)

         tAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACt                                                                                <  2:537304/93‑1 (MQ=255)

         tAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACt                                                                                >  2:5309469/1‑93 (MQ=255)

         tAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACt                                                                                >  2:4847032/1‑93 (MQ=255)

          aaTGATGGCACTTTCCGTCGGGTTGGCGAAGACCCTGAGGTGCATGTCGCTGTGCGGGGGATTTGCGCTACGCAGAAGAACCTGGTCGAACTg                                                                               >  2:878497/1‑93 (MQ=255)

          aaTGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTg                                                                               <  2:5810910/93‑1 (MQ=255)

           atgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTgg                                                                              <  2:4874965/93‑1 (MQ=255)

           atgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTgg                                                                              <  2:5678716/93‑1 (MQ=255)

           atgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTgg                                                                              <  1:2480277/93‑1 (MQ=255)

           atgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTgg                                                                              <  1:466194/93‑1 (MQ=255)

           atgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTgg                                                                              <  2:5488509/93‑1 (MQ=255)

            tgatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGt                                                                             <  2:4501274/93‑1 (MQ=255)

             gatgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTg                                                                            >  2:3243534/1‑93 (MQ=255)

               tgGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCa                                                                          >  1:5741007/1‑93 (MQ=255)

                ggCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCaa                                                                         >  1:4072124/1‑93 (MQ=255)

                  cACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCaaaa                                                                       <  2:313299/93‑1 (MQ=255)

                   aCTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCaaaag                                                                      <  2:273125/93‑2 (MQ=255)

                    cTTTCCGTCGGGTTGGCTAAGACCATGAGGTGCATGTCGATGTGCGGGTAATTTGCGCTACGCAGAAGCACGTCGTCGAACTGGGGCAAAAAg                                                                     <  2:823147/93‑1 (MQ=255)

                    cTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAg                                                                     >  1:2215862/1‑93 (MQ=255)

                    cTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAg                                                                     >  2:2166882/1‑93 (MQ=255)

                     tttCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAgg                                                                    >  2:3460314/1‑93 (MQ=255)

                     tgtCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGACTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAgg                                                                    <  1:4430604/91‑1 (MQ=255)

                      ttCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGc                                                                   <  2:290087/93‑1 (MQ=255)

                       tCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCa                                                                  <  1:2361686/93‑1 (MQ=255)

                       tCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCa                                                                  <  1:3907687/93‑1 (MQ=255)

                       tCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCa                                                                  <  1:138531/93‑1 (MQ=255)

                       tCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCa                                                                  <  2:5068558/93‑1 (MQ=255)

                        tcGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCAt                                                                 <  2:5167417/92‑1 (MQ=255)

                        ccGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCAt                                                                 >  1:1207954/1‑93 (MQ=255)

                        ccGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCAt                                                                 <  1:3685084/93‑1 (MQ=255)

                         cGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATg                                                                >  1:2879902/1‑93 (MQ=255)

                          gTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGt                                                               >  1:5995812/1‑93 (MQ=255)

                           tCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGtt                                                              >  1:4690389/1‑93 (MQ=255)

                           tCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGtt                                                              >  1:3855357/1‑93 (MQ=255)

                           tCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGtt                                                              <  1:2367890/93‑1 (MQ=255)

                             gggTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTcc                                                            >  1:1680296/1‑93 (MQ=255)

                             gggTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTcc                                                            >  2:373858/1‑93 (MQ=255)

                             gggTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTcc                                                            <  2:2452494/93‑1 (MQ=255)

                             gggTTGACGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTcc                                                            <  1:864451/93‑1 (MQ=255)

                              ggTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCg                                                           <  2:5452967/93‑1 (MQ=255)

                               gTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGt                                                          <  2:5912322/93‑1 (MQ=255)

                                 tGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGa                                                        <  2:4746963/93‑1 (MQ=255)

                                   gCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAg                                                      <  2:2100803/93‑1 (MQ=255)

                                   gCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAg                                                      <  2:3654952/93‑1 (MQ=255)

                                    cGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGa                                                     <  2:4005687/93‑1 (MQ=255)

                                      aaGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtc                                                   >  1:4889616/1‑93 (MQ=255)

                                       aGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtct                                                  >  2:5601215/1‑93 (MQ=255)

                                       aGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtct                                                  >  2:3051177/1‑93 (MQ=255)

                                        gACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtctc                                                 <  1:110543/93‑1 (MQ=255)

                                        gACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtctc                                                 >  2:4411809/1‑93 (MQ=255)

                                        gACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGAtctc                                                 >  1:3211557/1‑93 (MQ=255)

                                          ccATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCta                                               <  1:3994549/93‑1 (MQ=255)

                                          ccATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCta                                               <  1:3953647/93‑1 (MQ=255)

                                           cATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCtat                                              <  2:1160482/93‑1 (MQ=255)

                                           cATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCtat                                              >  2:163029/1‑93 (MQ=255)

                                              gAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCtattat                                           >  1:4202209/1‑93 (MQ=255)

                                               aGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATc                                          >  2:172698/1‑93 (MQ=255)

                                               aGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATc                                          <  2:5236872/93‑1 (MQ=255)

                                                ggTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCg                                         >  1:5985974/1‑93 (MQ=255)

                                                 gTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGt                                        <  1:1815696/93‑1 (MQ=255)

                                                  tGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTc                                       <  1:5087508/93‑1 (MQ=255)

                                                  tGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTc                                       >  1:5814763/1‑93 (MQ=255)

                                                  tGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTc                                       <  2:12748/93‑1 (MQ=255)

                                                    cATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTg                                     <  1:4750976/93‑1 (MQ=255)

                                                     aTGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGa                                    <  1:4349471/93‑1 (MQ=255)

                                                      tGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGaa                                   <  1:6153744/93‑1 (MQ=255)

                                                       gTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAAc                                  >  2:4062980/1‑93 (MQ=255)

                                                       gTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAAc                                  <  1:3139056/93‑1 (MQ=255)

                                                       gTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCCTCTGAAc                                  >  1:1669802/1‑93 (MQ=255)

                                                         cGATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACgt                                <  1:3458044/93‑1 (MQ=255)

                                                          gATGTGCGTGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACgtg                               <  1:5625796/93‑1 (MQ=255)

                                                          gATGTGCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACgtg                               >  1:762632/1‑93 (MQ=255)

                                                            tgtgCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGtt                             >  2:4952173/1‑93 (MQ=255)

                                                              tgCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGACCGTGTTGa                           <  1:515221/93‑1 (MQ=255)

                                                               gCGGGTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGAc                          <  2:6284470/93‑1 (MQ=255)

                                                                 gggTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGc                        >  2:77068/1‑93 (MQ=255)

                                                                   gTGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTc                      >  1:4587496/1‑93 (MQ=255)

                                                                    tGATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCa                     <  1:3028365/93‑1 (MQ=255)

                                                                     gATTTGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGGGTTGACGCTCaa                    <  1:2393129/93‑1 (MQ=255)

                                                                       tttGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGCACGTGTTGACGCTCAATc                  >  1:5429660/1‑93 (MQ=255)

                                                                       tttGCGCTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATc                  >  1:1250480/1‑93 (MQ=255)

                                                                          gcgcTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTgc               <  1:5199310/93‑1 (MQ=255)

                                                                             cTACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTgccgc            <  1:4920923/93‑1 (MQ=255)

                                                                              tACGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTgccgcc           >  2:2236417/1‑93 (MQ=255)

                                                                                cGCAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCcgccgc         <  2:4825982/93‑1 (MQ=255)

                                                                                  cAGAAGAACCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCCGCCGCTa       >  2:5511524/1‑93 (MQ=255)

                                                                                     aagaaCCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCCGCCGCTACGt    >  1:73202/1‑93 (MQ=255)

                                                                                     aagaaCCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCCGCCGCTACGt    <  2:5532511/93‑1 (MQ=255)

GCGTTTCCTTAATGATGGCACTTTCCGTCGGGTTGGCGAAGACCATGAGGTGCATGTCGATGTGCGGGTGATTTGCGCTACGCAGAAGAATCTGGTCGAACTGGTGCAAAAAGGCATGTTCCGTGAAGATCTCTATTATCGTCTGAACGTGTTGACGCTCAATCTGCCGCCGCTACGTGA  >  NC_000913/1385679‑1385858

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