BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	pseudogene (915‑916/1358 nt)	gatC ←	galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,361	0	C	.	98.4%	1273.1 / 9.6	305	pseudogene (918/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand): ref base C (0/5); new base . (170/130); total (170/135)
Fisher's exact test for biased strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,362	0	C	.	98.7%	1270.3 / 8.6	305	pseudogene (917/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand): ref base C (0/4); new base . (170/130); total (170/135)
Fisher's exact test for biased strand distribution p-value = 3.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 364 total aligned reads displayed.

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCC  >  NC_000913/2173222‑2173494

||

tttCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATaccac                                                                                                                                       <  1:3418437/140‑1 (MQ=255)

 ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATaccacc                                                                                                                                      <  1:1450982/140‑1 (MQ=255)

    cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATaccaccgcc                                                                                                                                   >  1:4508019/1‑137 (MQ=255)

        aTGCACGGCGACCGCAATCGCCCCGCAGACGCCGATGATGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATaccaccgccgtat                                                                                                                               <  1:2071423/140‑8 (MQ=255)

         tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGACAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                                                                                            >  1:6899727/1‑140 (MQ=255)

          gCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                           >  1:3982190/1‑140 (MQ=255)

                 gACCGCCATCGCCACTAAGAAGTCGCTGGTGGCAAGAGCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                    <  1:3468522/140‑1 (MQ=255)

                  aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                   <  1:3781928/140‑1 (MQ=255)

                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                  >  2:6914168/1‑140 (MQ=255)

                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGAGGGATAAAGAT‑CTGGCTTGCCGATAGCA‑‑CCTCCGTATGTCCCagcagc                                                                                                                  >  2:2071423/1‑140 (MQ=255)

                        atagCCACGAAGAAGCCGATGGGGGTAATAACGCAAAA‑CGCGAGCACCTGATAACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGAATAAAAAT‑CAGGATTGCCGAAACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                             <  2:1071296/137‑1 (MQ=255)

                        aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                             <  1:5722633/140‑1 (MQ=255)

                        aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                             >  1:4329788/1‑140 (MQ=255)

                        aTCGCCACGAAGAAGCCCATGGTGGCAAGATCGCCTAA‑CCGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                             <  2:1994140/140‑1 (MQ=255)

                                gaagaaGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                     >  1:519826/1‑140 (MQ=255)

                                  agaagCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                  <  2:4548866/140‑1 (MQ=255)

                                   aaagCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                  <  1:5435714/139‑1 (MQ=255)

                                   aaagCCGATGGTGGAAAGATCGCAAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                  <  2:6102842/139‑1 (MQ=255)

                                     agCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                >  1:7005345/1‑140 (MQ=255)

                                      gCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                               >  2:4997598/1‑140 (MQ=255)

                                      gCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                               <  1:4591351/140‑1 (MQ=255)

                                       ccGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCGGCTGCAACGTCGGATCAAGGCCAAt                                                                                              >  1:2812689/1‑140 (MQ=255)

                                        cGATGGTGGCAAGATCGCCGAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCTGGATCAAGGCCAATc                                                                                             <  2:6537854/140‑1 (MQ=255)

                                          aTGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                             >  1:4891471/1‑122 (MQ=255)

                                           tggtggCACGATCGCCTGA‑CGGCAGCACCTCGTTCTCCGGCACACAAAAAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                          <  2:3625416/140‑1 (MQ=255)

                                           tggtggCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                          >  1:1194431/1‑140 (MQ=255)

                                           tggtggCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                          <  2:864193/140‑1 (MQ=255)

                                            ggtggCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                         >  2:5078367/1‑140 (MQ=255)

                                            ggtggCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAATGCCGCATCAAGGCCAATCAGGa                                                                                         >  2:2942279/1‑140 (MQ=255)

                                             gtggCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                        <  1:5078367/140‑1 (MQ=255)

                                                     aTCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCTTATGTCCCAGCAGCAAAGCCGGATCAAGGCGAATCAGGAACTCCTGAc                                                                                >  2:2328064/1‑140 (MQ=255)

                                                     aTCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGAAGCAACGACGGCTCAAGGCCAATCAGGAACTCCTGAc                                                                                >  1:2408378/1‑140 (MQ=255)

                                                       cGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCATCAACGCCCGATTAACGCCAATAAGGAACTCCTGAccg                                                                              >  1:3082546/1‑140 (MQ=255)

                                                       cGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                              >  2:519618/1‑140 (MQ=255)

                                                        gCCAAA‑CGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                             <  1:7368998/140‑1 (MQ=255)

                                                        gCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTACGTCCCAGCCGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                             >  2:2592667/1‑140 (MQ=255)

                                                        cccAAG‑CGGCACCACCTGATCCCCCGGCAAACAAACAGGAATTAAATGGGTGCGTGGGAAAAAAAC‑CAGGCTCGCCTACACAA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                             <  2:1879545/139‑1 (MQ=255)

                                                         ccAAACCGGCAGCACCTTATTCCCCGCCACACAAACAGCAATTAAAATGGTGAGTGGCAAAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                             <  2:1350410/140‑1 (MQ=255)

                                                               cGGCAGCACCTGATTCCGCGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                       <  2:3390175/140‑1 (MQ=255)

                                                                 gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                     >  2:6184783/1‑140 (MQ=255)

                                                                   agcaCCTGATTCCCCCGCACACCAACAGCTATTAAAACGGTGAGTGGGATAAAAAC‑CAGGCTTGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                   <  1:1916365/140‑1 (MQ=255)

                                                                    gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                  <  1:5251342/140‑1 (MQ=255)

                                                                           aTTCCCCGGCACACAAACAGCAATTAAAATGGTGAGAGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                           <  2:6633443/140‑1 (MQ=255)

                                                                                ccGGCACACAAACAGAAATTAAAATGGTGAGTGGGAAAAAAAT‑CAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                      <  2:7143817/140‑1 (MQ=255)

                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                   >  2:223215/1‑125 (MQ=255)

                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                   <  2:6807763/123‑1 (MQ=255)

                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                  >  1:6797495/1‑140 (MQ=255)

                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                  >  1:784033/1‑140 (MQ=255)

                                                                                     acacaAACATCAATTAAAAGGGTGAGTGGGAAAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:4395197/140‑1 (MQ=255)

                                                                                     acacaAACAGGAATAAAAAGGGTGAGTGGGATAAAAAT‑AAGTCTTGCGGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:6980102/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTCAAATGGGGAATGGTATAAAAAC‑CAGGCTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  1:2942279/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTCAAATGGCGACTGGTATAAAAAA‑CAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  1:2231788/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTCAAATGGAGACGGGTATAAAAAA‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  1:2133977/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:7571348/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:5873390/140‑1 (MQ=255)

                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 >  1:4783157/1‑140 (MQ=255)

                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCACCGTACTTCGCCTGTAAACGACTACGAGc                                                 >  1:3170749/1‑140 (MQ=255)

                                                                                     acacaAACAGAAATTAAAATGGTGAGGGGAAAAAAAAT‑CAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:1239253/140‑1 (MQ=255)

                                                                                     acacaAACAGAAATTAAAATGGCGAGTGGGAAAAACAA‑CCGGCTTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                 <  2:2812689/140‑1 (MQ=255)

                                                                                            cAGAAATAAAAAGGGTGAGGGGGATAAAAAT‑CAGGCTGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                          <  1:1821852/140‑1 (MQ=255)

                                                                                                 ataaaaaTGGTGAGGGGGATAAAAAT‑CAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                     <  2:6698752/137‑1 (MQ=255)

                                                                                                 aTTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                     <  1:3890758/140‑1 (MQ=255)

                                                                                                 aTTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                     >  2:5418917/1‑140 (MQ=255)

                                                                                                 aTTAAAATGGTGAGTGGGAAAAAAAT‑CAGGTATGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                     <  2:3082546/140‑1 (MQ=255)

                                                                                                   tAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGactg                                   >  1:5612337/1‑137 (MQ=255)

                                                                                                     aaaTGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                 >  2:5643668/1‑140 (MQ=255)

                                                                                                       aTGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                               >  1:982205/1‑140 (MQ=255)

                                                                                                          gTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                            >  2:5959360/1‑140 (MQ=255)

                                                                                                           tGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                           <  2:163386/140‑1 (MQ=255)

                                                                                                           tGAGTGGGATAAAAAT‑CAGGCTGGCCGATCCCA‑‑CCGCGGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                           <  1:2592667/140‑1 (MQ=255)

                                                                                                            gagaGGGATAAAAAT‑GGGGCTGTCCGAT‑CC‑GCCCGCCGTCTGTCCCAGTAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                          <  2:2137759/140‑1 (MQ=255)

                                                                                                            gAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                >  1:4209949/1‑134 (MQ=255)

                                                                                                            gAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                          >  2:4878169/1‑140 (MQ=255)

                                                                                                            gAGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAACGACAACGAGCATGCTTAACGATGGGCGTTAAAc                          >  1:3648416/1‑140 (MQ=255)

                                                                                                             aGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                         >  1:1055664/1‑140 (MQ=255)

                                                                                                             aGTGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                         >  1:1679163/1‑140 (MQ=255)

                                                                                                               tGGGATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                       >  2:3259001/1‑140 (MQ=255)

                                                                                                                gggATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                      <  1:6004558/140‑1 (MQ=255)

                                                                                                                gggATAAAAAG‑AAGGCTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                      <  2:3858313/140‑1 (MQ=255)

                                                                                                                gggAAAAAAATACAGGCTTGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                       <  2:7379231/140‑1 (MQ=255)

                                                                                                                  gATAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                    >  1:1300880/1‑140 (MQ=255)

                                                                                                                    tAAAAAT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                  >  1:4437886/1‑140 (MQ=255)

                                                                                                                    aaaaaaT‑CAGGCTTGCAGATCCGA‑‑CCGCCGTATTTTCCACCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                  <  2:3448576/139‑1 (MQ=255)

                                                                                                                     aaaaaT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                 <  1:4101057/140‑1 (MQ=255)

                                                                                                                     aaaaaT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                 <  1:5905432/140‑1 (MQ=255)

                                                                                                                     aaaaaT‑CAGACTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                 <  1:6103562/140‑1 (MQ=255)

                                                                                                                      aaaaT‑CAGGCTTGCCCATAC‑‑CCCCGACGTTTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                <  1:2299682/140‑1 (MQ=255)

                                                                                                                       aaaT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg               >  2:1536004/1‑140 (MQ=255)

                                                                                                                       aaaT‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg               >  1:1520249/1‑140 (MQ=255)

                                                                                                                         aT‑CAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg             <  2:4123978/140‑1 (MQ=255)

                                                                                                                          t‑CAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCATGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt            <  2:423275/140‑1 (MQ=255)

                                                                                                                             aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGATTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt          >  1:7425877/1‑140 (MQ=255)

                                                                                                                             aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt          >  2:5535120/1‑140 (MQ=255)

                                                                                                                             aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGAGCGTTAAACCATCCATGATGGGttt          >  1:6540799/1‑140 (MQ=255)

                                                                                                                                cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt       <  1:921148/140‑1 (MQ=255)

                                                                                                                                cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt       <  2:3997086/140‑1 (MQ=255)

                                                                                                                                  ggCCGATAC‑‑CCCCGCCGCATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa     <  2:2448118/139‑1 (MQ=255)

                                                                                                                                   gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc    >  1:6531985/1‑140 (MQ=255)

                                                                                                                                   gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc    >  2:4706852/1‑140 (MQ=255)

                                                                                                                                    ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc   >  1:7291871/1‑140 (MQ=255)

||

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAA‑CGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCC  >  NC_000913/2173222‑2173494

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