BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	266,465	G→A	5.5%	A30V (GCC→GTC)	yafY ←	lipoprotein, inner membrane; degP regulator; CP4‑6 prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	266,465	0	G	A	5.5%	365.5 / 7.1	108	A30V (GCC→GTC)	yafY	lipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand): ref base G (50/52); new base A (3/3); total (53/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 111 total aligned reads displayed.

GCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAAT  >  NC_000913/266320‑266612

gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGcc                                                                                                                                                   >  2:182993/1‑148 (MQ=255)

 ccATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCt                                                                                                                                                  <  2:482090/148‑1 (MQ=255)

    tCGGGCTACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCtgaa                                                                                                                                               >  2:619035/1‑145 (MQ=255)

         cGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtca                                                                                                                                          <  4:604/148‑1 (MQ=255)

                cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                                                                                                                                >  3:428154/1‑135 (MQ=255)

                cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                                                                                                                                <  4:428154/135‑1 (MQ=255)

                cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGCAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagtaa                                                                                                                                   <  1:619035/148‑5 (MQ=18)

                 gggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                                                                                                                                  <  1:182993/148‑1 (MQ=255)

                  ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                                                                                                                                 <  1:783508/148‑1 (MQ=255)

                        tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCAc                                                                                                                                      <  1:765613/137‑1 (MQ=255)

                        tAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCg                                                                                                                           >  3:576670/1‑148 (MQ=255)

                         aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                                                                                                                                 >  2:473006/1‑141 (MQ=255)

                         aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                                                                                                                                 <  1:473006/141‑1 (MQ=255)

                         aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccat                                                                                                                                         <  1:444799/133‑1 (MQ=255)

                         aaaaACTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccat                                                                                                                                         >  2:444799/1‑133 (MQ=18)

                            aCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg                                                                                                                       <  3:803253/148‑1 (MQ=255)

                            aCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATg                                                                                                                       >  4:803253/1‑148 (MQ=255)

                              cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCg                                                                                                                     <  4:576670/148‑1 (MQ=255)

                                    ggTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgaga                                                                                                               >  3:417091/1‑148 (MQ=255)

                                     gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa                                                                                                              <  1:508617/148‑1 (MQ=255)

                                         gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          <  1:581552/148‑1 (MQ=255)

                                         gACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          >  2:581552/1‑148 (MQ=255)

                                           ccGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACagag                                                                                                        >  2:655309/1‑148 (MQ=255)

                                               aTAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggt                                                                                                    >  1:858910/1‑148 (MQ=255)

                                                tAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtggtg                                                                                                   >  1:754315/1‑148 (MQ=255)

                                                 aaGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGtggtgg                                                                                                  <  4:417091/148‑1 (MQ=255)

                                                      gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCg                                                                                                                     >  2:793757/1‑124 (MQ=255)

                                                      gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCg                                                                                                                     <  1:793757/124‑1 (MQ=255)

                                                      gtCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAAcc                                                                                             >  3:950877/1‑148 (MQ=255)

                                                               tgACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                                                                                                                                  <  2:880938/102‑1 (MQ=255)

                                                                 aCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagcag                                                                                  <  1:655309/148‑1 (MQ=255)

                                                                 aCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCAGCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCcagcag                                                                                  <  2:858910/148‑1 (MQ=255)

                                                                  ccATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                                                                                   >  4:1056411/1‑98 (MQ=255)

                                                                  ccATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                                                                                   <  3:1056411/98‑1 (MQ=255)

                                                                  ccATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagtaa                                                                                                                                   >  3:324474/1‑94 (MQ=18)

                                                                  ccATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGACCTGAAGGccatcagtaa                                                                                                                                   <  4:324474/98‑5 (MQ=18)

                                                                    aTTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa                                                                                                              <  2:663704/117‑1 (MQ=255)

                                                                      tcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGa                                                                                                                    >  1:1117565/1‑109 (MQ=255)

                                                                      tcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGa                                                                                                                    <  2:1117565/109‑1 (MQ=255)

                                                                      tcttcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCa                                                                             >  4:299990/1‑148 (MQ=255)

                                                                         tcGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAgtg                                                                          <  3:299990/148‑1 (MQ=255)

                                                                            cTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg                                                                                                                  >  4:467824/1‑105 (MQ=255)

                                                                                  tcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                                                                                   >  1:15835/1‑82 (MQ=255)

                                                                                  tcatcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc                                                                                                                                   <  2:15835/82‑1 (MQ=255)

                                                                                    atcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          <  3:756176/105‑1 (MQ=255)

                                                                                    atcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          >  4:756176/1‑105 (MQ=255)

                                                                                     tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGtt                                                                 <  2:242053/145‑1 (MQ=255)

                                                                                     tcCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGtt                                                                 >  1:242053/1‑145 (MQ=255)

                                                                                      ccTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          >  4:690313/1‑103 (MQ=255)

                                                                                      ccTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACag                                                                                                          <  3:690313/103‑1 (MQ=255)

                                                                                      ccTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCat                                                             >  1:260026/1‑148 (MQ=255)

                                                                                       cTGGTGGCTGTCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                                                                                        <  3:335763/72‑1 (MQ=255)

                                                                                       cTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                                                                                                                                        >  4:335763/1‑72 (MQ=255)

                                                                                        tggtggCTGACCACACCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCatta                                                           >  1:336601/1‑148 (MQ=255)

                                                                                         ggtggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAgtgg                                                                                                     >  1:222845/1‑105 (MQ=255)

                                                                                         ggtggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCattat                                                          <  2:260026/148‑1 (MQ=255)

                                                                                           tggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGATCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTAttt                                                        <  4:950877/148‑1 (MQ=255)

                                                                                             gCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCaga                                                                                 >  4:914316/1‑121 (MQ=255)

                                                                                             gCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCaga                                                                                 <  3:914316/121‑1 (MQ=255)

                                                                                                 aCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                                                                                                                                 <  2:187857/53‑1 (MQ=255)

                                                                                                 aCCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                                                                                                                                                 >  1:187857/1‑53 (MQ=255)

                                                                                                   cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCa                                                                                                <  2:941477/100‑1 (MQ=255)

                                                                                                   cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCa                                                                                                >  1:941477/1‑100 (MQ=255)

                                                                                                   cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg                                                                            <  1:834350/120‑1 (MQ=255)

                                                                                                   cacaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAg                                                                            >  2:834350/1‑120 (MQ=255)

                                                                                                    acaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt                                                                                                                                     <  2:705933/62‑1 (MQ=255)

                                                                                                    acaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt                                                                                                                                     >  1:705933/1‑62 (MQ=255)

                                                                                                    acaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCaa                                                                                               >  4:612079/1‑100 (MQ=255)

                                                                                                    acaTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCaa                                                                                               <  3:612079/100‑1 (MQ=255)

                                                                                                      accACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGaaa                                             <  2:336601/146‑1 (MQ=255)

                                                                                                      aTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg                                                                                                                  >  4:60473/1‑79 (MQ=255)

                                                                                                      aTCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg                                                                                                                  <  3:60473/79‑1 (MQ=255)

                                                                                                          cTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAg                                         >  4:219696/1‑148 (MQ=255)

                                                                                                                           gAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCaa                        <  3:219696/148‑1 (MQ=255)

                                                                                                                               tCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACg                                                                   <  2:88540/101‑1 (MQ=255)

                                                                                                                               tCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACg                                                                   >  1:88540/1‑101 (MQ=255)

                                                                                                                                ccTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgag                                                                                                                <  2:874798/55‑1 (MQ=255)

                                                                                                                                ccTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgag                                                                                                                >  1:874798/1‑55 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCCCTTCTGTCATCGCATGCGATGAGGAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACa                   <  1:322688/145‑1 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                                                                                                                         <  1:613634/43‑1 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCa                                                                                                                         >  2:613634/1‑43 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCagag                                                                                <  1:444517/84‑1 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCagag                                                                                >  2:444517/1‑84 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGttt                                                                       >  2:339081/1‑93 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGttt                                                                       <  1:339081/93‑1 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAg                                         <  3:473839/123‑1 (MQ=255)

                                                                                                                                   tAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACa                   >  2:322688/1‑145 (MQ=255)

                                                                                                                                      aTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCATAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGGCATAAATg             <  2:921368/148‑1 (MQ=255)

                                                                                                                                      aTTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATg             >  1:921368/1‑148 (MQ=255)

                                                                                                                                        ttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCagag                                                                                <  1:662349/79‑1 (MQ=255)

                                                                                                                                        ttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCagag                                                                                >  2:662349/1‑79 (MQ=255)

                                                                                                                                        ttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTc                                                   >  1:683690/1‑108 (MQ=255)

                                                                                                                                        ttACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTc                                                   <  2:683690/108‑1 (MQ=255)

                                                                                                                                            tGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                           <  4:328671/112‑1 (MQ=255)

                                                                                                                                            tGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                           >  3:328671/1‑112 (MQ=255)

                                                                                                                                             gATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCa                                                              >  1:998042/1‑92 (MQ=255)

                                                                                                                                             gATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCa                                                              <  2:998042/92‑1 (MQ=255)

                                                                                                                                             gATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGCAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGa      >  2:83806/1‑148 (MQ=255)

                                                                                                                                                 gCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACttt                             <  1:968599/121‑1 (MQ=255)

                                                                                                                                                 gCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAAt  <  2:754315/148‑1 (MQ=255)

GCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAAT  >  NC_000913/266320‑266612

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