BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,495	G→A	V22V (GTC→GTT)	smf ←	protein Smf

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,495	0	G	A	100.0%	1131.8 / NA	263	V22V (GTC→GTT)	smf	protein Smf
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (147/116); total (147/116)

Only 100 of 263 total aligned reads displayed.

AATGAAAGAAAGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCA  >  NC_000913/3433397‑3433591

aaTGAAAGAAAGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAAc                                                                                                >  1:5407272/1‑101 (MQ=255)

 aTGAAAGAAAGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAAcc                                                                                               <  2:5406942/101‑1 (MQ=255)

 aTGAAAGAAAGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAAcc                                                                                               >  1:5407275/1‑101 (MQ=255)

   gaaagaaaGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCat                                                                                             >  2:5407278/1‑101 (MQ=255)

    aaagaaagCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCata                                                                                            >  1:5407282/1‑101 (MQ=255)

    aaagaaagCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCata                                                                                            >  2:5407283/1‑101 (MQ=255)

       gaaagCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATatca                                                                                         >  1:5407285/1‑101 (MQ=255)

           gCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGc                                                                                     <  2:5406888/101‑1 (MQ=255)

            cGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGcc                                                                                    <  1:5406698/101‑1 (MQ=255)

            cGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGcc                                                                                    >  2:5407290/1‑101 (MQ=255)

             gTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCg                                                                                   >  1:5407292/1‑101 (MQ=255)

              ttGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGt                                                                                  <  1:5406854/101‑1 (MQ=255)

                 tgcctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACa                                                                               >  2:5407294/1‑101 (MQ=255)

                 tgcctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACa                                                                               <  1:5406975/101‑1 (MQ=255)

                 tgcctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACa                                                                               <  2:5406841/101‑1 (MQ=255)

                  gcctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACaa                                                                              >  1:5407299/1‑101 (MQ=255)

                  gcctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACaa                                                                              >  2:5407300/1‑101 (MQ=255)

                   cctgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAg                                                                             <  1:5407076/101‑1 (MQ=255)

                     tgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAgct                                                                           >  1:5407306/1‑101 (MQ=255)

                     tgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAgct                                                                           >  2:5407303/1‑101 (MQ=255)

                     tgccGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAgct                                                                           <  1:5406857/101‑1 (MQ=255)

                        cGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAgctgct                                                                        >  2:5407310/1‑101 (MQ=255)

                          cAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGa                                                                      >  1:5407313/1‑101 (MQ=255)

                          cAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGa                                                                      <  2:5406793/101‑1 (MQ=255)

                             tGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATAc                                                                   <  1:5406783/101‑1 (MQ=255)

                                aaGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCa                                                                <  1:5406804/101‑1 (MQ=255)

                                 aGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCAt                                                               >  1:5407318/1‑101 (MQ=255)

                                 aGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCAt                                                               <  1:5406693/101‑1 (MQ=255)

                                  gCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCAtt                                                              <  1:5406878/101‑1 (MQ=255)

                                  gCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCAtt                                                              >  2:5407319/1‑101 (MQ=255)

                                  gCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCAtt                                                              <  1:5406911/101‑1 (MQ=255)

                                   cccTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTa                                                             >  2:5407321/1‑101 (MQ=255)

                                   cccTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTa                                                             <  2:5406932/101‑1 (MQ=255)

                                    ccTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTaa                                                            <  2:5406605/101‑1 (MQ=255)

                                       gTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACg                                                         >  1:5407325/1‑101 (MQ=255)

                                        tttGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGc                                                        >  2:5407326/1‑101 (MQ=255)

                                         ttGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCa                                                       >  2:5407329/1‑101 (MQ=255)

                                          tgctgcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAg                                                      <  2:5406718/101‑1 (MQ=255)

                                              gcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGAt                                                                     >  2:5407335/1‑82 (MQ=255)

                                              gcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCaa                                                  >  1:5407333/1‑101 (MQ=255)

                                              gcAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCaa                                                  <  1:5406829/101‑1 (MQ=255)

                                               cAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCaaa                                                 >  1:5407336/1‑101 (MQ=255)

                                                 aTACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAgc                                                                            >  1:5407341/1‑72 (MQ=255)

                                                 aTACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAAtt                                               <  2:5406861/101‑1 (MQ=255)

                                                   aCAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTc                                             >  2:5407345/1‑101 (MQ=255)

                                                    cAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCt                                            >  1:5407346/1‑101 (MQ=255)

                                                     aaCCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTg                                           >  1:5407351/1‑101 (MQ=255)

                                                     aaCCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTg                                           >  2:5407349/1‑101 (MQ=255)

                                                     aaCCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTg                                           <  1:5406884/101‑1 (MQ=255)

                                                      aCCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGt                                          <  2:5406979/101‑1 (MQ=255)

                                                       ccGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAAtt                                               >  2:5407355/1‑95 (MQ=255)

                                                       ccGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTa                                         <  2:5406756/101‑1 (MQ=255)

                                                         gCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATc                                       <  1:5407016/101‑1 (MQ=255)

                                                         gCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATc                                       <  2:5406955/101‑1 (MQ=255)

                                                          cATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCg                                      >  1:5407358/1‑101 (MQ=255)

                                                          cATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCg                                      >  2:5407359/1‑101 (MQ=255)

                                                           aTCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGa                                     >  1:5407363/1‑101 (MQ=255)

                                                              aaTATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCa                                  <  2:5406920/101‑1 (MQ=255)

                                                              aaTATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCa                                  >  1:5407367/1‑101 (MQ=255)

                                                              aaTATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCa                                  <  1:5406921/101‑1 (MQ=255)

                                                              aaTATGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCa                                  >  2:5407365/1‑101 (MQ=255)

                                                               atatGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCAt                                 <  2:5407026/101‑1 (MQ=255)

                                                               atatGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCAt                                 >  1:5407370/1‑101 (MQ=255)

                                                                tatGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATc                                >  2:5407372/1‑101 (MQ=255)

                                                                tatGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATc                                >  1:5407374/1‑101 (MQ=255)

                                                                 atGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATcc                               >  1:5407377/1‑101 (MQ=255)

                                                                  tGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCt                              >  2:5407379/1‑101 (MQ=255)

                                                                  tGCGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCt                              <  1:5407049/101‑1 (MQ=255)

                                                                    cGACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTa                            <  1:5406711/101‑1 (MQ=255)

                                                                     gACTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTAt                           >  1:5407383/1‑101 (MQ=255)

                                                                      aCTGTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTAtc                          >  1:5407386/1‑101 (MQ=255)

                                                                        tGTTTTGCCACCCAGTGAGCGCTACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTAtctc                        >  1:5407389/1‑101 (MQ=255)

                                                                         gTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTcc                       >  2:5407391/1‑101 (MQ=255)

                                                                         gTTTTGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTcc                       <  2:5407306/101‑1 (MQ=255)

                                                                           tttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCt                     <  1:5406905/101‑1 (MQ=255)

                                                                           tttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCt                     <  2:5406810/101‑1 (MQ=255)

                                                                           tttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCt                     <  1:5406958/101‑1 (MQ=255)

                                                                           tttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCt                     >  1:5407394/1‑101 (MQ=255)

                                                                            ttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTg                    <  1:5406930/101‑1 (MQ=255)

                                                                            ttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTg                    >  2:5407396/1‑101 (MQ=255)

                                                                            ttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTg                    <  1:5407000/101‑1 (MQ=255)

                                                                            ttGCCACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTg                    <  1:5407053/101‑1 (MQ=255)

                                                                                cACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCAt                <  2:5406869/101‑1 (MQ=255)

                                                                                cACCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCAt                >  2:5407399/1‑101 (MQ=255)

                                                                                 aCCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTg                                           <  1:5406817/73‑1 (MQ=255)

                                                                                 aCCCAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATa               >  1:5407401/1‑101 (MQ=255)

                                                                                  cccAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATaa              <  2:5407140/101‑1 (MQ=255)

                                                                                   ccAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAag             >  2:5407405/1‑101 (MQ=255)

                                                                                   ccAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAag             <  1:5406948/101‑1 (MQ=255)

                                                                                   ccAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAag             <  1:5406903/101‑1 (MQ=255)

                                                                                    cAGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAagc            <  1:5406890/101‑1 (MQ=255)

                                                                                     aGTGAGCGATACGAACCATATCATCGCCGTACAAGCTGCNGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAagca           >  1:5407407/1‑101 (MQ=255)

                                                                                        gAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATa                                                                    <  1:5406734/41‑1 (MQ=255)

                                                                                        gAGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGcc        >  2:5407410/1‑101 (MQ=255)

                                                                                         aGCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCt       >  2:5407413/1‑101 (MQ=255)

                                                                                          gCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCtt      <  1:5406993/101‑1 (MQ=255)

                                                                                          gCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCtt      >  1:5407414/1‑101 (MQ=255)

                                                                                          gCGATACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCtt      <  2:5406778/101‑1 (MQ=255)

                                                                                             aTACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGc   >  2:5407417/1‑101 (MQ=255)

                                                                                              tACGAACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCa  >  1:5407420/1‑101 (MQ=255)

AATGAAAGAAAGCGTTGTGCCTGCCGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCA  >  NC_000913/3433397‑3433591

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`