BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,927	C→A	6.5%	T45T (ACG→ACT)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,927	0	C	A	6.5%	428.5 / 14.9	124	T45T (ACG→ACT)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base C (58/58); new base A (6/2); total (64/60)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 126 total aligned reads displayed.

GAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCA  >  NC_000913/858780‑859073

gAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTcgcg                                                                                                                                                   >  3:940256/1‑149 (MQ=255)

  aaTGCCCTTTTATCACCGCACAACCGCGACGTGTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTACTTTGAATGCGTTGATTTCAGCCACTCGCGt                                                                                                                                                  >  1:277239/1‑148 (MQ=255)

   aTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTg                                                                                                                                                 <  1:830685/148‑1 (MQ=255)

   aTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTg                                                                                                                                                 >  2:935232/1‑148 (MQ=255)

     ggCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGc                                                                                                                                               >  1:201215/1‑148 (MQ=255)

      gCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCa                                                                                                                                              >  2:523226/1‑148 (MQ=255)

       ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAAt                                                                                                                                            >  3:372307/1‑149 (MQ=255)

       ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAAt                                                                                                                                            >  3:262440/1‑149 (MQ=255)

         ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                           <  1:60624/148‑1 (MQ=255)

         ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                           >  1:909707/1‑148 (MQ=255)

         ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                           >  2:177877/1‑148 (MQ=255)

             aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                      >  3:49907/1‑149 (MQ=255)

             aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                      >  4:883787/1‑149 (MQ=255)

               cACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGc                                                                                                                                    >  4:64814/1‑149 (MQ=255)

                       aaCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTg                                                                                                                             <  1:807094/148‑1 (MQ=255)

                        aCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAt                                                                                                                           <  3:391705/149‑1 (MQ=255)

                         ccGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAt                                                                                                                           >  1:861822/1‑148 (MQ=255)

                            cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                        <  1:935232/148‑1 (MQ=255)

                             gACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAg                                                                                                                      <  3:557186/149‑1 (MQ=255)

                               cGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGa                                                                                                                     <  2:379680/148‑1 (MQ=255)

                                 tttAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATATCATCGCTTAGCTGATTAAAGATc                                                                                                                   >  1:102070/1‑148 (MQ=255)

                                   tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                                <  3:599164/149‑1 (MQ=255)

                                   tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctg                                                                                                                 <  2:34016/148‑1 (MQ=255)

                                     aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgct                                                                                                               >  2:811880/1‑148 (MQ=255)

                                       cTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                             >  1:524168/1‑148 (MQ=255)

                                       cTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                             >  2:694582/1‑148 (MQ=255)

                                       cTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                             <  1:544310/148‑1 (MQ=255)

                                         gTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACg                                                                                                          >  3:608094/1‑149 (MQ=255)

                                           tCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGt                                                                                                         >  2:596913/1‑148 (MQ=255)

                                           tCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGt                                                                                                         >  1:409373/1‑148 (MQ=255)

                                             gcgcTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGc                                                                                                      >  3:573425/1‑149 (MQ=255)

                                                  gctCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGc                                                                                                 <  3:499034/149‑1 (MQ=255)

                                                     ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                              <  3:757153/149‑1 (MQ=255)

                                                     ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                               >  1:101230/1‑148 (MQ=255)

                                                     ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                               >  1:258093/1‑148 (MQ=255)

                                                       gCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgcg                                                                                             <  2:409373/148‑1 (MQ=255)

                                                         aGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGccc                                                                                          >  4:755750/1‑149 (MQ=255)

                                                         aGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGcc                                                                                           <  1:523226/148‑1 (MQ=255)

                                                          gTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGccc                                                                                          >  2:575268/1‑148 (MQ=255)

                                                           tAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTAGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCg                                                                                         >  2:376466/1‑148 (MQ=255)

                                                              cATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgc                                                                                     <  4:254519/149‑1 (MQ=255)

                                                              cATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAg                                                                                      >  2:483958/1‑148 (MQ=255)

                                                               aTGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgcg                                                                                    >  4:332650/1‑149 (MQ=255)

                                                               aTGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgcg                                                                                    >  4:596802/1‑149 (MQ=255)

                                                               aTGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgc                                                                                     >  2:234691/1‑148 (MQ=255)

                                                               aTGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgc                                                                                     >  1:1956/1‑148 (MQ=255)

                                                                tcaCCTGCTTTGATATCGGCATAAGACAGCACCGTCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgcgc                                                                                   <  3:778149/147‑1 (MQ=255)

                                                                tGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAgcg                                                                                    <  1:503055/148‑1 (MQ=255)

                                                                 gACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGGGCGCAgctt                                                                                   <  2:909707/148‑3 (MQ=255)

                                                                     tGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGc                                                                               >  1:607492/1‑148 (MQ=255)

                                                                     tGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCt                                                                              >  3:276943/1‑149 (MQ=255)

                                                                        tttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                           >  4:507903/1‑149 (MQ=255)

                                                                         ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                                 <  1:254736/142‑1 (MQ=255)

                                                                         ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                                 >  2:254736/1‑142 (MQ=255)

                                                                         ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                           >  1:45591/1‑148 (MQ=255)

                                                                         ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCAt                                                                          <  4:358997/149‑1 (MQ=255)

                                                                         ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATATCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                           <  2:102070/148‑1 (MQ=255)

                                                                          tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                         <  4:327123/149‑1 (MQ=255)

                                                                           gATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCt                                                                        <  3:220191/149‑1 (MQ=255)

                                                                             tatCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCAt                                                                          >  4:839034/1‑145 (MQ=255)

                                                                              atCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                      <  2:449074/148‑1 (MQ=255)

                                                                                cGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACg                                                                    >  1:215781/1‑148 (MQ=255)

                                                                                 ggCATAAGACAGCACCGTCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCCGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGTCCGCAGCGCGTGCTTCATCTGACGGa                                                                  <  4:711156/149‑1 (MQ=255)

                                                                                   cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                                 >  1:205440/1‑148 (MQ=255)

                                                                                     tAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgc                                                              <  3:507903/149‑1 (MQ=255)

                                                                                        gACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                           >  3:204418/1‑149 (MQ=255)

                                                                                         aCAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                           >  2:870261/1‑148 (MQ=255)

                                                                                           aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCttt                                                        <  4:737596/149‑1 (MQ=255)

                                                                                           aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCtt                                                         <  2:258093/148‑1 (MQ=255)

                                                                                              aCCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTtgtg                                                     <  4:372307/149‑1 (MQ=255)

                                                                                               ccGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGccc                                                                                          >  1:627628/1‑111 (MQ=255)

                                                                                               ccGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGccc                                                                                          <  2:627628/111‑1 (MQ=255)

                                                                                                cGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGa                                                    <  1:234691/148‑1 (MQ=255)

                                                                                                   ccAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                                  <  1:417607/131‑1 (MQ=255)

                                                                                                         ggCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggcaggca                                          >  4:169117/1‑149 (MQ=255)

                                                                                                          gCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggcaggca                                          <  2:663910/148‑1 (MQ=255)

                                                                                                           cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCaa                                         <  1:956987/148‑1 (MQ=255)

                                                                                                             tCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCCGGCAAtg                                       >  1:108056/1‑148 (MQ=255)

                                                                                                               gCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAtgtg                                     <  1:452807/148‑1 (MQ=255)

                                                                                                                  cTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGc                                 >  3:415763/1‑149 (MQ=255)

                                                                                                                  cTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGc                                 <  3:415757/149‑1 (MQ=255)

                                                                                                                    gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTg                               <  3:169117/149‑1 (MQ=255)

                                                                                                                       cTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                           >  3:314669/1‑118 (MQ=255)

                                                                                                                       cTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                           <  4:314669/118‑1 (MQ=255)

                                                                                                                          tGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                         <  3:767212/149‑1 (MQ=255)

                                                                                                                                cGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCATGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                   <  4:566613/149‑1 (MQ=255)

                                                                                                                                cGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                   <  4:608094/149‑1 (MQ=255)

                                                                                                                                 gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg                  <  3:943035/149‑1 (MQ=255)

                                                                                                                                 gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg                  <  4:431313/149‑1 (MQ=255)

                                                                                                                                  ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc                 >  3:703707/1‑149 (MQ=255)

                                                                                                                                   tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc                 >  2:569915/1‑148 (MQ=255)

                                                                                                                                    tgacattgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa               >  4:724210/8‑149 (MQ=255)

                                                                                                                                      tttCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCAt              >  1:164404/1‑148 (MQ=255)

                                                                                                                                      tttCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATa             >  3:120484/1‑149 (MQ=255)

                                                                                                                                        attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCtt                                                                             <  4:300782/80‑1 (MQ=255)

                                                                                                                                        attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCtt                                                                             >  3:300782/4‑83 (MQ=255)

                                                                                                                                        attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac           >  4:355174/4‑149 (MQ=255)

                                                                                                                                        attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac           >  3:451110/4‑149 (MQ=255)

                                                                                                                                          tgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTTCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaa                          <  1:633887/131‑1 (MQ=255)

                                                                                                                                              aCTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACct      >  2:534157/1‑148 (MQ=255)

GAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCA  >  NC_000913/858780‑859073

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