BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,893,400	G→A	9.2%	L248L (CTG→CTA)	ygcN →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,893,400	0	G	A	9.2%	321.6 / 16.1	106	L248L (CTG→CTA)	ygcN	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (49/47); new base A (7/3); total (56/50)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGA  >  NC_000913/2893253‑2893546

cGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAgctgc                                                                                                                                                   <  2:676513/149‑1 (MQ=255)

 gAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAgctgct                                                                                                                                                  >  1:979217/1‑149 (MQ=255)

 gAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAgctgct                                                                                                                                                  >  1:1066611/1‑149 (MQ=255)

   aGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTGGTGTTCCGGCAAGCGAGCTGCTGa                                                                                                                                                <  2:189434/149‑1 (MQ=255)

    ggCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGAc                                                                                                                                               >  2:181144/1‑149 (MQ=255)

        cAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCg                                                                                                                                            >  3:1010576/1‑148 (MQ=255)

         aGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGc                                                                                                                                           <  4:1016314/148‑1 (MQ=255)

         aGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGc                                                                                                                                           <  3:733046/148‑1 (MQ=255)

              tGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaa                                                                                                                                      <  3:174793/148‑1 (MQ=255)

              tGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaa                                                                                                                                      <  4:365319/148‑1 (MQ=255)

              tGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaa                                                                                                                                      <  4:1010576/148‑1 (MQ=255)

                cTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAg                                                                                                                                    <  3:436085/148‑1 (MQ=255)

                    tCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCa                                                                                                                                >  4:1019748/1‑148 (MQ=255)

                       gCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCg                                                                                                                            >  1:258548/1‑149 (MQ=255)

                         ggCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGc                                                                                                                          >  2:1067644/1‑149 (MQ=255)

                          gCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAg                                                                                                                                    >  3:679567/1‑138 (MQ=255)

                          gCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAg                                                                                                                                    <  4:679567/138‑1 (MQ=255)

                               aTCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTgc                                                                                                                     >  3:1043068/1‑148 (MQ=255)

                                   gtgATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGccc                                                                                                                 >  4:394396/1‑148 (MQ=255)

                                    tgatgaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGc                                                                                                               <  2:63370/149‑1 (MQ=255)

                                      atgaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCtt                                                                                                             >  1:483569/1‑149 (MQ=255)

                                       tgaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaaa                                                                                                                                     >  2:322903/1‑124 (MQ=255)

                                       tgaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTaaa                                                                                                                                     <  1:322903/124‑1 (MQ=255)

                                        gaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTAt                                                                                                           <  2:1066611/149‑1 (MQ=255)

                                        gaCTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTAt                                                                                                           <  2:979217/149‑1 (MQ=255)

                                             aCCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCGCGCTCTCTTCCCTTACGCAAAATCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTa                                                                                                            <  1:990936/143‑1 (MQ=255)

                                             aCCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCaaaa                                                                                                      >  2:217042/1‑149 (MQ=255)

                                             aCCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGAGTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTa                                                                                                            >  2:990936/1‑143 (MQ=255)

                                              ccgggcggcGCATTTCTTTATACTAATCAACCATCGCTCTCGTTAGGGATTGTTTGCCCGGTCTCTTCGCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCaaaa                                                                                                      <  4:553720/145‑1 (MQ=255)

                                              cccGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCaaaaa                                                                                                     <  1:97725/149‑1 (MQ=255)

                                                cggcggCGCATTTCTTTATACTAGTCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACTCAAAGTCGTGTTCCGGCAAGCGAGCTGATGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAac                                                                                                    >  3:663868/1‑148 (MQ=255)

                                                cggcggCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGc                                                                                                               <  1:368339/137‑1 (MQ=255)

                                                cggcggCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGc                                                                                                               >  2:368339/1‑137 (MQ=255)

                                                cggcggCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAac                                                                                                    >  4:499220/1‑148 (MQ=255)

                                                cggcggCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAac                                                                                                    <  3:223642/148‑1 (MQ=255)

                                                     gcgcATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGaa                                                                                              <  2:815520/149‑1 (MQ=255)

                                                          tttCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCAc                                                                                          <  3:7889/148‑1 (MQ=255)

                                                          tttCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACt                                                                                         <  1:844656/149‑1 (MQ=255)

                                                           ttCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACt                                                                                         <  3:500235/148‑1 (MQ=255)

                                                             cTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGa                                                                                      <  1:181144/149‑1 (MQ=255)

                                                              tttATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAg                                                                                     >  1:577085/1‑149 (MQ=255)

                                                              tttATACTAATAAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAg                                                                                     >  1:564029/1‑149 (MQ=255)

                                                                    cTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGt                                                                               <  2:813145/149‑1 (MQ=255)

                                                                    cTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGt                                                                               <  2:500104/149‑1 (MQ=255)

                                                                     tAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGt                                                                               <  3:12036/148‑1 (MQ=255)

                                                                        tCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTgcgc                                                                           >  1:882203/1‑149 (MQ=255)

                                                                        tCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTgcg                                                                            >  3:529381/1‑148 (MQ=255)

                                                                         caacaaACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTgcgc                                                                           >  3:569337/1‑148 (MQ=255)

                                                                            caaACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAGGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATc                                                                        >  4:162363/1‑148 (MQ=255)

                                                                              aaCGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTg                                                                      >  3:846726/1‑148 (MQ=255)

                                                                               aCGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTgg                                                                     <  4:425092/148‑1 (MQ=255)

                                                                               aCGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTgg                                                                     >  3:445325/1‑148 (MQ=255)

                                                                                  ctctCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGcc                                                                 <  2:1038718/149‑1 (MQ=255)

                                                                                   tctcGTTAGGGATTGTTTGCCTGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCa                                                                >  1:34868/1‑149 (MQ=255)

                                                                                   tctcGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCa                                                                >  2:135261/1‑149 (MQ=255)

                                                                                       gTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAg                                                            >  1:424336/1‑149 (MQ=255)

                                                                                        ttAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATc                                                                        <  2:142758/136‑1 (MQ=255)

                                                                                        ttAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATc                                                                        >  1:142758/1‑136 (MQ=255)

                                                                                         tAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAgg                                                           <  4:446944/148‑1 (MQ=255)

                                                                                          aGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGt                                                          <  4:569337/148‑1 (MQ=255)

                                                                                           gggATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTg                                                         >  4:5342/1‑148 (MQ=255)

                                                                                            ggATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTgg                                                        >  4:780188/1‑148 (MQ=255)

                                                                                            ggATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTgg                                                        <  3:393370/148‑1 (MQ=255)

                                                                                              aTTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCtt                                                     >  2:995585/1‑149 (MQ=255)

                                                                                                 gTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGca                                                  <  1:995585/149‑1 (MQ=255)

                                                                                                  tttGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGcac                                                 >  1:910005/1‑149 (MQ=255)

                                                                                                   ttGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAAt                                                                                             >  4:331581/1‑104 (MQ=255)

                                                                                                   ttGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGcc                                                                 <  3:1045214/132‑1 (MQ=255)

                                                                                                    tGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAg                                               <  2:13649/149‑1 (MQ=255)

                                                                                                     gCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATg                                                                                 >  2:691149/1‑114 (MQ=255)

                                                                                                     gCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATg                                                                                 <  1:691149/114‑1 (MQ=255)

                                                                                                      cccGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGt                                              <  3:886669/148‑1 (MQ=255)

                                                                                                       ccGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTa                                             <  3:74846/148‑1 (MQ=255)

                                                                                                         gCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATg                                           <  4:1043068/148‑1 (MQ=255)

                                                                                                          ctctctTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGTTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGcc                                         >  2:435698/1‑149 (MQ=255)

                                                                                                           tctctTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCg                                        >  2:453440/1‑149 (MQ=255)

                                                                                                            ctctTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCg                                        >  4:43821/1‑148 (MQ=255)

                                                                                                             tctTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGt                                      <  1:147515/149‑1 (MQ=255)

                                                                                                                tCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGc                                    >  3:529666/1‑148 (MQ=255)

                                                                                                                  ccTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAAt                                 >  1:620327/1‑149 (MQ=255)

                                                                                                                  ccTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAAt                                 <  2:10858/149‑1 (MQ=255)

                                                                                                                      aCGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGc                             >  2:276885/1‑149 (MQ=255)

                                                                                                                          aaaGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGt                         >  2:19179/1‑149 (MQ=255)

                                                                                                                          aaaGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGt                         >  2:831433/1‑149 (MQ=255)

                                                                                                                            aGTCGTGTTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTa                        >  3:714654/1‑148 (MQ=255)

                                                                                                                             gTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAAc                      >  2:120735/1‑149 (MQ=255)

                                                                                                                              tCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAAc                      <  4:846726/148‑1 (MQ=255)

                                                                                                                                gtgtTCCGGCAAGCGAGCTACTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACggc                   >  2:489456/1‑149 (MQ=255)

                                                                                                                                 tgtTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGTCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACggc                   >  4:1008053/1‑148 (MQ=255)

                                                                                                                                       ggCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCGATACGCCGGTAACGgctggctg             >  4:742160/1‑148 (MQ=255)

                                                                                                                                       ggCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGgctgc            >  1:947108/1‑149 (MQ=255)

                                                                                                                                        gCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTg                                                                                        >  2:432516/1‑72 (MQ=255)

                                                                                                                                        gCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTg                                                                                        <  1:432516/72‑1 (MQ=255)

                                                                                                                                             cGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCTCTTATCAAAATCACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGGTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCtggtgt      <  2:1052098/149‑2 (MQ=255)

                                                                                                                                             cGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCtggtg       >  4:1044476/1‑148 (MQ=255)

                                                                                                                                             cGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCtggtg       >  3:877428/1‑148 (MQ=255)

                                                                                                                                               aGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGc    <  1:276885/149‑1 (MQ=255)

                                                                                                                                                gctgctGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCg   <  1:217042/149‑1 (MQ=255)

                                                                                                                                                gctgctGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCg   <  1:135261/149‑1 (MQ=255)

                                                                                                                                                  tgctgACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGa  <  3:578053/148‑1 (MQ=255)

CGAAGGCGCAGCGTTGCTGTTCAGCGGCAGGATCTGTGATGACTTACCCGGCGGCGCATTTCTTTATACTAATCAACAAACGCTCTCGTTAGGGATTGTTTGCCCGCTCTCTTCCCTTACGCAAAGTCGTGTTCCGGCAAGCGAGCTGCTGACTCGCTTTAAAGCGCATCCGGCAGTGCGCCCGCTTATCAAAAACACGGAATCACTGGAGTATGGTGCGCATCTGGTGCCAGAAGGTGGCTTGCACAGTATGCCGGTGCAATACGCCGGTAACGGCTGGCTGCTGGTGGGCGA  >  NC_000913/2893253‑2893546

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