BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,762	G→A	15.2%	A339V (GCC→GTC)	ypfG ←	DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,762	0	G	A	15.2%	342.8 / 42.4	143	A339V (GCC→GTC)	ypfG	DUF1176 family protein
Reads supporting (aligned to +/- strand): ref base G (65/56); new base A (10/12); total (75/68)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 149 total aligned reads displayed.

CCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAA  >  NC_000913/2581663‑2581860

ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACATCGTTGGCCAgg                                                                                                    <  1:1358363/100‑1 (MQ=255)

ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAgg                                                                                                    >  2:801777/1‑100 (MQ=255)

ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAgg                                                                                                    >  2:407059/1‑100 (MQ=255)

ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAgg                                                                                                    <  1:880112/100‑1 (MQ=255)

   gatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGt                                                                                                 <  2:111897/100‑1 (MQ=255)

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

     tgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCt                                                                                               >  1:1401403/1‑100 (MQ=255)

       cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCCGGCCTCTg                                                                                              >  1:93946/1‑99 (MQ=255)

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

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

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

       cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTgg                                                                                             <  2:10885/100‑1 (MQ=255)

         aCGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                                                           <  1:2303005/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:1437445/100‑1 (MQ=255)

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

              gACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCGACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                      >  1:2011002/1‑100 (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)

                     cTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTAt                                                                               >  1:2379071/1‑100 (MQ=255)

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

                      ttATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTtctc                                                                              >  2:1220197/1‑100 (MQ=255)

                         tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTCACGGGTGATCGACAAAGTGGCCGCGGCGTCTGGCTCATGTCAGTTaatccc                                                                           >  1:2195273/1‑95 (MQ=255)

                          ncGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                          >  1:96685/2‑100 (MQ=255)

                          ccGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                          >  2:1004249/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)

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

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

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

                                    gCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTg                                                                >  2:442410/1‑100 (MQ=255)

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

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

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

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

                                       tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCGGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                             <  1:1310983/100‑1 (MQ=255)

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

                                         ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcggcg                                                           <  1:1541035/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)

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

                                                    aCATGTCCCTACATTTTGTTTAACGGGTGATCCACAACGGGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                <  1:1056849/100‑1 (MQ=255)

                                                    aCATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                <  1:2498359/100‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                  tttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAc                                  <  2:15862/100‑1 (MQ=255)

                                                                   ttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAcc                                 <  2:508854/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:438935/100‑1 (MQ=255)

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

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

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

                                                                           cGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgc                         <  2:398079/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                  >  2:770860/1‑100 (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)

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

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

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

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

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

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

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

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

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

                                                                                           tGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAAt         >  2:567256/1‑100 (MQ=255)

                                                                                             ggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGc       <  2:265081/100‑1 (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:1261694/1‑100 (MQ=255)

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

                                                                                                  ggCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACaa  >  1:113839/1‑100 (MQ=255)

CCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAA  >  NC_000913/2581663‑2581860

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