BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,961	Δ1 bp	coding (1386/1458 nt)	b3876 ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,959	0	C	.	100.0%	949.7 / NA	206	coding (1388/1458 nt)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (104/102); total (104/102)

Only 100 of 222 total aligned reads displayed.

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  2:5186560/1‑93 (MQ=255)

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  2:4777963/1‑93 (MQ=255)

  tttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                        >  2:866798/1‑92 (MQ=255)

   ttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgt                                                                                       >  1:4140773/1‑91 (MQ=255)

      ataataAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttt                                                                                   <  1:6276500/93‑1 (MQ=255)

      ataataAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttt                                                                                   >  2:5163818/1‑93 (MQ=255)

       taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                  >  2:1305831/1‑93 (MQ=255)

       taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                  >  1:1571871/1‑93 (MQ=255)

       taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                  >  1:43029/1‑93 (MQ=255)

       taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                  >  1:5145254/1‑93 (MQ=255)

        aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                 >  2:4806089/1‑93 (MQ=255)

         ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttc                                                                                >  1:487856/1‑93 (MQ=255)

          taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttct                                                                               >  2:5560479/1‑93 (MQ=255)

           aaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttctt                                                                              <  2:5205875/93‑1 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:212886/93‑1 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            >  1:6100573/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  2:2672500/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  1:2831422/1‑93 (MQ=255)

                cATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTTTAACGTTACGGAAG‑CCGTTTTTTCTTCTGTc                                                                         >  1:4858326/1‑93 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGc                                                                       >  1:5256587/1‑93 (MQ=255)

                    ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTATGGAAG‑CCGTCTTTTCTTCTGTCGCCa                                                                     >  2:3134692/1‑93 (MQ=255)

                    ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCa                                                                     <  2:3529320/93‑1 (MQ=255)

                    ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCa                                                                     <  1:105491/93‑1 (MQ=255)

                     tttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAt                                                                    >  1:5798890/1‑93 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             <  1:4957116/93‑1 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             <  1:4786938/93‑1 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             >  2:3723291/1‑93 (MQ=255)

                             tAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAt                                                            >  1:1076359/1‑93 (MQ=255)

                              aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                           <  1:1876088/93‑1 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          >  1:3800363/1‑93 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          >  2:648270/1‑93 (MQ=255)

                                ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                         <  2:3649496/93‑1 (MQ=255)

                                ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                         <  2:5760870/93‑1 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      >  2:1958024/1‑93 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      >  2:5525066/1‑93 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      <  1:6153059/93‑1 (MQ=255)

                                    gccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTTCg                                                     >  2:3351667/1‑93 (MQ=255)

                                     ccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgt                                                    <  1:3907556/93‑1 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   <  2:4683445/93‑1 (MQ=255)

                                       gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                  >  2:1275063/1‑93 (MQ=255)

                                        ccacccacCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttgg                                                 >  2:2599114/1‑93 (MQ=255)

                                        ccacccacCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttgg                                                 <  2:5976182/93‑1 (MQ=255)

                                        ccacccacCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttgg                                                 <  1:1902894/93‑1 (MQ=255)

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                <  1:5015006/93‑1 (MQ=255)

                                           cccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggttg                                              <  1:711206/93‑1 (MQ=255)

                                             caccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAt                                            >  2:2949098/1‑93 (MQ=255)

                                              accaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAtt                                           >  2:4189056/1‑93 (MQ=255)

                                               ccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAttt                                          <  1:2656026/93‑1 (MQ=255)

                                               ccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAttt                                          <  1:1441110/93‑1 (MQ=255)

                                               ccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAttt                                          <  1:3121021/93‑1 (MQ=255)

                                                caGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTc                                         <  2:3957060/93‑1 (MQ=255)

                                                 aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTcc                                        <  2:289212/93‑1 (MQ=255)

                                                 aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTcc                                        >  1:2912468/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       >  1:4009428/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGACTATTGACATTAATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       <  2:3240629/93‑1 (MQ=255)

                                                   gATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTc                                      <  2:5497066/93‑1 (MQ=255)

                                                   gATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTc                                      >  2:4924178/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  2:5370949/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  2:1093873/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  2:1588195/1‑93 (MQ=255)

                                                     tGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTCTTGCGTTGGTTGATTTCCTCaa                                    >  1:2014982/1‑93 (MQ=255)

                                                      ggAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                   <  1:959143/93‑1 (MQ=255)

                                                       gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                  <  1:5116336/93‑1 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 <  2:922193/93‑1 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 >  2:2517555/1‑93 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 <  1:1848795/93‑1 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 <  1:3409012/93‑1 (MQ=255)

                                                          aTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATc                               <  2:1789114/93‑1 (MQ=255)

                                                           tAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCa                              >  2:1440644/1‑93 (MQ=255)

                                                            aGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                             <  1:1928439/93‑1 (MQ=255)

                                                             gCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                            >  1:1308346/1‑93 (MQ=255)

                                                             gCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                            <  1:2953822/93‑1 (MQ=255)

                                                              cGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATg                           <  2:2226024/93‑1 (MQ=255)

                                                              cGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATg                           >  1:3047131/1‑93 (MQ=255)

                                                              cGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATg                           >  1:352844/1‑93 (MQ=255)

                                                              cGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATg                           >  2:1371346/1‑93 (MQ=255)

                                                                  tATTGACATTTATAACGTTCCGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                       <  1:5942250/93‑1 (MQ=255)

                                                                  tATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                       <  2:3999871/93‑1 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      <  2:1229230/93‑1 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     <  1:3346875/93‑1 (MQ=255)

                                                                     tGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                    >  1:868720/1‑93 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  1:2236043/1‑93 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  1:3653350/1‑93 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  2:4150087/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  1:5901756/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  1:6062164/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  1:5307914/93‑1 (MQ=255)

                                                                            tataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCa             <  1:2474317/93‑1 (MQ=255)

                                                                             ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTTCTCCATAATCAATGCAAACCGTTTTTCAt            >  2:561956/1‑93 (MQ=255)

                                                                              taACGTTTCGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:4904557/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:2146986/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:3281398/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           >  1:416712/1‑93 (MQ=255)

                                                                               aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTg          >  1:6189029/1‑93 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         <  1:4554079/93‑1 (MQ=255)

                                                                                 cGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAg        <  1:1834655/93‑1 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtg      >  2:1012794/1‑93 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtg      >  2:4129648/1‑93 (MQ=255)

                                                                                    tACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt     >  2:4436686/1‑93 (MQ=255)

                                                                                       ggAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTa  >  1:3648/1‑93 (MQ=255)

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

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