BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	422,537	T→G	5.6%	L7R (CTG→CGG)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	422,537	0	T	G	5.6%	378.6 / 13.3	108	L7R (CTG→CGG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand): ref base T (46/56); new base G (2/4); total (48/60)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 108 total aligned reads displayed.

AATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGT  >  NC_000913/422391‑422682

aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGc                                                                                                                                                 >  4:365140/1‑149 (MQ=255)

  tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCg                                                                                                                                               >  4:183521/1‑149 (MQ=255)

    cTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCgg                                                                                                                                              >  2:361958/1‑148 (MQ=255)

        tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcc                                                                                                                                          >  2:381328/1‑148 (MQ=255)

         ccGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccat                                                                                                                                        <  4:500199/149‑3 (MQ=255)

          cGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGCTCCTGCATTAGGCTATGGCAAGGTGGTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGccccc                                                                                                                                       <  4:388422/149‑1 (MQ=255)

            gCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCTGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAt                                                                                                                                     >  3:659122/1‑149 (MQ=255)

            gCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCa                                                                                                                                      <  2:346631/148‑1 (MQ=255)

              tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAtt                                                                                                                                    <  1:381328/148‑1 (MQ=255)

                 tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                                                                                                                                <  3:284739/149‑1 (MQ=255)

                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaaa                                                                                                                             <  3:877884/149‑1 (MQ=255)

                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaaa                                                                                                                             >  4:743317/1‑149 (MQ=255)

                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                                                                                                                              <  2:641665/148‑1 (MQ=255)

                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGcca                                                                                                                                         <  1:164749/137‑2 (MQ=255)

                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGcca                                                                                                                                         >  2:164749/1‑136 (MQ=255)

                     aaaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaaa                                                                                                                             >  2:887133/1‑148 (MQ=255)

                        gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa                                                                                                                         <  3:365140/149‑1 (MQ=255)

                            tCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa                                                                                                                     >  3:664692/1‑149 (MQ=255)

                                 ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccat                                                                                                                                        <  3:596420/125‑3 (MQ=255)

                                 ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccat                                                                                                                                        >  4:596420/1‑123 (MQ=255)

                                   ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa                                                                                                              <  4:615161/149‑1 (MQ=255)

                                           ctccCCCAAGTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa                                                                                                               <  1:748601/140‑1 (MQ=255)

                                           ctccCCCAAGTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa                                                                                                               >  2:748601/1‑140 (MQ=255)

                                            tccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt                                                                                                     >  4:819749/1‑149 (MQ=255)

                                             cccccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt                                                                                                     <  2:903437/148‑1 (MQ=255)

                                              ccccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTac                                                                                                   >  4:184718/1‑149 (MQ=255)

                                               cccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGTCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTaca                                                                                                  >  4:763680/1‑149 (MQ=255)

                                               cccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc                                                                                                          <  3:427607/141‑1 (MQ=255)

                                               cccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc                                                                                                          >  4:427607/1‑141 (MQ=255)

                                                ccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTacac                                                                                                 >  3:681371/1‑149 (MQ=255)

                                                ccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTacac                                                                                                 <  3:446052/149‑1 (MQ=255)

                                                ccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTacac                                                                                                 <  3:184718/149‑1 (MQ=255)

                                                ccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTaca                                                                                                  >  1:893985/1‑148 (MQ=255)

                                                   aaTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACtgt                                                                                              <  3:371451/149‑1 (MQ=255)

                                                    atatCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACtgt                                                                                              <  1:910048/148‑1 (MQ=255)

                                                          ttCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGa                                                                                        >  1:671754/1‑148 (MQ=255)

                                                          ttCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGa                                                                                        <  1:887133/148‑1 (MQ=255)

                                                          ttCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGAc                                                                                       >  4:549116/1‑149 (MQ=255)

                                                           tCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGAc                                                                                       >  1:7378/1‑148 (MQ=255)

                                                           tCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGAc                                                                                       <  1:650199/148‑1 (MQ=255)

                                                           tCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGCGGCTGAc                                                                                       >  1:122977/1‑148 (MQ=255)

                                                              gatgatGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACggg                                                                                    <  1:293747/148‑1 (MQ=255)

                                                              gatgatGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGc                                                                                   <  4:664692/149‑1 (MQ=255)

                                                                 gatgaGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGaa                                                                                <  3:454345/149‑1 (MQ=255)

                                                                  atgaGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGaa                                                                                >  1:353045/1‑148 (MQ=255)

                                                                    gaGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGa                                                                              <  1:336680/148‑1 (MQ=255)

                                                                      gTGATCCTGCCTTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTAGCTGACGGGCGAAGAccc                                                                           <  4:436901/149‑1 (MQ=255)

                                                                      gTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGAccc                                                                           <  4:728833/149‑1 (MQ=255)

                                                                       tGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGAccc                                                                           <  2:134848/148‑1 (MQ=255)

                                                                          tCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCg                                                                       >  4:79255/1‑149 (MQ=255)

                                                                           ccTGCATTAGGCTATGGCAAGGTGATCAGTTTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCg                                                                       >  1:230458/1‑148 (MQ=255)

                                                                             tGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCa                                                                     >  2:638860/1‑148 (MQ=255)

                                                                                aTTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAgcgc                                                                 >  3:44499/1‑149 (MQ=255)

                                                                                aTTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAgcg                                                                  <  2:319189/148‑1 (MQ=255)

                                                                                  tAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAt                                                               <  3:379847/149‑1 (MQ=255)

                                                                                  tAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAC‑AAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAtt                                                              <  4:659122/149‑1 (MQ=255)

                                                                                      cTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtg                                                           <  3:285624/149‑1 (MQ=255)

                                                                                       tATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtgc                                                          >  3:228626/1‑149 (MQ=255)

                                                                                        aTGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtg                                                           >  4:726420/1‑147 (MQ=255)

                                                                                        aTGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtg                                                           <  3:726420/147‑1 (MQ=255)

                                                                                           gCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtgctgc                                                       >  2:457959/1‑148 (MQ=255)

                                                                                           gCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCg                                                      >  3:646133/1‑149 (MQ=255)

                                                                                             aaGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAt                                                                                                      <  4:818021/99‑1 (MQ=255)

                                                                                             aaGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAt                                                                                                      >  3:818021/1‑99 (MQ=255)

                                                                                               ggTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTAc                                                   <  1:361958/148‑1 (MQ=255)

                                                                                               ggTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACa                                                  <  4:681371/149‑1 (MQ=255)

                                                                                               ggTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACa                                                  <  4:228626/149‑1 (MQ=255)

                                                                                                     cAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAaca                                             >  2:332632/1‑148 (MQ=255)

                                                                                                      aGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                           >  4:860853/1‑149 (MQ=255)

                                                                                                       gATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                           <  1:457959/148‑1 (MQ=255)

                                                                                                        aTTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGa                                          <  1:332632/148‑1 (MQ=255)

                                                                                                          tttCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGATGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATa                                        <  2:671754/148‑1 (MQ=255)

                                                                                                                cacaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa                                                                                                              <  1:894940/72‑1 (MQ=255)

                                                                                                                cacaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa                                                                                                              >  2:894940/1‑72 (MQ=255)

                                                                                                                 acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAc                                                                                                                            <  4:754300/57‑1 (MQ=255)

                                                                                                                 acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAc                                                                                                                            >  3:754300/1‑57 (MQ=255)

                                                                                                                   aGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGaaa                               <  2:353045/148‑1 (MQ=255)

                                                                                                                      ggAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGa                                                                              <  4:867266/98‑1 (MQ=255)

                                                                                                                      ggAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGa                                                                              >  3:867266/1‑98 (MQ=255)

                                                                                                                         aTTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAt                                                            >  2:304434/1‑113 (MQ=255)

                                                                                                                         aTTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAt                                                            <  1:304434/113‑1 (MQ=255)

                                                                                                                          ttATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGt                        <  2:7378/148‑1 (MQ=255)

                                                                                                                           tATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTAc                      >  4:245069/1‑149 (MQ=255)

                                                                                                                              gatgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGt                            >  3:64617/1‑140 (MQ=255)

                                                                                                                              gatgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGt                            <  4:64617/140‑1 (MQ=255)

                                                                                                                              gatgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGa                   <  3:819749/149‑1 (MQ=255)

                                                                                                                               ctgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACGAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt                  <  3:79255/148‑1 (MQ=255)

                                                                                                                               atgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGa                   <  2:230458/148‑1 (MQ=255)

                                                                                                                               atgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt                  >  3:813356/1‑149 (MQ=255)

                                                                                                                                tgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt                  <  1:638860/148‑1 (MQ=255)

                                                                                                                                tgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATg                 <  3:549116/149‑1 (MQ=255)

                                                                                                                                     aaTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATaa            <  4:44499/149‑1 (MQ=255)

                                                                                                                                     aaTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATa             >  2:422437/1‑148 (MQ=255)

                                                                                                                                      aTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAag           >  3:676298/1‑149 (MQ=255)

                                                                                                                                      aTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAag           >  4:713655/1‑149 (MQ=255)

                                                                                                                                         cATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAagca         <  1:286720/148‑1 (MQ=255)

                                                                                                                                           tGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtgct                                                         <  4:584806/98‑1 (MQ=255)

                                                                                                                                           tGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtgct                                                         >  3:584806/1‑98 (MQ=255)

                                                                                                                                             gCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAgcgc     <  2:582120/148‑1 (MQ=255)

                                                                                                                                                ccTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGt  >  1:285790/1‑148 (MQ=255)

AATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGT  >  NC_000913/422391‑422682

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