BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,333,152	A→G	L16P (CTA→CCA)	yjdB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,333,152	0	A	G	99.1%	703.6 / ‑5.3	226	L16P (CTA→CCA)	yjdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (122/102); total (124/102)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

acCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggt                                                                                         >  1:951173/1‑93 (MQ=255)

 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  1:1936790/93‑1 (MQ=255)

     tttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgcg                                                                                    >  2:3001944/1‑93 (MQ=255)

     tttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAAATCAAAGAGGGTCTTTTTAGTAGgcg                                                                                    >  1:2473495/1‑93 (MQ=255)

        aaaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGGCTTTTTGGTAGGCGCtt                                                                                 >  2:2821229/1‑93 (MQ=255)

          aaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCa                                                                               >  1:1042725/1‑93 (MQ=255)

          aaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCa                                                                               >  1:921059/1‑93 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  2:4233488/93‑1 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  2:3410138/93‑1 (MQ=255)

           aaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCaa                                                                              <  1:3226443/93‑1 (MQ=255)

            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             <  1:2376461/93‑1 (MQ=255)

             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATCGCCAGGCGAGCAAATTCAAAGAGGGACTTTTAGGTAGGCGCCTAAACa                                                                            >  2:3359849/1‑93 (MQ=255)

             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            >  1:3839175/1‑93 (MQ=255)

             aGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACa                                                                            >  1:3777561/1‑93 (MQ=255)

              ggCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACAt                                                                           >  2:1273830/1‑93 (MQ=255)

                cAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCa                                                                         <  1:319914/93‑1 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        <  1:2917620/93‑1 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        <  1:2839947/93‑1 (MQ=255)

                  atatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACg                                                                       >  2:2250608/1‑93 (MQ=255)

                  atatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGATCTTCAACATCACg                                                                       >  1:1623662/1‑93 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      <  2:4338488/93‑1 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      <  2:2777905/93‑1 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    <  1:177716/93‑1 (MQ=255)

                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  1:3539797/93‑1 (MQ=255)

                        aGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGttt                                                                 >  2:3806653/1‑93 (MQ=255)

                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTca                                                               >  1:1261414/1‑92 (MQ=255)

                           cAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCa                                                              <  1:3629453/93‑1 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  2:3201967/93‑1 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  2:182346/93‑1 (MQ=255)

                              aTAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCg                                                           >  1:3377703/1‑93 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  2:4031892/1‑93 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  2:654296/1‑93 (MQ=255)

                                  agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACa                                                       <  2:1193764/93‑1 (MQ=255)

                                   gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaa                                                      >  1:2030029/1‑93 (MQ=255)

                                   gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACaa                                                      >  1:1954311/1‑93 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    >  2:811459/1‑93 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGCGTTTCCATCGAACAAAg                                                    >  2:2723251/1‑93 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGAAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGCTGTTTCCATCGCACAAAg                                                    >  2:3783261/1‑93 (MQ=255)

                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   >  1:2414705/1‑93 (MQ=255)

                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   <  2:4647714/93‑1 (MQ=255)

                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   >  2:501672/1‑93 (MQ=255)

                                      agaaaaaGCGTCGAACAATAGCCAGGCGAGTAAATGCAAAGAGGGTCTTTTTGGTAGGCGCTTCACCATCACGGTGTTTCCATCGAACAAAGt                                                   <  2:2509777/91‑1 (MQ=255)

                                         aaaaGGGGCCAACAATAGCCAGGTGAGGGAGTTCAAAGAGGGTCTTTTTGGTAGGCCATTCAACATCACGGCGTTTCCATCGAACAAAGTgcg                                                >  2:2149697/1‑93 (MQ=255)

                                         aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcg                                                >  2:1014418/1‑93 (MQ=255)

                                         aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcg                                                >  2:438112/1‑93 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              <  1:1701233/93‑1 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  1:4293926/1‑93 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              <  2:3873716/93‑1 (MQ=255)

                                             gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            <  2:512646/93‑1 (MQ=255)

                                             gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            <  2:3485830/93‑1 (MQ=255)

                                             gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            >  1:143348/1‑93 (MQ=255)

                                              cGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCatat                                           >  1:1842240/1‑93 (MQ=255)

                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCATTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          <  1:2127901/93‑1 (MQ=255)

                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        <  1:488673/93‑1 (MQ=255)

                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  1:1003568/93‑1 (MQ=255)

                                                   aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCg                                      >  1:2759425/1‑93 (MQ=255)

                                                    acaaTAGCCAGGCGAGTAAATTCAAAGAGGGGCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGGATATGCTCgc                                     >  2:2507862/1‑93 (MQ=255)

                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgca                                    >  2:753717/1‑93 (MQ=255)

                                                      aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaa                                   >  2:153696/1‑93 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  >  1:3759511/1‑93 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  >  2:1997335/1‑93 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  2:1615266/93‑1 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 >  1:645960/1‑93 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 >  2:1899047/1‑93 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 >  2:1586876/1‑93 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTAGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  2:2023956/93‑1 (MQ=255)

                                                          gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagca                               >  1:2888847/1‑93 (MQ=255)

                                                            cAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagcaaa                             >  1:2093592/1‑93 (MQ=255)

                                                            cAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagcaaa                             >  1:3107486/1‑93 (MQ=255)

                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            >  2:1694737/1‑93 (MQ=255)

                                                             aGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCaaaa                            <  1:82153/93‑1 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           <  2:3745495/93‑1 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  1:4536271/1‑93 (MQ=255)

                                                                cGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTa                         <  2:2827054/93‑1 (MQ=255)

                                                                cGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTa                         <  2:3952197/93‑1 (MQ=255)

                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTaa                        >  2:3624749/1‑93 (MQ=255)

                                                                  aGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAg                       >  1:4157572/1‑93 (MQ=255)

                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      <  2:2852112/93‑1 (MQ=255)

                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGCTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      <  1:3044249/93‑1 (MQ=255)

                                                                    tAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGcc                     >  1:4266737/1‑93 (MQ=255)

                                                                      aaTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCac                   >  2:2604451/1‑92 (MQ=255)

                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  <  1:2908129/93‑1 (MQ=255)

                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  >  1:4557694/1‑93 (MQ=255)

                                                                       aTTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAc                  <  2:5927/93‑1 (MQ=255)

                                                                        ttCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAcc                 >  1:1663845/1‑93 (MQ=255)

                                                                           aaaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTa              >  2:2102946/1‑93 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             <  2:3194416/93‑1 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             >  2:1475008/1‑93 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  2:3310397/93‑1 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  1:2700190/93‑1 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  2:695347/93‑1 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          <  2:3263186/93‑1 (MQ=255)

                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         <  1:4004610/93‑1 (MQ=255)

                                                                                   tCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTa      >  1:4222908/1‑93 (MQ=255)

                                                                                    ttttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  1:1715168/92‑1 (MQ=255)

                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  1:2727961/93‑1 (MQ=255)

                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  2:2957014/93‑1 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    <  2:4464382/93‑1 (MQ=255)

                                                                                      ttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAgg   >  2:1676886/1‑93 (MQ=255)

                                                                                       tttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:851870/93‑1 (MQ=255)

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

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: `‑`