BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,463,652	C→G	6.7%	G353A (GGC→GCC)	treC ←	trehalose‑6‑P hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,463,652	0	C	G	6.7%	368.8 / 15.1	105	G353A (GGC→GCC)	treC	trehalose‑6‑P hydrolase
Reads supporting (aligned to +/- strand): ref base C (49/49); new base G (3/4); total (52/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

Only 100 of 113 total aligned reads displayed.

AATAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAA  >  NC_000913/4463505‑4463800

aaTAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGcatgcc                                                                                                                                                     >  1:80110/1‑149 (MQ=255)

aaTAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGcatgcc                                                                                                                                                     <  2:80110/149‑1 (MQ=255)

 aTAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGcatgcca                                                                                                                                                    <  2:37301/149‑1 (MQ=255)

      tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                               >  4:7227/1‑149 (MQ=255)

      tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                               >  3:651035/1‑149 (MQ=255)

      tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                               >  4:188159/1‑149 (MQ=255)

      tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                               >  1:585744/1‑149 (MQ=255)

       cgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagg                                                                                                                                              <  3:130149/149‑2 (MQ=255)

        gtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcaggg                                                                                                                                             <  1:386786/149‑3 (MQ=255)

        gtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcagca                                                                                                                                             >  2:278299/1‑149 (MQ=255)

         tcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCAc                                                                                                                                            >  4:184576/1‑149 (MQ=255)

         tcgGCATCACGCCCATCGTTGCGCAGCGCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCAc                                                                                                                                            >  4:680431/1‑149 (MQ=255)

           ggCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCa                                                                                                                                          >  4:352040/1‑149 (MQ=255)

            gCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagggaa                                                                                                                                           >  1:273491/1‑143 (MQ=255)

            gCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagggaa                                                                                                                                           <  2:273491/147‑5 (MQ=255)

             cATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATc                                                                                                                                        <  2:334821/149‑1 (MQ=255)

              aTCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCg                                                                                                                                       >  3:491308/1‑149 (MQ=255)

               tCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGc                                                                                                                                      >  2:381793/1‑149 (MQ=255)

               tCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGc                                                                                                                                      <  2:577437/149‑1 (MQ=255)

                     ccATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAt                                                                                                                                >  4:74689/1‑149 (MQ=255)

                     ccATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAt                                                                                                                                >  4:50553/1‑149 (MQ=255)

                      cATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAtt                                                                                                                               >  1:465515/1‑149 (MQ=255)

                       aTCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttt                                                                                                                              >  3:685594/1‑149 (MQ=255)

                        tCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAtttt                                                                                                                             >  1:294952/1‑149 (MQ=255)

                         cGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                            <  4:251654/149‑1 (MQ=255)

                         cGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                            <  4:609674/149‑1 (MQ=255)

                           ttGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcaggga                                                                                                                                            <  4:633436/131‑4 (MQ=255)

                           ttGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcaggga                                                                                                                                            >  3:633436/1‑128 (MQ=255)

                             gcgcAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCg                                                                                                                        >  2:476416/1‑149 (MQ=255)

                              cgcAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGc                                                                                                                       <  3:193208/149‑1 (MQ=255)

                               gcAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCa                                                                                                                      >  3:6175/1‑149 (MQ=255)

                                cAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAg                                                                                                                     >  1:449558/1‑149 (MQ=255)

                                  gCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAACACCATCGCCAGCATTTTTGCCGCAGGc                                                                                                                   >  2:224315/1‑149 (MQ=255)

                                    tCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                            <  4:139230/138‑1 (MQ=255)

                                    tCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                            >  3:139230/1‑138 (MQ=255)

                                    tCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCAc                                                                                                                 >  1:371934/1‑149 (MQ=255)

                                      ggCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCAcgc                                                                                                               <  2:294952/149‑1 (MQ=255)

                                       gCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCAcgcg                                                                                                              <  3:352040/149‑1 (MQ=255)

                                         aaaCATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGt                                                                                                            >  2:641934/1‑149 (MQ=255)

                                         aaaCATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGt                                                                                                            >  4:42196/1‑149 (MQ=255)

                                          aaCATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTa                                                                                                           >  2:248928/1‑149 (MQ=255)

                                             atatTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTc                                                                                                        >  3:218426/1‑149 (MQ=255)

                                             atatTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTc                                                                                                        <  1:211707/149‑1 (MQ=255)

                                                ttGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                     >  4:605154/1‑149 (MQ=255)

                                                ttGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                     >  3:157299/1‑149 (MQ=255)

                                                     gCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTtcat                                                                                                <  1:278299/149‑1 (MQ=255)

                                                          ctACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCa                                                                                           >  1:285184/1‑149 (MQ=255)

                                                            aCGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCaaa                                                                                         <  2:640946/149‑1 (MQ=255)

                                                            aCGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCaaa                                                                                         <  4:171070/149‑1 (MQ=255)

                                                                cgcgATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCga                                                                                     <  1:476416/149‑1 (MQ=255)

                                                                   gATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAtttt                                                                                                                             <  1:429728/106‑1 (MQ=255)

                                                                   gATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGaa                                                                                  <  2:629853/149‑1 (MQ=255)

                                                                     tAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAAc                                                                                >  1:65196/1‑149 (MQ=255)

                                                                          aGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgc                                                                           >  4:472429/1‑149 (MQ=255)

                                                                            tAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgcgc                                                                         <  2:585744/149‑1 (MQ=255)

                                                                             aaTGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGcgcg                                                                        <  2:449558/149‑1 (MQ=255)

                                                                               tGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGc                                                                      <  3:74689/149‑1 (MQ=255)

                                                                                gcgcGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCt                                                                     >  4:145868/1‑149 (MQ=255)

                                                                                       aaaTGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTg                                                              <  2:371934/149‑1 (MQ=255)

                                                                                          tGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGtt                                                           <  3:50553/149‑1 (MQ=255)

                                                                                             gggTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTaca                                                        <  3:7227/149‑1 (MQ=255)

                                                                                                   gTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGc                                                                      <  2:505407/129‑1 (MQ=255)

                                                                                                   gTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGc                                                                      >  1:505407/1‑129 (MQ=255)

                                                                                                      aTGCCAATCTCTTCGCCCTGTTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                               >  1:550571/1‑149 (MQ=255)

                                                                                                      aTGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                               <  4:252691/149‑1 (MQ=255)

                                                                                                      aTGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                               <  3:322838/149‑1 (MQ=255)

                                                                                                      aTGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                               <  4:617586/149‑1 (MQ=255)

                                                                                                       tGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAg                                              <  3:98759/149‑1 (MQ=255)

                                                                                                        gCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAgg                                             >  4:133940/1‑149 (MQ=255)

                                                                                                        gCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAgg                                             <  2:285184/149‑1 (MQ=255)

                                                                                                          cAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTTACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCa                                           >  2:569691/1‑149 (MQ=255)

                                                                                                          cAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGaa                                                  >  4:664339/1‑142 (MQ=255)

                                                                                                          cAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCa                                           >  3:61249/1‑149 (MQ=255)

                                                                                                           aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCa                                           >  3:651003/1‑148 (MQ=255)

                                                                                                           aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCa                                           <  4:651003/148‑1 (MQ=255)

                                                                                                           aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCAt                                          <  2:628312/149‑1 (MQ=255)

                                                                                                           aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCAt                                          >  2:466605/1‑149 (MQ=255)

                                                                                                           aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGACGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGt                                                               <  3:401748/128‑1 (MQ=255)

                                                                                                             tctcTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTc                                        >  1:181119/1‑149 (MQ=255)

                                                                                                             tctcTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTc                                        <  3:145868/149‑1 (MQ=255)

                                                                                                                ctTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCAt                                     <  4:218426/149‑1 (MQ=255)

                                                                                                                 ttCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATg                                    <  3:605154/149‑1 (MQ=255)

                                                                                                                    gCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTa                                 <  3:522235/149‑1 (MQ=255)

                                                                                                                     cccTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTAc                                <  1:420118/149‑1 (MQ=255)

                                                                                                                       cTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGt                              >  3:494020/1‑149 (MQ=255)

                                                                                                                        tGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCgag                                                                                    <  3:555294/94‑1 (MQ=255)

                                                                                                                        tGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCgag                                                                                    >  4:555294/1‑94 (MQ=255)

                                                                                                                          gTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTtgt                           <  2:97296/149‑1 (MQ=255)

                                                                                                                                ataCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTc                     <  4:494020/149‑1 (MQ=255)

                                                                                                                                 taCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGc                                            >  2:417922/1‑125 (MQ=255)

                                                                                                                                 taCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGc                                            <  1:417922/125‑1 (MQ=255)

                                                                                                                                 taCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTcc                    >  3:124432/1‑149 (MQ=255)

                                                                                                                                    ggCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCATAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttg                 <  1:641934/149‑1 (MQ=255)

                                                                                                                                    ggCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttg                 <  3:42196/149‑1 (MQ=255)

                                                                                                                                    ggCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttg                 <  3:680431/149‑1 (MQ=255)

                                                                                                                                      cGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttgtt               <  3:184576/149‑1 (MQ=255)

                                                                                                                                               catgccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGc      <  4:491308/149‑1 (MQ=255)

                                                                                                                                                  gccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATc                                                                 >  3:143700/1‑87 (MQ=255)

                                                                                                                                                   ccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGaa  >  2:573924/1‑149 (MQ=255)

                                                                                                                                                   ccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGaa  >  3:308850/1‑149 (MQ=255)

AATAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAA  >  NC_000913/4463505‑4463800

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