BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,567,459	T→C	N192S (AAC→AGC)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,567,459	0	T	C	97.7%	449.9 / ‑1.3	133	N192S (AAC→AGC)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (54/75); total (56/77)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 133 total aligned reads displayed.

AAACTGCGGTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCAC  >  NC_000913/3567362‑3567556

tccgatctgTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTg                                                                                                 <  2:906076/92‑1 (MQ=255)

aaaCTGCGGTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTg                                                                                                 <  2:92287/100‑1 (MQ=255)

     gCGGTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTc                                                                                            <  2:2089372/100‑1 (MQ=255)

      cGGTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTcc                                                                                           <  2:510158/100‑1 (MQ=255)

        gTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAg                                                                                         <  1:470886/100‑1 (MQ=255)

         ttCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGc                                                                                        <  1:2782623/100‑1 (MQ=255)

          tCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGcc                                                                                       <  2:425329/100‑1 (MQ=255)

           cGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGccc                                                                                      >  1:1738303/1‑100 (MQ=255)

           cGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGccc                                                                                      <  2:2519597/100‑1 (MQ=255)

           cGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGccc                                                                                      >  2:1159303/1‑100 (MQ=255)

            ggTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGCTTTGTCCGCTGAATTCCAGcc                                                                                       >  2:1129003/1‑98 (MQ=255)

            ggTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCa                                                                                     <  1:2679795/100‑1 (MQ=255)

             gcgATCGCGCGAGCGTAGGTTGCACCGGCCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGTTGAATTCCAGc                                                                                        <  2:2872609/94‑1 (MQ=255)

             gTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCAt                                                                                    >  1:1811322/1‑100 (MQ=255)

             gTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCAt                                                                                    <  1:592577/100‑1 (MQ=255)

              tGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATg                                                                                   <  1:2824482/100‑1 (MQ=255)

              tGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATg                                                                                   <  2:2435991/100‑1 (MQ=255)

              tGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATg                                                                                   <  1:1431607/100‑1 (MQ=255)

                 tctcGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAat                                                                                >  2:67448/1‑100 (MQ=255)

                  ctcGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAata                                                                               <  1:67449/100‑1 (MQ=255)

                   tcGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAatat                                                                              >  2:1154252/1‑100 (MQ=255)

                    cgcgAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATAtt                                                                             <  1:536151/100‑1 (MQ=255)

                     gagAGCGTAGGTTAGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTa                                                                            >  1:2056157/3‑100 (MQ=255)

                      cgAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTaa                                                                           <  2:31553/100‑1 (MQ=255)

                      cgAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTaa                                                                           <  1:1154263/100‑1 (MQ=255)

                       gAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTaaa                                                                          <  1:2371404/100‑1 (MQ=255)

                       gAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTaaa                                                                          <  1:2774151/100‑1 (MQ=255)

                          cGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGaa                                                                       <  1:1932541/100‑1 (MQ=255)

                             aGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGa                                                                    >  1:2209335/1‑100 (MQ=255)

                             aGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGa                                                                    <  2:524111/100‑1 (MQ=255)

                              ggTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGAc                                                                   >  2:1450838/1‑100 (MQ=255)

                              ggTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGAc                                                                   <  1:2671478/100‑1 (MQ=255)

                                  ggACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATg                                                               <  1:1412116/100‑1 (MQ=255)

                                   gACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATgg                                                              <  2:1485829/100‑1 (MQ=255)

                                     cTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATgg                                                              <  2:1164074/98‑1 (MQ=255)

                                     cTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCa                                                            >  1:2082369/1‑100 (MQ=255)

                                       gACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAAt                                                          >  1:47633/1‑100 (MQ=255)

                                       gACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAAt                                                          >  1:2514474/1‑100 (MQ=255)

                                       gACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAAt                                                          >  1:1219535/1‑100 (MQ=255)

                                          cgccgTAATGTGATCGGCATAGTACAGACCGGCCTTCAAGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAACGAATGACCATGGCAATTg                                                        <  2:2971987/99‑1 (MQ=255)

                                           gccgTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGa                                                      >  1:289789/1‑100 (MQ=255)

                                           gccgTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGa                                                      >  2:2307178/1‑100 (MQ=255)

                                           gccgTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGa                                                      <  1:889556/100‑1 (MQ=255)

                                             cgTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATg                                                    <  2:1637148/100‑1 (MQ=255)

                                              gTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGt                                                   >  1:2891149/1‑100 (MQ=255)

                                                 aTGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcat                                                <  1:2179951/100‑1 (MQ=255)

                                                    tgATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattca                                             <  1:422888/100‑1 (MQ=255)

                                                    tgATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattca                                             <  1:2664306/100‑1 (MQ=255)

                                                    tgATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattca                                             >  1:2469906/1‑100 (MQ=255)

                                                     gATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattcat                                            >  1:2653077/1‑100 (MQ=255)

                                                     gATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattcat                                            >  2:26134/1‑100 (MQ=255)

                                                     gATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGtcattcat                                            >  1:1466018/1‑100 (MQ=255)

                                                      aTCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatg                                           <  1:1991765/100‑1 (MQ=255)

                                                      aTCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatg                                           <  1:1009574/100‑1 (MQ=255)

                                                       tCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgt                                          >  1:1055657/1‑100 (MQ=255)

                                                        cGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtg                                         >  1:113831/1‑100 (MQ=255)

                                                        cGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtg                                         <  2:2005066/100‑1 (MQ=255)

                                                         ggCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtga                                        >  1:771338/1‑100 (MQ=255)

                                                         ggCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtga                                        <  1:342746/100‑1 (MQ=255)

                                                         ggCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtga                                        >  1:1685449/1‑100 (MQ=255)

                                                         ggCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtga                                        <  2:1494734/100‑1 (MQ=255)

                                                         ggCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtga                                        >  1:1424614/1‑100 (MQ=255)

                                                             tAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtg                                         <  2:1879416/95‑1 (MQ=255)

                                                               gTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCTATTGGATGTCATTCAGTTGATGTGc                                   >  2:32482/1‑99 (MQ=255)

                                                               gTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCa                                  <  1:1559479/100‑1 (MQ=255)

                                                               gTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCa                                  >  2:207891/1‑100 (MQ=255)

                                                                tACAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCAt                                 >  1:2491771/1‑100 (MQ=255)

                                                                 aCAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATa                                <  2:745642/100‑1 (MQ=255)

                                                                 aCAGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATa                                >  2:2242877/1‑100 (MQ=255)

                                                                   aGACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATaaa                              <  2:527461/100‑1 (MQ=255)

                                                                    gACCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATaaaa                             >  2:2882616/1‑100 (MQ=255)

                                                                     aCCGGCCTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAAc                            <  1:885429/100‑1 (MQ=255)

                                                                         gCCTTCAGGAAAGTGATTTGTCCGGTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCACGTGATGTGCATAAAACATGc                        <  2:1287837/100‑1 (MQ=255)

                                                                          ccTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGcc                       <  1:511775/100‑1 (MQ=255)

                                                                           cTTCATGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCt                      >  2:1202776/1‑100 (MQ=255)

                                                                           cTTCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCt                      >  1:628677/1‑100 (MQ=255)

                                                                            ttCAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCtt                     <  1:212679/100‑1 (MQ=255)

                                                                              cAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGa                   <  1:2608368/100‑1 (MQ=255)

                                                                              cAGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGa                   >  2:812252/1‑100 (MQ=255)

                                                                               tccgatctgaTTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGAt                  <  2:1345046/92‑1 (MQ=255)

                                                                               aGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGAt                  <  2:631839/100‑1 (MQ=255)

                                                                               aGGAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGAt                  <  2:1738288/100‑1 (MQ=255)

                                                                                ggAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATa                 >  2:2920845/1‑100 (MQ=255)

                                                                                 gAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAg                <  2:2514452/100‑1 (MQ=255)

                                                                                 gAAAGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAg                >  2:287847/1‑100 (MQ=255)

                                                                                  aaaGAGATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAgg               <  2:1702510/100‑1 (MQ=255)

                                                                                  aaaGAGATTTGTCCGCTGAATTCCAGCAAATGAATATTACAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAgg               <  1:1202787/100‑1 (MQ=255)

                                                                                    agagATTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGcc             <  1:812401/100‑1 (MQ=255)

                                                                                     gagaTTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCa            <  2:2532213/100‑1 (MQ=255)

                                                                                       gaTTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAgatcgga          >  1:1345059/1‑94 (MQ=255)

                                                                                        aTTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGt         <  2:1962882/100‑1 (MQ=255)

                                                                                        aTTTGTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGt         >  1:960573/1‑100 (MQ=255)

                                                                                        aTTTGTCCGCTGAATTCCAGCCCAGGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGt         <  2:2655668/100‑1 (MQ=255)

                                                                                            gTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTtgtg     <  2:2556024/100‑1 (MQ=255)

                                                                                            gTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTtgtg     <  2:467520/100‑1 (MQ=255)

                                                                                            gTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTtgtg     <  2:2222041/100‑1 (MQ=255)

                                                                                            gTCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTtgtg     <  2:2555959/100‑1 (MQ=255)

                                                                                             tCCGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGc    <  2:1055647/100‑1 (MQ=255)

                                                                                              ccGCTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGca   >  2:1830855/1‑100 (MQ=255)

                                                                                               cGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGcac  >  2:37765/1‑100 (MQ=255)

AAACTGCGGTTCGGTGATCTCGCGAGCGTAGGTTGGACTGACCGCCGTAATGTGATCGGCATAGTACAGACCGGCCTTCAGGAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCAC  >  NC_000913/3567362‑3567556

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