BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	422,540	C→A	5.6%	P8Q (CCG→CAG)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	422,540	0	C	A	5.6%	375.2 / 9.8	107	P8Q (CCG→CAG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand): ref base C (45/56); new base A (2/4); total (47/60)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

Only 100 of 109 total aligned reads displayed.

TCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGTCAGCC  >  NC_000913/422393‑422687

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)

       ccGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccat                                                                                                                                             <  3:684488/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)

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

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

                      gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAGGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCAACTGCCGGTGCCCCCATTTGTTAAAcaaa                                                                                                                              <  3:401408/149‑1 (MQ=255)

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

                          tCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagca                                                                                                                           >  1:582120/1‑148 (MQ=255)

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

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

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

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

                                         ctccCCCAAGTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa                                                                                                                    <  1:748601/140‑1 (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)

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

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

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

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

                                              ccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTaca                                                                                                       <  2:736810/148‑1 (MQ=255)

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

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

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

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

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

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

                                                         tCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGAc                                                                                            >  1:7378/1‑148 (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)

                                                                         ccTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCg                                                                            >  2:286720/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                                                                <  3:726420/147‑1 (MQ=255)

                                                                                      aTGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTAtg                                                                >  4:726420/1‑147 (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:228626/149‑1 (MQ=255)

                                                                                                gATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGaa                                                    <  4:256781/149‑1 (MQ=255)

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

                                                                                                    aGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                                >  4:860853/1‑149 (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                                                                                                                                 >  3:754300/1‑57 (MQ=255)

                                                                                                               acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAc                                                                                                                                 <  4:754300/57‑1 (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                                 <  4:64617/140‑1 (MQ=255)

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

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

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

                                                                                                                             atgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGa                        <  2:230458/148‑1 (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                >  4:713655/1‑149 (MQ=255)

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

                                                                                                                                      gCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAagca              <  4:676298/149‑1 (MQ=255)

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

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

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

                                                                                                                                                  ccGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGTCAGcc  <  3:860853/149‑1 (MQ=255)

TCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAAGCAGCGCAGTCAGCC  >  NC_000913/422393‑422687

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