BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,186,032	G→C	5.6%	P63P (CCG→CCC) ‡	pepT →	peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,186,032	0	G	C	5.6%	366.2 / 9.2	107	P63P (CCG→CCC) ‡	pepT	peptidase T
Reads supporting (aligned to +/- strand): ref base G (54/47); new base C (3/3); total (57/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 109 total aligned reads displayed.

CTGGATACCCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATC  >  NC_000913/1185883‑1186181

ctGGATACCCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCg                                                                                                                                                       <  1:1337529/150‑1 (MQ=255)

      aCCCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAAc                                                                                                                                                 >  2:409208/1‑150 (MQ=255)

       cccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTaaag                                                                                                                                                >  1:899079/1‑148 (MQ=255)

       cccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  2:312910/1‑150 (MQ=255)

       cccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  1:1201128/1‑150 (MQ=255)

       cccAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACg                                                                                                                                                >  2:3622/1‑150 (MQ=255)

          aaTCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTAGCCGGTAACGTcg                                                                                                                                             <  1:519294/150‑2 (MQ=255)

          aaTCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTcc                                                                                                                                             >  1:839029/1‑150 (MQ=255)

            tCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTAGCCGGTAACGTcgcc                                                                                                                                           <  2:242652/150‑4 (MQ=255)

             cAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTAGCCGGTAACGTcgcca                                                                                                                                          <  1:380437/150‑5 (MQ=255)

                aaGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCg                                                                                                                                       >  2:285080/1‑150 (MQ=255)

                     ggggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATc                                                                                                                                  >  2:851312/1‑150 (MQ=255)

                       ggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATccc                                                                                                                                >  1:1182638/1‑150 (MQ=255)

                        gTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCg                                                                                                                               >  1:829748/1‑150 (MQ=255)

                        gTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCg                                                                                                                               <  1:675063/150‑1 (MQ=255)

                          gagaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGc                                                                                                                             >  2:637134/1‑150 (MQ=255)

                          gagaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATGGCATCTGCTGACAGCTCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGGGAGAAGGTCACTTTGATGTCGACGTTACCGGCTAACGTCCCTGTCGATATCCCGGc                                                                                                                             <  1:1263231/150‑1 (MQ=255)

                             aCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGAt                                                                                                                          >  1:1440229/1‑150 (MQ=255)

                                ggTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTgg                                                                                                                       >  1:637268/1‑150 (MQ=255)

                                 gTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGc                                                                                                                      >  2:1351025/1‑150 (MQ=255)

                                   tCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCtt                                                                                                                    <  1:779181/150‑1 (MQ=255)

                                   tCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCtt                                                                                                                    <  1:456563/150‑1 (MQ=255)

                                    cccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCttt                                                                                                                   >  1:541805/1‑150 (MQ=255)

                                    cccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCttt                                                                                                                   >  1:187204/1‑150 (MQ=255)

                                      cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTAt                                                                                                                 <  2:1312577/150‑1 (MQ=255)

                                       aGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTAtt                                                                                                                >  2:1292265/1‑150 (MQ=255)

                                        gCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTAttt                                                                                                               >  1:133888/1‑150 (MQ=255)

                                            ggAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCa                                                                                                           <  1:526330/150‑1 (MQ=255)

                                            ggAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCa                                                                                                           <  1:3622/150‑1 (MQ=255)

                                             gAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                          >  1:1275835/1‑150 (MQ=255)

                                             gAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                          <  1:358147/150‑1 (MQ=255)

                                               aGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATAAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgt                                                                                                        >  1:1242596/1‑150 (MQ=255)

                                                ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgtg                                                                                                       <  2:1365909/150‑1 (MQ=255)

                                                ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCAtgtg                                                                                                       >  2:598740/1‑150 (MQ=255)

                                                 gCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTgg                                                                                                      >  2:1117512/1‑150 (MQ=255)

                                                  ccAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                     <  1:899266/150‑1 (MQ=255)

                                                  ccAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                     <  1:893569/150‑1 (MQ=255)

                                                  ccAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGa                                                                                                     >  2:48190/1‑150 (MQ=255)

                                                      tGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAcc                                                                                                 <  1:637134/150‑1 (MQ=255)

                                                      tGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAcc                                                                                                 <  1:581245/150‑1 (MQ=255)

                                                        gAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTc                                                                                               >  2:546248/1‑150 (MQ=255)

                                                           gTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAcc                                                                                            >  1:349095/1‑150 (MQ=255)

                                                            ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCg                                                                                           <  2:53018/150‑1 (MQ=255)

                                                             tATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCgg                                                                                          >  2:805962/1‑150 (MQ=255)

                                                                 gCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTg                                                                                      >  2:1288713/1‑150 (MQ=255)

                                                                  cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                     >  2:681259/1‑150 (MQ=255)

                                                                  cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                     <  1:797734/150‑1 (MQ=255)

                                                                  cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                     >  2:75308/1‑150 (MQ=255)

                                                                  cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                     <  2:779990/150‑1 (MQ=255)

                                                                      tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTg                                                                                                                        >  2:426318/1‑111 (MQ=255)

                                                                      tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTg                                                                                                                        <  1:426318/111‑1 (MQ=255)

                                                                      tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCg                                                                                 >  1:1392789/1‑150 (MQ=255)

                                                                       gctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCgg                                                                                >  1:388588/1‑150 (MQ=255)

                                                                        ctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTtctc                                                                                                            <  2:966101/121‑1 (MQ=255)

                                                                        ctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTtctc                                                                                                            >  1:966101/1‑121 (MQ=255)

                                                                        ctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGc                                                                               >  2:1562636/1‑150 (MQ=255)

                                                                          gAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaa                                                                             >  1:1179442/1‑150 (MQ=255)

                                                                             agagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaaaaa                                                                          >  2:1522133/1‑150 (MQ=255)

                                                                              gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAt                                                                         <  2:305654/150‑1 (MQ=255)

                                                                                gcagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgt                                                                       >  1:459065/1‑150 (MQ=255)

                                                                                gcagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgt                                                                       >  1:41468/1‑150 (MQ=255)

                                                                                 cagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgtg                                                                      <  2:51067/150‑1 (MQ=255)

                                                                                      cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                 >  1:1370994/1‑150 (MQ=255)

                                                                                       gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                <  1:1016729/150‑1 (MQ=255)

                                                                                       gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                >  2:1236955/1‑150 (MQ=255)

                                                                                       gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCg                                                                <  2:133888/150‑1 (MQ=255)

                                                                                           agaTGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCaaa                                                            <  2:653151/150‑1 (MQ=255)

                                                                                               ggggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgt                                                        <  2:187204/150‑1 (MQ=255)

                                                                                                 ggCTTATAAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgttg                                                      <  2:1242596/150‑1 (MQ=255)

                                                                                                    ttATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTGTTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                   >  2:321373/1‑150 (MQ=255)

                                                                                                     tATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCACGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaaa                                                  >  2:284886/1‑150 (MQ=255)

                                                                                                        cAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTg                                                                                      <  2:194203/111‑1 (MQ=255)

                                                                                                        cAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTg                                                                                      >  1:194203/1‑111 (MQ=255)

                                                                                                           tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                            >  2:1300857/1‑150 (MQ=255)

                                                                                                           tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                            <  1:1351025/150‑1 (MQ=255)

                                                                                                            gtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgc                                           >  2:405726/1‑150 (MQ=255)

                                                                                                                ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                       >  1:1607378/1‑150 (MQ=255)

                                                                                                                ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                       <  1:1074192/150‑1 (MQ=255)

                                                                                                                ccTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTg                                       >  2:1452585/1‑150 (MQ=255)

                                                                                                                  ttAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGc                                     <  1:1300857/150‑1 (MQ=255)

                                                                                                                  ttAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGc                                     <  2:839029/150‑1 (MQ=255)

                                                                                                                   tAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCg                                    <  2:1274610/150‑1 (MQ=255)

                                                                                                                   tAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCg                                    <  2:1370994/150‑1 (MQ=255)

                                                                                                                        gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATAtt                               >  2:911879/1‑150 (MQ=255)

                                                                                                                          gaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgc                             >  1:1262541/1‑150 (MQ=255)

                                                                                                                          gaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTgc                             >  2:1043008/1‑150 (MQ=255)

                                                                                                                              ggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                         <  1:1224679/150‑1 (MQ=255)

                                                                                                                              ggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                         <  1:1117512/150‑1 (MQ=255)

                                                                                                                                cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTggg                       <  1:1521694/150‑1 (MQ=255)

                                                                                                                                cACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTggg                       <  1:312910/150‑1 (MQ=255)

                                                                                                                                       ccagggACGTTAGCCGGTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGc                >  1:104607/6‑150 (MQ=255)

                                                                                                                                       aTGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGc                <  2:1179442/150‑1 (MQ=255)

                                                                                                                                        tGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCg               <  2:459065/150‑1 (MQ=255)

                                                                                                                                        tGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCg               <  2:896606/150‑1 (MQ=255)

                                                                                                                                         agggACGTTAGCCGGTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGa              >  1:229146/4‑150 (MQ=255)

                                                                                                                                          gggACGTTAGCCGGTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGAt             >  1:159739/3‑150 (MQ=255)

                                                                                                                                               gTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGt        <  1:1522133/150‑1 (MQ=255)

                                                                                                                                                 tACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGttt      <  1:1292265/150‑1 (MQ=255)

                                                                                                                                                  aCCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt     <  2:829748/150‑1 (MQ=255)

                                                                                                                                                     ggCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATc  >  1:1183250/1‑150 (MQ=255)

CTGGATACCCAATCAAAAGCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATC  >  NC_000913/1185883‑1186181

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