BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,734,543	(C)_8→7	100%	coding (492/732 nt)	yfiH ←	UPF0124 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,734,536	0	C	.	81.8%	282.7 / 72.1	99	coding (499/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand): ref base C (11/6); new base . (42/39); total (54/45)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 106 total aligned reads displayed.

GCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCC  >  NC_000913/2734437‑2734630

gCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑c                                                                                                >  1:1338658/1‑100 (MQ=255)

   gggCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGTACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑cccc                                                                                             <  1:496167/100‑1 (MQ=255)

   gggCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑cccc                                                                                             <  2:400434/100‑1 (MQ=255)

    ggCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑ccccc                                                                                            >  2:881567/1‑100 (MQ=255)

       aaGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑ccccccca                                                                                         <  2:2111341/100‑2 (MQ=255)

        aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑cccccccac                                                                                        >  2:333507/1‑98 (MQ=255)

        aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑cccccccac                                                                                        <  1:126178/100‑3 (MQ=255)

           tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACtt                                                                                     >  2:2113578/1‑100 (MQ=255)

           tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTCCACTTGCTTTACCGTCGACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTc                                                                                    <  2:1810713/100‑1 (MQ=255)

            gATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTc                                                                                    <  1:654739/99‑1 (MQ=255)

              tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGaa                                                                                 >  2:1069544/1‑100 (MQ=255)

              tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGaa                                                                                 >  2:664429/1‑100 (MQ=255)

               aaaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAAc                                                                                <  1:1845524/100‑1 (MQ=255)

                aaTATCCGCCAAATACTTATCACCAAGCTGAATGAAAGCTGCACTTGCTTTCGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCACCCCCCCTTCg                                                                                   >  1:1278043/1‑97 (MQ=255)

                  tatCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAAcgcg                                                                             <  1:31658/100‑1 (MQ=255)

                   atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAAcgcg                                                                             >  2:1444676/1‑100 (MQ=255)

                    tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAAcgcgcg                                                                           <  1:1435523/100‑1 (MQ=255)

                     ccGCCAGATCCTTATCACCATGCTGAATGAACGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAAc                                                                                >  1:1760906/1‑94 (MQ=255)

                     ccGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGt                                                                          >  2:2239811/1‑100 (MQ=255)

                       gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTg                                                                         <  1:1075868/100‑1 (MQ=255)

                       gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTg                                                                         >  1:880121/1‑99 (MQ=255)

                          aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGAcc                                                                     <  2:637375/100‑1 (MQ=255)

                           atctcTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCa                                                                    <  2:671207/96‑1 (MQ=255)

                             tACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCaa                                                                   >  2:165576/1‑100 (MQ=255)

                              aCTTATCACCATTCTGAATGAAAGCTTCACTTGCTTTAGCGTCTACTGCCATAACCGCCTCGCGAACCT‑‑CCCCCCCACGTCGAACGCGCGTGGACCAAtt                                                                 <  1:2438431/100‑1 (MQ=255)

                               cTTATCACCATGCTGAATGAAAGCTGCGCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCTCACTTCGAACGCGCGTGGACCAATTg                                                                <  2:1580105/100‑1 (MQ=255)

                               cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGCCCAAtt                                                                 >  1:2149460/1‑99 (MQ=255)

                               cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCaagat                                                                >  1:671204/1‑97 (MQ=255)

                               cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                >  1:2244844/1‑100 (MQ=255)

                                 tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGcc                                                              >  2:2343495/1‑100 (MQ=255)

                                 tATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                               >  1:983658/1‑99 (MQ=255)

                                  aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCg                                                             <  1:1208879/100‑1 (MQ=255)

                                  aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCg                                                             <  2:2120566/100‑1 (MQ=255)

                                  aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCg                                                             <  1:1208861/100‑1 (MQ=255)

                                    cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCgg                                                            <  1:980896/100‑1 (MQ=255)

                                    cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                           <  2:2475869/100‑1 (MQ=255)

                                    cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAAACT‑‑CCCCCACACTTCGAACGCGCGTGGACCAATTGCCgg                                                            >  2:109954/1‑99 (MQ=255)

                                    cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCGCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCgg                                                            <  2:1880385/100‑1 (MQ=255)

                                      ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGGCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGccc                                                         <  1:1392973/100‑1 (MQ=255)

                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGccc                                                         <  2:2370511/100‑1 (MQ=255)

                                        aTGCTGAATGAAAGCGGCCCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAGTTGCCGGcn                                                          <  1:332255/98‑2 (MQ=255)

                                          tcTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                      <  1:2462021/99‑1 (MQ=255)

                                          gCTGAATGAACGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                     <  1:1701414/100‑1 (MQ=255)

                                           cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    >  1:2445121/1‑100 (MQ=255)

                                           cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    >  1:1330276/1‑100 (MQ=255)

                                           cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    >  2:584330/1‑100 (MQ=255)

                                           cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCCAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    >  1:2087576/1‑100 (MQ=255)

                                           cTGAATGAAAGCTGCACTTGCTTTAGCGTCTAATGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTGCGAACGCGCGTGGACCAATTTCCGGCCCTaa                                                      >  2:2462038/1‑99 (MQ=255)

                                            tgaatgaaAGCTCCACTTGCTTTAACGTCTACTGCCATAAACGCATCGCGAACCT‑‑CCTCCCCACGTCGAACTCCCGTGCACCAATTCCCGGCCCTAAcc                                                    >  1:672401/1‑99 (MQ=255)

                                             gaatgaaAGCTGCACTTGCTTTAGCGTTTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                   >  2:284405/1‑100 (MQ=255)

                                             gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                                  <  2:2502521/100‑1 (MQ=255)

                                              aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                 >  1:145740/1‑100 (MQ=255)

                                                tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                               >  2:1674210/1‑100 (MQ=255)

                                                 gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                               >  2:2387309/1‑100 (MQ=255)

                                                  aaaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgag                                             <  2:628246/100‑1 (MQ=255)

                                                   aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                            >  2:2272391/1‑100 (MQ=255)

                                                   aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                            >  2:2307153/1‑100 (MQ=255)

                                                    aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                           <  2:688000/100‑1 (MQ=255)

                                                       tGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAatat                                        >  1:2053321/1‑100 (MQ=255)

                                                         cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                      >  1:1874124/1‑100 (MQ=255)

                                                         cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                      >  1:107350/1‑100 (MQ=255)

                                                            ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                   <  1:1046768/100‑1 (MQ=255)

                                                             tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCg                                  <  2:471987/100‑1 (MQ=255)

                                                               cTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                <  2:2523452/100‑1 (MQ=255)

                                                               cTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                <  1:463977/100‑1 (MQ=255)

                                                               cTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                >  2:1605181/1‑100 (MQ=255)

                                                                tttAGCGTCTACTGCCATAAACGCCTCGCGAATTT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAt                               <  1:929831/100‑1 (MQ=255)

                                                                   aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAt                            >  2:2057953/1‑100 (MQ=255)

                                                                   aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAt                            >  1:432369/1‑100 (MQ=255)

                                                                    gCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATc                           >  2:567287/1‑100 (MQ=255)

                                                                     cGTCTACTGCCATAAACGCCTCGCGACCCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                          <  1:548775/100‑1 (MQ=255)

                                                                     cGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                          <  2:1676853/100‑1 (MQ=255)

                                                                        cTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                       >  2:922985/1‑100 (MQ=255)

                                                                        cTACTGCCATAAACGCCTAGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAACAGGCGAGAATATTTTCCGGATTATCAGCa                       >  1:1262103/1‑100 (MQ=255)

                                                                         tACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa                      >  1:459791/1‑100 (MQ=255)

                                                                         tACTGACACACACGCGTCGCCAACCT‑‑CCCCCCCACTTCGAACGCGCGTAGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                       <  1:1118091/99‑1 (MQ=255)

                                                                          aCTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAGCGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                     <  1:1498518/100‑1 (MQ=255)

                                                                             gCCATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                  <  1:2358702/100‑1 (MQ=255)

                                                                              ccATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                 <  1:2503988/100‑1 (MQ=255)

                                                                              ccATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                 >  2:1063296/1‑100 (MQ=255)

                                                                              ccATAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                 <  1:1063292/100‑1 (MQ=255)

                                                                                aTAAACGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCCCGTGGACCAATCGCCGCCCCTACCCCGGCGCGAATACTTTCCGGCTTATCAGCAacac                   >  1:2285758/1‑96 (MQ=255)

                                                                                 tAAACGCCTCGCGAGCC‑‑GCCCCCCCACTTCGAACTCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa              >  1:1777406/1‑100 (MQ=255)

                                                                                  aaaCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa               >  1:1283273/1‑100 (MQ=255)

                                                                                  aaaCGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa             <  1:1726060/100‑1 (MQ=255)

                                                                                   aaCGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGACCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc            >  2:2379340/1‑100 (MQ=255)

                                                                                     cGCCTCGCGAACCT‑CCCCCCCCCCTTCGAACGCGCGTGGACCAACTGCCGGCCCTAACCAGGCGAGAATAATTTCCGGATTATCCGCAAAACAGGaaaac           >  1:1119442/1‑98 (MQ=255)

                                                                                     cGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAcc           >  2:254332/1‑100 (MQ=255)

                                                                                     cGCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg          >  2:648191/1‑100 (MQ=255)

                                                                                      gCCTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg          >  1:1001951/1‑99 (MQ=255)

                                                                                       ccTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc        <  2:880124/100‑1 (MQ=255)

                                                                                       ccTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc        >  1:2176449/1‑100 (MQ=255)

                                                                                        cTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc            >  1:2540576/1‑95 (MQ=255)

                                                                                        cTCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtct       >  2:1469655/1‑100 (MQ=255)

                                                                                         tCGCGAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc      <  2:454435/100‑1 (MQ=255)

                                                                                          cgcgAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc            >  1:895966/1‑93 (MQ=255)

                                                                                          cgcgAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct     <  1:2343481/100‑1 (MQ=255)

                                                                                          cgcgAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct     <  2:2049396/100‑1 (MQ=255)

                                                                                             gAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc  <  1:109953/100‑1 (MQ=255)

                                                                                             gAACCT‑‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc  <  1:389547/100‑1 (MQ=255)

GCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCC  >  NC_000913/2734437‑2734630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`