BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,676,370	C→G	5.2%	A661P (GCG→CCG)	yhjG ←	AsmA family protein YhjG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,676,370	0	C	G	5.2%	512.8 / 13.8	153	A661P (GCG→CCG)	yhjG	AsmA family protein YhjG
Reads supporting (aligned to +/- strand): ref base C (72/73); new base G (2/6); total (74/79)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

Only 100 of 165 total aligned reads displayed.

CAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTCCG  >  NC_000913/3676223‑3676519

  aaGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACgccgccgcc                                                                                                                                                   >  1:1413828/1‑150 (MQ=255)

     cacaCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgt                                                                                                                                                >  1:1821166/1‑150 (MQ=255)

      acacTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgtg                                                                                                                                               <  2:907592/150‑1 (MQ=255)

      acacTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgtg                                                                                                                                               <  1:1464698/150‑1 (MQ=255)

      acacTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgtg                                                                                                                                               <  1:597200/150‑1 (MQ=255)

       cacTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgtgt                                                                                                                                              >  2:514209/1‑150 (MQ=255)

         cTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTa                                                                                                                                            <  1:758061/150‑1 (MQ=255)

           ttGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACACCGGCGGCGGCGTTAct                                                                                                                                          <  2:195915/150‑2 (MQ=255)

             gCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACACCGGCGGCGGCGTTActgg                                                                                                                                        <  1:336760/150‑4 (MQ=255)

              ccAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGc                                                                                                                                       <  1:1053178/150‑1 (MQ=255)

               caacaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCg                                                                                                                                      >  1:1647677/1‑150 (MQ=255)

                aacaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGt                                                                                                                                     <  1:445524/150‑1 (MQ=255)

                  caaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTg                                                                                                                                   <  1:1573035/150‑1 (MQ=255)

                   aaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACACCAGCGGCGGCGTTActgg                                                                                                                                        <  2:2014393/144‑4 (MQ=255)

                    aTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGcc                                                                                                                                 <  2:1728742/150‑1 (MQ=255)

                        aaaTCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTg                                                                                                                             <  1:2106483/150‑1 (MQ=255)

                        aaaTCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTg                                                                                                                             >  2:954330/1‑150 (MQ=255)

                        aaaTCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTAGCCAGTg                                                                                                                             >  2:2233444/1‑150 (MQ=255)

                          aTCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGcc                                                                                                                           <  1:1904757/150‑1 (MQ=255)

                           tCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCg                                                                                                                          <  1:2446872/150‑1 (MQ=255)

                            cTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgc                                                                                                                         <  1:789171/150‑1 (MQ=255)

                              gTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgcag                                                                                                                       >  1:2414322/1‑150 (MQ=255)

                               tttAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgcagc                                                                                                                      <  1:886998/150‑1 (MQ=255)

                                  aaTTGCGGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAAc                                                                                                                   <  2:1640317/150‑1 (MQ=255)

                                  aaTTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAAc                                                                                                                   >  1:2086938/1‑150 (MQ=255)

                                  aaTTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAAc                                                                                                                   >  2:910115/1‑150 (MQ=255)

                                   aTTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACa                                                                                                                  <  1:2413726/150‑1 (MQ=255)

                                     tGCCGTATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGc                                                                                                                <  1:46197/150‑1 (MQ=255)

                                     tGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGc                                                                                                                <  1:2385692/150‑1 (MQ=255)

                                      gCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGca                                                                                                               <  1:137507/150‑1 (MQ=255)

                                      gCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGca                                                                                                               <  1:209181/150‑1 (MQ=255)

                                      gCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGca                                                                                                               >  1:1627604/1‑150 (MQ=255)

                                      gCCGGATGTCAGACATCCGGAAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGca                                                                                                               <  1:2288327/150‑1 (MQ=255)

                                       ccGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCTGTGTTACCAGCGTTGCCAGTGCCGAAGCAACAGcac                                                                                                              <  2:1856197/150‑1 (MQ=255)

                                       ccGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcac                                                                                                              >  2:304722/1‑150 (MQ=255)

                                        cGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcacc                                                                                                             >  2:2105549/1‑150 (MQ=255)

                                         ggATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcacca                                                                                                            >  2:1896323/1‑150 (MQ=255)

                                          gATGTCTGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcaccac                                                                                                           >  1:1839464/1‑150 (MQ=255)

                                             gTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACggg                                                                                                        <  1:1839995/150‑1 (MQ=255)

                                              tctgACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGc                                                                                                       >  1:2001018/4‑150 (MQ=255)

                                              tCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGc                                                                                                       <  1:970422/150‑1 (MQ=255)

                                               cAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCa                                                                                                      >  2:2317363/1‑150 (MQ=255)

                                                    aTCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACACCGGCGGCGGCGTTAct                                                                                                                                          >  2:732439/1‑108 (MQ=255)

                                                     tCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCACCAGCACCACGGGCAATCAGc                                                                                                >  2:359953/1‑150 (MQ=255)

                                                     tCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGc                                                                                                <  1:1093859/150‑1 (MQ=255)

                                                     tCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGc                                                                                                >  1:1039950/1‑150 (MQ=255)

                                                      ccGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCAGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCg                                                                                               >  1:1553044/1‑150 (MQ=255)

                                                       cGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCgg                                                                                              >  1:2490449/1‑150 (MQ=255)

                                                        ggCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGc                                                                                                       >  1:1711158/1‑140 (MQ=255)

                                                         gCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTc                                                                                            <  2:633636/150‑1 (MQ=255)

                                                          cAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTcc                                                                                           >  1:1828811/1‑150 (MQ=255)

                                                            aTTATCATCACTTCTTCATCTGCGCCAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCgg                                                                                         >  2:1334299/1‑150 (MQ=255)

                                                              tATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCt                                                                                       <  1:1158192/150‑1 (MQ=255)

                                                              tATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCt                                                                                       <  2:2170003/150‑1 (MQ=255)

                                                                 catcaCTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTc                                                                                    <  2:1471048/150‑1 (MQ=255)

                                                                  atcaCTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCGGCGGTCCGGCTTTCa                                                                                   <  1:715095/150‑1 (MQ=255)

                                                                    caCTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCAcc                                                                                 >  2:565396/1‑150 (MQ=255)

                                                                    caCTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCAcc                                                                                 >  1:2148390/1‑150 (MQ=255)

                                                                      cttcttCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGGAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCAcccc                                                                               <  2:1108542/150‑1 (MQ=255)

                                                                        tcttcATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAg                                                                             >  1:682420/1‑150 (MQ=255)

                                                                         cttcATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGc                                                                            <  2:400008/150‑1 (MQ=255)

                                                                          ttcATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGcc                                                                           <  1:802146/150‑1 (MQ=255)

                                                                             aTCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGc                                                                        <  2:1214940/150‑1 (MQ=255)

                                                                                tGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGa                                                                     >  1:1271874/1‑150 (MQ=255)

                                                                                tGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTGTCACCCCAGCCTGCGGa                                                                     >  2:1871590/1‑150 (MQ=255)

                                                                                  cGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGAtt                                                                   >  1:1759122/1‑150 (MQ=255)

                                                                                      aaaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCTGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTCCACCCCAGCCTGCGGATTTTTa                                                               <  2:1553044/150‑1 (MQ=255)

                                                                                      aaaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTa                                                               <  1:395899/150‑1 (MQ=255)

                                                                                      aaaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTa                                                               <  2:1413828/150‑1 (MQ=255)

                                                                                       aaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTaa                                                              <  1:910115/150‑1 (MQ=255)

                                                                                         aTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAAc                                                            <  2:2396368/150‑1 (MQ=255)

                                                                                            gTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTc                                                         <  2:1828811/150‑1 (MQ=255)

                                                                                             tCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGGCGGTGTTACCAGGGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTcc                                                        <  2:2195252/150‑1 (MQ=255)

                                                                                                 acacTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCg                                                    >  2:469922/1‑150 (MQ=255)

                                                                                                     tGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGcaca                                                <  2:2414322/150‑1 (MQ=255)

                                                                                                      gATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACAt                                               >  2:1330357/1‑150 (MQ=255)

                                                                                                         tAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCCGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGGTTTCACCCCCGCCTGCGGATTTTTAGGCGTCCCCCGcccctccg                                            >  2:1569402/1‑143 (MQ=255)

                                                                                                           gCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCCCCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAgc                                          >  2:299415/1‑150 (MQ=255)

                                                                                                              tCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAgcggc                                       >  1:234146/1‑150 (MQ=255)

                                                                                                                cccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGa                                     >  1:401094/1‑150 (MQ=255)

                                                                                                                cccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCGTTCCCCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGa                                     >  2:1293451/1‑150 (MQ=255)

                                                                                                                 ccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATAAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGaa                                    >  1:19663/1‑150 (MQ=255)

                                                                                                                  cTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAAc                                   <  1:1334299/150‑1 (MQ=255)

                                                                                                                  cTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAAc                                   <  2:732913/150‑1 (MQ=255)

                                                                                                                  cTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAAc                                   <  1:268929/150‑1 (MQ=255)

                                                                                                                       gAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAgt                              <  2:576389/150‑1 (MQ=255)

                                                                                                                       gAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAgt                              >  1:1100980/1‑150 (MQ=255)

                                                                                                                          ggTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGtgtg                           >  1:1730197/1‑150 (MQ=255)

                                                                                                                          ggTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGtgtg                           <  2:1410246/150‑1 (MQ=255)

                                                                                                                                  cAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCg                   >  1:1213817/1‑150 (MQ=255)

                                                                                                                                  cAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCg                   <  1:1232812/150‑1 (MQ=255)

                                                                                                                                     tGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAAt                <  2:1839464/150‑1 (MQ=255)

                                                                                                                                     tGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGA‑TTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAAtt               >  2:27587/1‑150 (MQ=255)

                                                                                                                                      gCCATTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAAtt               <  2:1421187/150‑1 (MQ=255)

                                                                                                                                          gTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCtt           <  2:2001018/150‑1 (MQ=255)

                                                                                                                                               gccgccgccGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTActc      >  2:1636861/1‑150 (MQ=255)

                                                                                                                                                  gccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTcc   <  2:78480/150‑1 (MQ=255)

                                                                                                                                                  gccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTcc   <  2:1039950/150‑1 (MQ=255)

                                                                                                                                                   ccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTCCg  >  1:1456030/1‑150 (MQ=255)

                                                                                                                                                   ccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTCCg  >  2:1955092/1‑150 (MQ=255)

CAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTCCG  >  NC_000913/3676223‑3676519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`