BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,765	T→C	15.3%	D338G (GAC→GGC)	ypfG ←	DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,765	0	T	C	15.3%	316.6 / 47.3	143	D338G (GAC→GGC)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand): ref base T (67/54); new base C (10/12); total (77/66)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 147 total aligned reads displayed.

GATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGC  >  NC_000913/2581666‑2581864

 atgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTc                                                                                                    >  1:265078/1‑100 (MQ=255)

    cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATTCGGCCTACATGTCCCGCCAATTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                 <  1:1409594/100‑1 (MQ=255)

    cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                 >  2:950344/1‑100 (MQ=255)

    cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                 <  1:348866/100‑1 (MQ=255)

    cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                 >  1:2096994/1‑100 (MQ=255)

    cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                 <  1:1920409/100‑1 (MQ=255)

      aCGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                                                               <  1:2303005/100‑1 (MQ=255)

      aCGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                                                               <  2:1494916/100‑1 (MQ=255)

       cGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGccc                                                                                              <  1:1302082/100‑1 (MQ=255)

         cTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCAt                                                                                            >  2:840442/1‑100 (MQ=255)

         cTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCac                                                                                            <  2:384321/100‑2 (MQ=255)

         cTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCac                                                                                            <  2:384303/100‑2 (MQ=255)

          tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATg                                                                                           >  1:725386/1‑100 (MQ=255)

           gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgt                                                                                          <  2:1118858/100‑4 (MQ=255)

           gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgt                                                                                          <  2:1699783/100‑4 (MQ=255)

           gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgt                                                                                          <  2:77618/100‑4 (MQ=255)

           gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGAGCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                          >  1:2271091/1‑100 (MQ=255)

           gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCGACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                          >  1:2011002/1‑100 (MQ=255)

            aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                         <  1:950341/100‑1 (MQ=255)

            aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                         <  1:1832012/100‑1 (MQ=255)

            aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgtt                                                                                         <  1:114986/100‑5 (MQ=255)

                 tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgttgtgga                                                                                    <  1:9178/100‑10 (MQ=37)

                 tCTTATCCGGCATACAGCGTCTTATCCGGACTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTAGCCCATGCCAGt                                                                                      <  2:1678982/98‑1 (MQ=255)

                  cTTATCCGCCCTACAGCGTCTTATACGGCCTAAATGTCCCGCCATTTTGTTTAACGAGAGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                          >  1:941899/1‑93 (MQ=255)

                    tATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                 >  2:51438/1‑100 (MQ=255)

                       ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                              >  2:1004249/1‑100 (MQ=255)

                       ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                              >  1:99685/1‑100 (MQ=255)

                        cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAgg                                                                             >  1:2114229/1‑100 (MQ=255)

                         ggCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGt                                                                            <  2:503281/100‑1 (MQ=255)

                           ccTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                          >  2:553568/1‑100 (MQ=255)

                            cTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTggg                                                                         <  2:257639/100‑1 (MQ=255)

                             tACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTggg                                                                         <  1:1354082/99‑1 (MQ=255)

                             tACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGt                                                                        >  2:2370472/1‑100 (MQ=255)

                               cAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTc                                                                      >  1:2384768/1‑100 (MQ=255)

                               cAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTc                                                                      >  2:2309840/1‑100 (MQ=255)

                               cAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTc                                                                      <  1:1425719/100‑1 (MQ=255)

                                   gTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                  >  2:1922993/1‑100 (MQ=255)

                                   gTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgc                                                                   <  1:645148/99‑1 (MQ=255)

                                    tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                 <  1:1334154/100‑1 (MQ=255)

                                    tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                 <  2:1488045/100‑1 (MQ=255)

                                      ttATCCGGCCTACATGTCCCGCCATTTTTTTTTACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCTGTTATCGCAGGTGGGTTCTgcggcg                                                               <  1:1541006/100‑1 (MQ=255)

                                      ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcggcg                                                               <  2:1140085/100‑1 (MQ=255)

                                      ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcggcg                                                               <  2:1799148/100‑1 (MQ=255)

                                       tATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGt                                                              <  2:1451525/100‑1 (MQ=255)

                                            ggCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCGGCGGCGTAACGc                                                         <  2:1706559/100‑1 (MQ=255)

                                             gCCTACATGTCCCGCCATTTTGTTTAACGGGTGATGCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                            <  1:483698/96‑1 (MQ=255)

                                                tACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                     >  2:1045640/1‑100 (MQ=255)

                                                    tGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCg                                                 <  1:670728/100‑1 (MQ=255)

                                                    tGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCg                                                 >  2:495790/1‑100 (MQ=255)

                                                     gTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGa                                                >  2:1831141/1‑100 (MQ=255)

                                                      tCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                               <  2:2384783/100‑1 (MQ=255)

                                                      tCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                               <  1:159949/100‑1 (MQ=255)

                                                        ccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaaaa                                             >  1:1447785/1‑100 (MQ=255)

                                                        ccGCCATTTTGTTTAACGGCTGATCCAGAACGTGGGCCACGCGTCTGGTCGATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTCACGCACCCGGCGgaa                                              >  1:93046/1‑99 (MQ=255)

                                                           ccATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                          <  1:375969/100‑1 (MQ=255)

                                                            cATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                         <  1:1593794/100‑1 (MQ=255)

                                                                ttGTTTAACGGGTGATCAACAACGTAGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTACTGCAGCGTAACGAACCAGGCGAAAACGCTGAcc                                     >  1:466110/1‑100 (MQ=255)

                                                                 tGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCg                                    <  2:1610475/100‑1 (MQ=255)

                                                                   tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCgcc                                  <  1:968115/100‑3 (MQ=255)

                                                                   tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTc                                  <  2:130392/100‑1 (MQ=255)

                                                                   tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTc                                  <  2:438935/100‑1 (MQ=255)

                                                                    ttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCa                                 >  2:2125735/1‑100 (MQ=255)

                                                                       aCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                              <  2:949829/100‑1 (MQ=255)

                                                                         gggTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                            <  2:1139403/100‑1 (MQ=255)

                                                                          ggTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcgc                           >  1:575139/1‑100 (MQ=255)

                                                                           gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGcc                          >  1:537530/1‑100 (MQ=255)

                                                                             gATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAg                        >  2:122348/1‑100 (MQ=255)

                                                                              aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgc                       >  2:1986946/1‑100 (MQ=255)

                                                                              aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgc                       <  1:2452924/100‑1 (MQ=255)

                                                                               tCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcg                      >  1:720364/1‑100 (MQ=255)

                                                                                ccACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTATCGCACCGGGCGAAAACGCTGACCGTCAAAGCGCCAgcgc                     >  1:1820811/1‑100 (MQ=255)

                                                                                aactggcatGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcgc                     >  1:689940/9‑100 (MQ=38)

                                                                                 cacaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcgcg                    >  1:2098960/1‑100 (MQ=255)

                                                                                  ctggcatGGGCCAGACGCCTGGCCCATGCCAGTAATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAgcgcgc                   >  1:1710200/7‑100 (MQ=38)

                                                                                   caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                  <  1:1220193/100‑1 (MQ=255)

                                                                                   caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                  <  1:619359/100‑1 (MQ=255)

                                                                                   caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                  >  2:2314762/1‑100 (MQ=255)

                                                                                   caACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                  >  1:1123354/1‑100 (MQ=255)

                                                                                    ggcGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCt                 >  2:1094865/3‑100 (MQ=255)

                                                                                     gcatGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                >  2:801641/4‑100 (MQ=255)

                                                                                     gcatGGGCCAGACGCCTGGCCCACGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                >  2:392630/4‑100 (MQ=255)

                                                                                     aCGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                <  2:810574/100‑1 (MQ=255)

                                                                                      catGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGa               >  2:636147/3‑100 (MQ=255)

                                                                                      cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGa               >  2:416278/1‑100 (MQ=255)

                                                                                       atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa              >  2:1713942/2‑100 (MQ=255)

                                                                                       atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa              >  1:376690/2‑100 (MQ=255)

                                                                                       atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa              >  2:1463020/2‑100 (MQ=255)

                                                                                         gggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATg            >  2:108475/1‑100 (MQ=255)

                                                                                          ggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc           <  2:265081/100‑1 (MQ=255)

                                                                                           gCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGcc          >  1:1307714/1‑100 (MQ=255)

                                                                                           gCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAACGcc          >  1:2580299/1‑100 (MQ=255)

                                                                                            ccAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAATCGCCAGCGCGCCTGAATGCca         <  1:108474/100‑1 (MQ=255)

                                                                                            ccAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCCGCGCGCCTGAATGCca         >  2:2157773/1‑100 (MQ=255)

                                                                                            ccAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCca         >  2:193801/1‑100 (MQ=255)

                                                                                            ccAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCca         <  2:725389/100‑1 (MQ=255)

                                                                                             cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCcac        <  2:1990668/100‑1 (MQ=255)

                                                                                             cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCcac        >  1:1050520/1‑100 (MQ=255)

                                                                                                   tCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc  <  1:2132406/100‑1 (MQ=255)

                                                                                                   tCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc  <  1:1006175/100‑1 (MQ=255)

                                                                                                   tCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc  <  1:1006966/100‑1 (MQ=255)

GATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGC  >  NC_000913/2581666‑2581864

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