BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999	A→G	T555A (ACT→GCT)	thrA →	aspartokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999	0	A	G	99.7%	956.2 / NA	294	T555A (ACT→GCT)	thrA	aspartokinase I
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (167/126); minor base T (1/0); total (168/126)
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 303 total aligned reads displayed.

CTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTC  >  NC_000913/1909‑2087

cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCt                                                                                        >  1:1399947/1‑93 (MQ=255)

cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCt                                                                                        >  2:2926167/1‑93 (MQ=255)

cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCt                                                                                        <  1:2151411/93‑1 (MQ=255)

cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCt                                                                                        >  1:912946/1‑93 (MQ=255)

  ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTc                                                                                      >  2:3961025/1‑93 (MQ=255)

   gctcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCCTcc                                                                                     <  1:2514656/90‑1 (MQ=255)

   gcgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTcc                                                                                     <  2:6682139/93‑1 (MQ=255)

     gcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccag                                                                                   <  1:3265674/93‑1 (MQ=255)

      cAAGCCAGAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagc                                                                                  >  1:5402085/1‑93 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                 >  2:3811141/1‑93 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                 >  1:351321/1‑93 (MQ=255)

        aGCCAAAGATCCGTTTAATCTCTGTCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcca                                                                                <  1:3387225/93‑1 (MQ=255)

        aGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcca                                                                                >  1:1460719/1‑93 (MQ=255)

          ccAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAgg                                                                              >  2:5096876/1‑93 (MQ=255)

            aaaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCa                                                                            >  2:3328251/1‑93 (MQ=255)

            aaaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCa                                                                            >  1:1031736/1‑93 (MQ=255)

             aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAg                                                                           >  1:4446837/1‑93 (MQ=255)

             aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAg                                                                           >  2:2505406/1‑93 (MQ=255)

             aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAg                                                                           >  1:1831987/1‑93 (MQ=255)

             aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAg                                                                           >  1:2687745/1‑93 (MQ=255)

              agagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGt                                                                          <  2:6161739/93‑1 (MQ=255)

               gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTg                                                                         >  1:4861481/1‑93 (MQ=255)

               gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTg                                                                         >  1:4158504/1‑93 (MQ=255)

               gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTg                                                                         <  2:2791181/93‑1 (MQ=255)

                agCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTgg                                                                        >  1:3564047/1‑93 (MQ=255)

                 gCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGc                                                                       <  1:3019935/93‑1 (MQ=255)

                 gCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGc                                                                       >  2:4266986/1‑93 (MQ=255)

                   cGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCgg                                                                     >  1:2814808/1‑93 (MQ=255)

                    gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGa                                                                    <  1:3718318/93‑1 (MQ=255)

                    gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGa                                                                    >  2:2204961/1‑93 (MQ=255)

                     tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGAt                                                                   <  1:5506421/93‑1 (MQ=255)

                      ttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATc                                                                  <  2:2421134/93‑1 (MQ=255)

                      ttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATc                                                                  >  1:6317994/1‑93 (MQ=255)

                       tAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCa                                                                 >  2:2391719/1‑93 (MQ=255)

                         aTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAat                                                               <  1:3359478/93‑1 (MQ=255)

                           ctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAatat                                                             >  1:4895668/1‑93 (MQ=255)

                                   cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACtt                                                     >  1:6733857/1‑93 (MQ=255)

                                    ttAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTc                                                    <  1:2925085/93‑1 (MQ=255)

                                      aaTTCGCCTGGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCt                                                  <  2:1364146/93‑1 (MQ=255)

                                      aaTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCt                                                  <  1:5758626/93‑1 (MQ=255)

                                      aaTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCt                                                  <  1:6498521/93‑1 (MQ=255)

                                       aTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTg                                                 <  2:191839/93‑1 (MQ=255)

                                         tCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcg                                               <  1:6269896/93‑1 (MQ=255)

                                          cGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcgc                                              <  1:2313208/93‑1 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                             >  2:6291620/1‑93 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                             >  2:1070534/1‑93 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                             >  1:3822653/1‑93 (MQ=255)

                                             cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                           >  1:1565190/1‑93 (MQ=255)

                                             cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                           <  2:5318782/93‑1 (MQ=255)

                                              tCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAg                                          <  2:5007345/93‑1 (MQ=255)

                                               cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                         >  2:6228877/1‑93 (MQ=255)

                                               cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                         >  1:5009170/1‑93 (MQ=255)

                                               cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCATGCGCGAAgg                                         >  2:5517020/1‑93 (MQ=255)

                                                gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGt                                        <  2:1659042/93‑1 (MQ=255)

                                                gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGt                                        >  2:4592183/1‑93 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                       <  1:6074415/93‑1 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                       <  2:1179855/93‑1 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                       <  2:3021637/93‑1 (MQ=255)

                                                  gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGttt                                      <  1:3990449/93‑1 (MQ=255)

                                                   aaaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGGTTc                                     >  1:3914862/1‑93 (MQ=255)

                                                    aaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTcc                                    >  1:4104592/1‑93 (MQ=255)

                                                    aaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTcc                                    <  2:266688/93‑1 (MQ=255)

                                                     aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCa                                   >  1:4214974/1‑93 (MQ=255)

                                                      gAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTTCAc                                  >  1:5853337/1‑93 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACg                                 <  2:3785398/93‑1 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACg                                 <  1:1094160/93‑1 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACg                                 <  1:3153173/93‑1 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGGTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGGCGTCCTTCGCGAAGGTTTCCACg                                 >  1:4653372/1‑93 (MQ=255)

                                                        atatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGt                                <  1:1984698/93‑1 (MQ=255)

                                                          atcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTg                              <  2:767683/93‑1 (MQ=255)

                                                           tcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGt                             <  1:5227677/93‑1 (MQ=255)

                                                              tcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCAc                          >  1:240933/1‑93 (MQ=255)

                                                              actgttgAGCCCGGTCATTGTTGCCTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCTCGTTGTCAc                          <  1:1291099/92‑1 (MQ=255)

                                                                tgctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGc                        >  1:1615685/1‑93 (MQ=255)

                                                                tgctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGc                        >  2:5155063/1‑93 (MQ=255)

                                                                  ctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCg                      <  2:3576265/93‑1 (MQ=255)

                                                                    gAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGaa                    >  2:2257811/1‑93 (MQ=255)

                                                                    gAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGaa                    <  2:1098656/93‑1 (MQ=255)

                                                                      aCCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACa                  >  2:2997785/1‑93 (MQ=255)

                                                                       cccGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaa                 <  2:6568005/93‑1 (MQ=255)

                                                                          ggTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa              >  2:3146588/1‑93 (MQ=255)

                                                                          ggTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa              <  1:6598708/93‑1 (MQ=255)

                                                                           gTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAg             <  2:1535948/93‑1 (MQ=255)

                                                                            tCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAgg            >  1:2192963/1‑93 (MQ=255)

                                                                              aTTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGCACAAAAAGGcc          >  2:2129703/1‑93 (MQ=255)

                                                                              aTTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc          >  2:3514524/1‑93 (MQ=255)

                                                                              aTTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc          >  2:608880/1‑93 (MQ=255)

                                                                              aTTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc          <  1:363835/93‑1 (MQ=255)

                                                                                tgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAGGGCCaa        <  1:5556039/93‑1 (MQ=255)

                                                                                tgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCaa        <  2:4438130/93‑1 (MQ=255)

                                                                                 gttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAac       <  2:5581594/93‑1 (MQ=255)

                                                                                 gttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAac       <  2:4070819/93‑1 (MQ=255)

                                                                                 gttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGCTTTCCACGTTGTCACGCCGAACAAAAAGGCCAac       >  2:3385134/1‑93 (MQ=255)

                                                                                    gACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACAcc    <  2:6317994/93‑1 (MQ=255)

                                                                                     aCTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCt   >  1:4332074/1‑93 (MQ=255)

                                                                                     aCTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCt   <  1:5113945/93‑1 (MQ=255)

                                                                                     aCTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTACCACGTTGTCACGCCGAACAAAAAGGCCAACACCt   <  2:180816/93‑1 (MQ=255)

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc  >  2:5662548/1‑93 (MQ=255)

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc  <  2:5831462/93‑1 (MQ=255)

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc  <  2:4520434/93‑1 (MQ=255)

CTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTC  >  NC_000913/1909‑2087

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