BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,530	T→C	P312P (CCT→CCC)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,530	0	T	C	99.5%	670.4 / ‑4.9	209	P312P (CCT→CCC)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (105/102); total (106/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 216 total aligned reads displayed.

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

taAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCgg                                                                                         <  1:4361529/93‑1 (MQ=255)

taAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCgg                                                                                         >  1:4130670/1‑93 (MQ=255)

  aGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTc                                                                                       <  2:3337586/93‑1 (MQ=255)

    aTCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTg                                                                                     >  2:417309/1‑93 (MQ=255)

     tCGATAAGTCATTTGTGCAAATGGGGAGTGTCGACGAAATCTCCGGTCATATTGTTGACAATATTGTCGAACTAGCGCGTAAGCCCGGCCTgt                                                                                    >  1:1445357/1‑92 (MQ=255)

      cGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAg                                                                                   <  1:3403484/93‑1 (MQ=255)

       gATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGGCCGGTCTGAga                                                                                  >  2:1507118/1‑92 (MQ=255)

       gATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGt                                                                                  <  2:1859951/93‑1 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTATCGCGTAAGCCCGGTCTGAGTa                                                                                 >  1:2090876/1‑93 (MQ=255)

        aTAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTa                                                                                 <  2:2594053/93‑1 (MQ=255)

          aaGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               >  1:3820752/1‑93 (MQ=255)

          aaGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               >  1:3639697/1‑93 (MQ=255)

          aaGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               >  2:4444780/1‑93 (MQ=255)

          aaGTCATTTAGGCAAAGGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATc                                                                               <  1:758991/93‑1 (MQ=255)

           aGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCg                                                                              >  2:3066963/1‑93 (MQ=255)

           aGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCg                                                                              <  1:1311605/93‑1 (MQ=255)

              cATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTgg                                                                           >  2:1138953/1‑93 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          <  2:4630138/93‑1 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          <  1:2033620/93‑1 (MQ=255)

               aTTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGc                                                                          >  2:4722062/1‑93 (MQ=255)

                tttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCg                                                                         >  2:3715068/1‑93 (MQ=255)

                tttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCg                                                                         >  2:1311028/1‑93 (MQ=255)

                 ttGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCgg                                                                        >  2:3034749/1‑93 (MQ=255)

                  tgtgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGa                                                                       >  2:1069848/1‑93 (MQ=255)

                    tgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAg                                                                     <  2:645594/93‑1 (MQ=255)

                    tgCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAg                                                                     <  1:3537829/93‑1 (MQ=255)

                     gCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgg                                                                    <  2:2690806/93‑1 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAgggga                                                                 >  2:2168068/1‑92 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGt                                                                 <  2:2535106/93‑1 (MQ=255)

                        aaTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGt                                                                 >  1:4553212/1‑93 (MQ=255)

                         aTGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTa                                                                <  1:4161466/93‑1 (MQ=255)

                          tGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGCAGGGGTAg                                                               >  2:4262429/1‑93 (MQ=255)

                          tGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAg                                                               <  1:3336706/93‑1 (MQ=255)

                          tGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAg                                                               >  2:717603/1‑93 (MQ=255)

                          tGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAg                                                               <  2:739369/93‑1 (MQ=255)

                           ggCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGa                                                              <  1:1691365/93‑1 (MQ=255)

                            gCGAGTGTCGACGAAATCTCCGGTCATATTGGGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGaa                                                             <  1:1036835/93‑1 (MQ=255)

                              gAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAc                                                           >  1:2544928/1‑93 (MQ=255)

                                gtgtCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAAccc                                                         <  2:4313609/93‑1 (MQ=255)

                                    cgacgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGACGGGGTAGAAACACAGGa                                                     >  1:1108307/1‑93 (MQ=255)

                                     gacgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAg                                                    >  2:2709248/1‑93 (MQ=255)

                                       cgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCa                                                  >  2:4017028/1‑93 (MQ=255)

                                       cgaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCa                                                  <  2:550220/93‑1 (MQ=255)

                                        gaAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAg                                                 >  1:61314/1‑93 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAgg                                                <  1:429853/93‑1 (MQ=255)

                                         aaaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAAGGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAgg                                                >  2:1545058/1‑93 (MQ=255)

                                          aaTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGAAACCCAGGAGCAGGc                                               >  2:1803722/1‑93 (MQ=255)

                                           aTCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCg                                              <  1:2439021/93‑1 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             >  2:1453478/1‑93 (MQ=255)

                                            tctcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCgg                                             <  1:1351436/93‑1 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           <  2:3187879/93‑1 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           >  1:2799527/1‑93 (MQ=255)

                                              tcCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAt                                           <  2:228273/93‑1 (MQ=255)

                                               ccGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAtt                                          <  2:1211026/93‑1 (MQ=255)

                                                cGGTCATATTGTGGACAATATTGTCGCACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAtgt                                         >  1:211181/1‑91 (MQ=255)

                                                cGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGAttt                                         >  1:920407/1‑93 (MQ=255)

                                                 ggTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTa                                        <  1:231277/93‑1 (MQ=255)

                                                 ggTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTa                                        <  1:3167050/93‑1 (MQ=255)

                                                  gTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTaa                                       <  2:1385517/93‑1 (MQ=255)

                                                   tCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAAt                                      >  2:4025807/1‑93 (MQ=255)

                                                    cATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATg                                     >  1:312171/1‑93 (MQ=255)

                                                    cATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATg                                     >  1:4240812/1‑93 (MQ=255)

                                                     atatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGa                                    <  2:3868487/93‑1 (MQ=255)

                                                     atatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGa                                    <  2:3626836/93‑1 (MQ=255)

                                                      tatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGAt                                   >  1:4712171/1‑93 (MQ=255)

                                                      tatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGAt                                   >  2:140929/1‑93 (MQ=255)

                                                      tatTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGAt                                   <  2:4265220/93‑1 (MQ=255)

                                                       atTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATc                                  <  2:2984306/93‑1 (MQ=255)

                                                       atTGTGGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATc                                  <  2:706154/93‑1 (MQ=255)

                                                         tgtgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCgg                                <  1:3455590/93‑1 (MQ=255)

                                                          gtgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGt                               <  2:3497214/93‑1 (MQ=255)

                                                           tgGACAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTa                              <  1:2992672/93‑1 (MQ=255)

                                                              aCAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAg                           <  1:3770859/93‑1 (MQ=255)

                                                              aCAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAg                           <  1:1451302/93‑1 (MQ=255)

                                                               cAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAgg                          >  1:3982618/1‑93 (MQ=255)

                                                               cAATATTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAgg                          >  1:3738997/1‑93 (MQ=255)

                                                                 atatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCg                        >  1:2342103/1‑93 (MQ=255)

                                                                  tatTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGt                       >  1:1594994/1‑93 (MQ=255)

                                                                   atTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGtt                      <  1:2966291/93‑1 (MQ=255)

                                                                   atTGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGtt                      <  2:1534901/93‑1 (MQ=255)

                                                                     tGTCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGAAAACCAGGAGCAGGCGGATTTAATGATCGGTAAAGGGGttta                    >  1:4724903/1‑91 (MQ=255)

                                                                       tCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACt                  >  1:2646594/1‑93 (MQ=255)

                                                                       tCGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACt                  <  2:3748057/93‑1 (MQ=255)

                                                                        cGAACTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtt                 <  1:1892746/93‑1 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtttt               <  1:2276502/93‑1 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtttt               <  1:4314444/93‑1 (MQ=255)

                                                                          aaCTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACtttt               >  1:1712967/1‑93 (MQ=255)

                                                                            cTAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGGAGAAACCCAGGAGCAGGGGGATTTCATGATCGGTAAGGGAGTTCACGTTTTg             >  2:4178170/1‑93 (MQ=255)

                                                                            cTAGCGCGTAAGACCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTg             >  2:2126433/1‑93 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGc            >  1:1212668/1‑93 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGc            >  2:169442/1‑93 (MQ=255)

                                                                             tAGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGCAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGAGTTCACTTTTTGc            >  1:3565721/1‑93 (MQ=255)

                                                                              aGCGCGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCa           >  2:1044936/1‑93 (MQ=255)

                                                                               gcgcGTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAg          <  2:2085611/93‑1 (MQ=255)

                                                                                  cgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGGTTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc       >  1:874444/1‑93 (MQ=255)

                                                                                  cgTAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc       <  2:2792924/93‑1 (MQ=255)

                                                                                    tAAGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTa     >  1:1985068/1‑93 (MQ=255)

                                                                                      aGCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAtt   <  1:3319277/93‑1 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  >  1:1738112/1‑93 (MQ=255)

                                                                                       gCCCGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt  >  1:671683/1‑93 (MQ=255)

TAAGATCGATAAGTCATTTGTGCAAATGGCGAGTGTCGACGAAATCTCCGGTCATATTGTGGACAATATTGTCGAACTAGCGCGTAAGCCTGGTCTGAGTATCGTGGCGGAAGGGGTAGAAACCCAGGAGCAGGCGGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTT  >  NC_000913/332440‑332619

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